arXiv:1503.01263v3 [cond-mat.stat-mech] 18Jun2015 


Inverse scattering-theory approach to the exact large-n solutions of 0{n) models on 
films and semi-infinite systems bounded by free surfaces 


Sergei B. Rutkevich and H. W. Diehl 
Fakultat fur Phvsik, Universitdt Duisburq-Essen, D-A7058 Duisburq, Germany 

(Dated: June 19, 2015) 

The 0{n) (f>‘^ model on a strip bounded by a pair of planar free surfaces at separation L can be 
solved exactly in the large-n limit in terms of the eigenvalues and eigenfunctions of a self-consistent 
one-dimensional Schrodinger equation. The scaling limit of a continuum version of this model is 
considered. It is shown that the self-consistent potential can be eliminated in favor of scattering data 
by means of appropriately extended methods of inverse scattering theory. The scattering data (Jost 
function) associated with the self-consistent potential are determined for the L — oo semi-infinite 
case in the scaling regime for all values of the temperature scaling field t = {T — Tc)/Tc above and 
below the bulk critical temperature Tc- These results are used in conjunction with semiclassical and 
boundary-operator expansions and a trace formula to derive exact analytical results for a number 
of quantities such as two-point functions, universal amplitudes of two excess surface quantities, the 
universal amplitude difference associated with the thermal singularity of the surface free energy, and 
potential coefficients. The asymptotic behaviors of the scaled eigenenergies and eigenfunctions of 
the self-consistent Schrodinger equation as function of a; = are determined for x —>■ —oo. 

In addition, the asymptotic x —>■ —oo forms of the universal finite-size scaling functions 0(a:) and 
■i?(a:) of the residual free energy and the Casimir force are computed exactly to order l/x, including 
their 2 :“'^ln|a;| anomalies. 


I. INTRODUCTION 

In the vicinity of critical points of systems undergoing 
continuous phase transitions fluctuations occur on length 
scales ranging from microscopic separations up to the cor¬ 
relation length If the dimensionality of the system is 
below the upper critical dimension d* above which the 
Ginzburg-Levanyuk criterion mm holds arbitrarily close 
to the critical point, the fluctuations on all such scales 
affect the long-distance behavior in a nontrivial fashion 
[5]. The renormalization group (RG) has provided an 
appropriate conceptual framework for dealing with such 
problems involving many length scales and led to the de¬ 
velopment of powerful calculational tools for their quan¬ 
titative investigation [Ml- 

When such fluctuations in a near-critical medium are 
confined by external boundaries (walls) or macroscopic 
bodies immersed into it, effective forces are induced 
between the walls and these objects. The theory of 
such fluctuation-induced critical forces (critical “Gasimir 
forces” m) is substantially harder than the theory of 
bulk critical behavior because it involves a number of 
additional challenges. First and foremost, beyond bulk 
critical behavior, boundary and finite-size critical behav¬ 
iors must be treated in an adequate manner. Second, 
difficult dimensional crossovers are typically encountered, 
which perturbative RG approaches do not normally cap¬ 
ture [TB1125| . A particularly demanding case is that of 
three-dimensional systems whose large-scale physics can 
be represented by a (j)'^ model with an 0{n) symmetric 
Hamiltonian in a geometry that involves a crossover to 
two-dimensional behavior, such as a slab of size oo^ x L 
whose width L is finite. Provided the boundary condi¬ 
tions along the finite direction do not explicitly break 
the 0(n) symmetry, we know from the Mermin-Wagner 


theorem pSIES] that for finite L the system cannot ex¬ 
hibit long-range order at temperatures T > 0. Thus the 
low-temperature behavior strongly interferes with the di¬ 
mensional crossover. 

This combination of hard problems has hampered the 
design and application of satisfactory analytical theories. 
Exact solutions of appropriate models can — and have — 
provided helpful guidance and benchmarks for approxi¬ 
mations. An example is Danchev’s exact solution of the 
0(n) (j)'^ model on a cylinder of circumference L, i.e., a 
slab of thickness L subject to periodic boundary condi¬ 
tions uniisD]. Its incompatibility with the n-dependence 
obtained for the critical Casimir force by naive extrap¬ 
olation of the e-expansion results of [Ml HI] to d = 3 
dimensions [n in] strongly hinted at a breakdown of 
the e expansion [Ml Hi] ■ 

Unfortunately, real experimental systems of finite size 
usually involve free rather than periodic boundary condi¬ 
tions (pbc). In the present paper we shall be concerned 
with the exact n —)■ oo solution of the 0{n) (f)'^ theory on 
a slab of size oo^ x L subject to free boundary conditions 
along the finite direction (called z-direction henceforth). 
Unlike the case of pbc, where the n ^ oo limit leads to a 
translation-invariant constraint Gaussian model equiva¬ 
lent to the spherical model m, the breakdown of trans¬ 
lation invariance along the z-direction due to the free 
boundary conditions (fbc) entails that the exact n —> oo 
solution involves a self-consistent Schrodinger equation 
with a z-dependent potential V (z). Upon solving this 
self-consistency problem by numerical means, very pre¬ 
cise results for the temperature-dependent scaling func¬ 
tions of the excess surface free energy and the Gasimir 
force could be obtained in [32l [33] • 

The aim of this paper is to explore the potential of 
inverse scattering-theory methods [Ml [M] for obtaining 
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exact analytical information about the self-consistent po¬ 
tential V(z) and the above-mentioned scaling functions 
for temperatures T ^ both for the semi-infinite case 
L = oo and for finite L. 


To put things in perspective, recall that the available 
exact analytical knowledge about n —>■ oo solutions for 
fbc is rather scarce. Bray and Moore I3S1I37] managed 
to determine the self-consistent potential V{z) for the 
semi-infinite case precisely at the critical point T = T^. 
Building on this work, we computed in [35] the correc¬ 
tion to y(z; L = 00 , t = 0) linear in t oc (T — Tc)/Tc and 
combined it with results deduced via boundary-operator 
and operator-product expansions to work out a few other 
exact analytical properties. These were accurately con¬ 
firmed by the numerical results of (33] [33] , along with the 
exact analytical information about the low-temperature 
behavior of the Casimir force derived in the latter one of 
these two papers. 

The basic idea of our subsequently developed approach 
based on inverse scattering-theory methods is to elim¬ 
inate the potential V{z) in favor of scattering data. 
Since the self-consistency equation for V(z) corresponds 
to a stationarity condition for the free-energy density, 
one can exploit the latter to determine the scattering 
data from the corresponding variational equations. Pro¬ 
ceeding in this way in the semi-inhnite case enables us 


to get the scattering data for all temperatures t = 0. 
In this manner, the determination of V(z) can be by¬ 
passed, although V{z;L = oo,t) can in principle be re¬ 
constructed as the solution to an integral equation anal¬ 
ogous to those of Gelfand, Levitan, and Marchenko 
[331 [35] . From Bray and Moore’s solution [351 (SZ] it fol¬ 
lows that the potential becomes singular at the boundary 
planes [38] . Although singular potentials were considered 
in some inverse scattering-theory investigations [551ISH] j 
the particular kind of boundary singularities that the self- 
consistent potential v{z) exhibits at c? = 3 has not yet 
been investigated and requires appropriate modifications 
of the established inverse scatterin theory. The neces¬ 
sary extensions of the latter are described in a separate, 
accompanying paper [40] (referred to as II henceforth). 
Here the results required from II will simply be stated 
and applied. 

The remainder of this paper is organized as follows. 
Although our ultimate interest is in the n —>■ oo solution 
of the continuum (j)'^ model on a strip of size oo^ x L, 
we start in the next section with a discretized version of 
it, namely, the lattice model called “model B” in the 
numerical analyses of [331 [ 53 ] . The chosen discretization 
serves to avoid any ultraviolet (UV) (bulk and surface) 
singularities. We then recall the exact n —>■ oo solution of 
the model in terms of the eigenenergies and eigenstates of 
a self-consistent one-dimensional Schrodinger equation, 
discuss the simplihcations that can be achieved by tak¬ 
ing the limit g —>■ oo of the coupling constant and the 
addition of appropriate irrelevant interactions, and re¬ 
formulate the self-consistency equation via Green’s func¬ 
tions. In Sec. Ill we turn to the analysis of the continuum 


limit of the model and its self-consistency equation, recall 
known properties of the self-consistent scaling solutions 
for the potential and required background on the thermal 
singularities of the surface free energy and excess energy, 
their logarithmic anomalies, and an associated universal 
amplitude difference. In Sec. |IV[ we first provide the nec¬ 
essary background on inverse scattering theory, next use 
it to reformulate the self-consistency equations in terms 
of scattering data rather than the potential, and then de¬ 
termine the scattering data (scattering amplitudes and 
phase shifts) for temperatures t > 0 , t = 0 , and t < 0 . 

In Sec. [V] these data are exploited to obtain exact re¬ 
sults for the asymptotic large-distance scaling forms of 
the two-point order-parameter correlation function at, 
above, and below T^. In Sec.JV^the exact scattering data 
are used in conjunction with a trace formula, semiclassi- 
cal expansions, and perturbative asymptotic solutions of 
the Schrodinger equation to determine the limiting value 
vq of the regular part of the self-consistent potential at 
the boundaries, universal amplitudes associated with the 
thermal singularities of the surface excess energy, and the 
surface excess squared order parameter for t < 0. 

Sec. jVIIj deals with the asymptotic behavior of the 
universal scaling functions 0(a;) and 'd{x) of the in¬ 
dependent part of the surface free energy per boundary 
area and the critical Gasimir force in the low-temperature 
scaling limit x = t(L/^+)^/'^ —>■ —oo, where = 1 is 
the critical exponent of the bulk correlation length ^ 
and its nonuniversal amplitude for t —>■ 0-1-. This re¬ 
quires the determination of the x —)■ —oo behavior of the 
eigenenergies and the eigenfunctions of the self-consistent 
Schrodinger equation. Subsequently, the universal ampli¬ 
tudes of the logarithmic anomalies x~^ In |a:| and the x~^ 
terms are computed. The hrst agree with the results pre¬ 
viously derived from a nonlinear sigma model in [33] ; the 
latter is new. Section jVIII] contains a brief summary of 
our main results and concluding remarks. Finally, there 
are 7 appendixes describing technical details of some of 
the required calculations. 


II. THE LATTICE MODEL AND ITS n ^ oo 
LIMIT 


A. Definition of the lattice model 


We consider a lattice model on a three-dimensional 
slab with the 0(n)-symmetric Hamiltonian 




2 = 1 




, ( 2 . 1 ) 


where j = (ji,j 2 , ja) e with 1 < A < Ni, i = 1,2,3, 
labels the sites of a finite simple cubic (sc) lattice whose 
lattice constant we denote as a [41]. For the sake of 
brevity, we will write j\\ = (^ 1 ,^ 2 ), js = j, and N 3 = N 
henceforth. Each 4>j is an n-vector spin G K” of uncon¬ 
strained length, and the represent orthonormal unit 
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vectors pointing along the principal directions i of the 
lattice. 

Along the first two directions we choose pbc: 

= 4 >j for i = 1,2. (2.2a) 

Along the third one, we impose Dirichlet boundary con¬ 
ditions, requiring 


and residual free energy densities by 

/b(r) ^ hm /w(T)/lV, 

N^oo 

r^{T) = fN- NMT), 

MT) = r^{T)/2, 

r^%T) = r^{T) - 2MT), (2.5) 

respectively, and the Casimir force 


(pj.. j = 0 for j = 0, A^ -I- 1. 


(2.2b) 


In conjunction with the chosen interactions in "Hi, the 
prb (2.2a) imply that the model has translation invari¬ 


ance with respect to lattice translations along the di¬ 
rections i = 1,2. By contrast, the Dirichlet boundary 
conditions (2.2b) break the corresponding discrete trans¬ 


lation invariance along the direction (z = 3) normal to 
the boundary planes j = 0, N. 

The latter breakdown of translation invariance gener- 
ically occurs for models with free surfaces perpendicu¬ 
lar to the 3-direction. As a simple and natural gen¬ 
eralization one might want to consider analogs of our 
model (2.1) where the strength of the coupling between 


<pj and 4>ji on nearest-neighbor (NN) sites has been 
changed from the uniform value AT = 1 to different ones 
Ki and K 2 in the layers j = 1 and j = N, respectively 
(cf. the standard semi-infinite lattice models reviewed in 
[32 113] )• However, this is unnecessary. Since a phase 
with long-range surface order is ruled out in the bulk 
limit A^i—>00,1=1,2,3, for all temperatures T > 0 and 
arbitrary finite values of Ki/K and K 2 IK by extensions 
of the Mermin-Wagner theorem [26H28], only the ordi¬ 
nary surface transition [421144] remains for this general¬ 
ized three-dimensional model in the semi-infinite case. It 
is well established gMI] that the universal surface crit¬ 
ical behavior at this transition is described by a contin¬ 
uum field theory satisfying large-scale Dirichlet bound¬ 
ary conditions. Deviations from this boundary condition 
induced by modified surface interactions (which entail 
Robin boundary conditions for the continuum theory) are 
irrelevant in the RG sense. As is expounded in jSSj) this 
irrelevance can be checked explicitly in the large-n theory 
for a three-dimensional film of finite thickness L = Na. 
For the sake of simplicity, we will therefore restrict our¬ 
selves to the above-defined model with K = Ki = K 2 


and the boundary conditions (2.2). 


For later use, we introduce the partition function of 
the model 


/?.Fc(T,iV) = = -^fNiT). (2.6) 

Near the bulk critical temperature Tc, the asymptotic 
behaviors of the singular parts of the above A^-dependent 
functions on long length-scales are described by familiar 
scaling forms. Let us introduce a dimensionless tempera¬ 
ture variable t = const (T/Tc — 1) and fix the proportion¬ 
ality constant by absorbing in it the nonuniversal ampli¬ 
tude of the bulk correlation length [52] 

^(^)~e+(T/^c-l)-^ T>T,, (2.7) 

in the disordered phase, defining 

t = ( 2 . 8 ) 

and the scaling variable 

a; = (2.9) 

Following the conventions of [ini[311[33], we write the 
scaling forms of the residual free energy and the Casimir 
force as 


f^%T) - Six). 

(2.10) 

and 


PTciT, L)~N-‘^i?(x). 

(2.11) 

The scaling function z?(a:) is related to 0(a;) via (see, e.g., 

[I1I11IM1I31I23]) 

'&{x) = {d — 1) 0(x) — — Q'{x). 

(2.12) 

Furthermore, the value of the function 0(a;) at cc = 0, 
which measures the strength of the critical Casimir force, 
is the so-called Casimir amplitude 

o' 

© 

II 

0 

<1 

(2.13) 


Z = 


/ 




(2.3) 


and define the n —> 00 limit of the free energy per number 
of boundary sites and number of components by 


Jn = - lim 


lim 


InZ 


n-foo Ni,N 2 ^oo n Ni N 2 


(2.4) 


We also introduce the associated bulk, excess, surface. 


Since we have fixed the scale of x and no other (e.g., 
surface-related) macroscopic lengths are present, the 
functions 0 and d are universal [53] . 


B. Large-n self-consistency equations 

The techniques for deriving the equations that govern 
the n —> 00 limit of models such as ours are well estab¬ 
lished (see, e.g., dn], IM], |5S], Appendix B of [S5], and 
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[35]). For the case of the lattice model defined above, 
they have been explicitly given in |33j . Hence we can 
be brief and list them in the slightly different notation 
preferred here. 

Let us disregard for the moment the possibility that 
the symmetry is spontaneously broken in the bulk limit 
Ni ^ oo, i = 1,2,3, focusing on the disordered pha se. 
Then, in the limit n —>■ oo, the lattice model (2.11 is 


equivalent to n copies of a constrained Gaussian model 
for a one-component field <i>j with the Hamiltonian 


Upon taking the limits Ni,N 2 


oo, the self- 


consistency equation implied by Eqs. (2.14) and (2.15) 
becomes 


N 


^-^3= -I Mj)? = -I (jW2(H)|j) 


(2.16) 

where W 2 (A) is a particular one of the Watson integrals 


J E [ E(‘J’l+e. - + ^3^ - I iV, - rf 


2 = 1 


(2.14) 

Here and Vj is a self-consistent potential sat¬ 

isfying the constraint 


^3=r+l{^lJ=T+^-{^l^). 


(2.15) 


iu.(A).r^.../ 

Jo 3T Jq 


^dqd 


0 X + 4J2i=iSin\qi/2) 

(2.17) 

The index ly = 1,..., iV in Eq. (2.16) labels the eigenval¬ 


ues and orthonormalized eigenstates fi,{j) = {j\v) of 
the discrete Sturm-Liouville problem 


N 


E H,Mf) = vMj) + mAj) - Uj - 1) - Uj + 1)] = MAj), j = 1 ,..., iv. 

i'=i 


(2.18) 


The latter must satisfy the Dirichlet boundary conditions 

^(0) = f,(iV + l) = 0. (2.19) 

We choose them real-valued, so that their orthonormality 
relations become 


N 


AW) = '^fAj)W{j) = A 

3 = 3 


( 2 . 20 ) 


The coefficients Hj ji in Eq. (2.18) correspond to the el¬ 
ements of an X iV tridiagonal matrix “Hamiltonian”, 
namely 


/2 + Ui -1 \ 

-1 

-1 

V -1 2 + VnJ 


( 2 . 21 ) 


In order that the free energy Jn be well defined, H must 
be positive dehnite, i.e., we must have e,^ > 0 for all 
13=1,. ..,N. 

Note that the trace of H is simply related to the trace 
of the diagonal potential matrix V = (V) Sjji). One has 


trH = 2Af-htrV. (2.22) 


Further, by symmetry, the self-consistent potential Vj 
must be even under reflection about the midplane, i.e.. 


U, = Vn+1-j. (2.23) 


Straightforward evaluation of the free energy (2.4) 
gives 


N N 

/^ = 2 E ^AeA - f Y.(^3 - A" + In, (2-24) 

where is a trivial background term independent of 
Vj, which can be eliminated by a shift and is henceforth 
dropped. The function Ud{X) means the antiderivative 


= /' V ■ ■ i" ^ + 4 |:si„^(,./ 2 )| 

(2.25) 

of the Watson integral (2.17). The explicit form of its 
derivative for d = 2, WWAW in te rms of a complete el¬ 
liptic integral is given in Eq. (A1) of Appendix Inte¬ 
grating tins result yields the explicit form of U 2 (Xd given 
in Eq. (A3) (cf. Eq. (48) of [SB], and [H]). We shall not 
work with the explicit forms of these functions, but will 
make use of some of their properties below. We postpone 
a discussion of these properties for the time being. 

As has been shown in |33] . the free energy /jv (with 
omitted) is given by the global maximum of the 
functional /Ar[V] of the potential V = {Vj) defi ned b y 
the first two terms on the right-hand side of Eq. (2.24). 
That is 


/w = nmx/Af[V] (2.26) 
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with 

MV) = I tr{C/ 2 [H(V)]} - A tr[(V - f)^], (2.27) 


where we have explicitly indicated the dependence of the 
matrix Hamiltonian (2.21) on V. That the self-consistent 
potential corresponds to the global maximum of /Ar(V) 
is a consequence of two facts: (i) fwiV) is concave in 
V because it is a difference of a concave and a convex 
function of V [SD]; (ii) the self-consistency Eq. (2.16) is 
equivalent to the extremum condition 


dfNiV) 

dVj 


(2.28) 


Taking the bulk limit N ^ oo in Eq. (2.16) shows that 
the bulk critical value of f is given by 


Tc = -|bE3(0). (2.29) 

6 

Eollowing 1331, define a temperature variable 
t into which the amplitude = g/{ 2 ATT) of the bulk 
correlation length ^ ~ ^_|_/(f —fc) for deviations f — fc —)■ 
0-1- is absorbed by 


247r,„ 

t = - (t — Tc) = 247r 

9 


- +V3(0) 
9 6 


(2.30) 


Owing to this definition, the bulk critical point is located 
at t = 0. 


C. Simplifications of the equations 


Substitution of the solution V = V(t, N) of Eq. (2.32) 
for given t and N that maximizes the functional (2.31) 
gives us the g = oo analog /^^(t) = f^\v{t,N),t] of 
the free energy /^v in Eq. (2.26). Taking its temperature 
derivative yields the useful equation 


^^y^ = ^trV(f,7V). (2.33) 

Let us also list some known results |33j for bulk quanti¬ 
ties that will be needed below. The bulk value of the self- 
consistent potential, Vh{t) = limAr_>oo L)v(i,-^), maxi¬ 
mizes the bulk analog 




^C/3[Vb] + Eb 


t 

Stt 


kE3(0) 

2 


(2.34) 


of the function /(°°)[V,t], where we temporarily restrict 
ourselves to the disordered bulk phase t > 0. When t < 0 
and IV = oo, one must allow for spontaneous symmetry 
breaking. Since the corresponding bulk results can be 
found in textbooks such as [SS] and elsewhere (see, e.g., 
o [ini EH), there is no need to rederive known bulk re¬ 
sults for t < 0 here. Instead, we shall incorporate them 
into the t ^ 0 results for 14 and given in Eqs. (2.37) 
and (|2.38|) below. 


The necessary condition dfl^''[Vh,t]/dVh = 0 yields 


W3M)-W3iO) = - — , t> 0 . 

47r 


(2.35) 


The function IT 3 (A) is known to behave as [59] 


1. due to the limit g ^ oo 


WM) = w"3(o) 


1 

47r 


yX-f 0(A) 


(2.36) 


Eurther simplifications can be achieved by taking the 
limit g —>■ oo. In this limit, the model ( |2.1| ) goes over into 
a layered spherical model in which ^ ■ ($4 4/ {N 1 N 2 ) = 

(4)0j) fulfills a separate constraint for each layer j = 
1,..., TV. As we know from [321 [SS], taking the limit g —>■ 
00 causes a suppression of corrections to scaling, making 
it easier to extract the universal large-scale behavior. To 
take this limit, we add a contribution regular in t and 
define 


fl^\v,t)= lim 


g^oo 


/»(V) + JV|(A_„.,(0) 


2n 


= itrt/2[H(V)] + i[£-144(0) 


trV, 

(2.31) 


where we have explicitly indicated the t-dependence of 
the limiting function, and the superscript reminds us 
that g has been set to 00 . 

The g = 00 analog of the self-consistency condition 


±-W3(0) = -(j|4E2(H)|j). (2.32) 


2.16) follows from df}^\Y,t)/dVj = 0; it reads 


for small A > 0. Using this together with the fact that 
14 for t < 0 corresponds to the inverse transverse bulk 
susceptibility (which vanishes on the coexistence line), 
one concludes that 




(2.37) 


where 0(t) is the Heaviside step function. The familiar 

result [IS 13111331 EH] 




(2.38) 


for the leading thermal singularity of /b can be recovered 
by integrating Eq. (2.36) to obtain 173 (A) for small A and 
substituting Eq. (2.37) into Eq. (2.34). 


2 . due to an appropriate choice of U 2 


The above Eq. (2.31) for the free-energy function and 


the self-consistency condition (2.32) can be simplified fur¬ 


ther. This is because the functions 144 and U 2 are ex¬ 
pected to contain contributions that are irrelevant in the 
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sense that their omission does not change the asymptotic 
large-scale behavior such as the leading thermal singular¬ 
ities of the bulk and surface free energies, and the scaling 
functions 0(a;) and 'd{x) of the residual free energy and 
the Casimir force. Hence we should be able to replace 
W 2 and U 2 by appropriate simplified functions W 2 and 
U 2 that differ from W 2 and U 2 hy such irrelevant con¬ 
tributions. To justify the choice of W 2 and U 2 we are 
going to make below, we need some properties of their 
exact counterparts W 2 and U 2 , which are established in 

Appendix [X] and w ill now be discussed. _ 

From Eq. ( |2.17 ) or the closed-form expression 
one sees that W 2 (A) is analytic in the complex A-plane 
except for the branch cut [—8,0]. Arguments given in 
Appendix [A] show that it can be written for small |A| as 


W^ 2 (A) — —W 2 {X) In A -|- i? 2 (A), 


(2.39) 


where both ?i' 2 (A) and i? 2 (A) are analytical functions at 
small enough |A|. The former one is given by the spectral 
function 


u;2 (A) = [fF2 (A - iO) - 1^2 (A + iO)] 

ZTTl 

=-ImVF2(A -I- iO), 


(2.40) 


which for real A characterizes the singularity across the 
branch cut. Integration of Eq. (2.39) shows that U 2 can 
be written as 

;72(A) = ^t]^[l + AA(A)]+H(A) (2.41) 

Att 

with regular functions A(A) and B{X). 

Let us introduce the analogs 

^2(A) = l-i|^, (2.42) 


W2(A) = 


47r 


(2.43) 


and 

C/2(A) = A-h —A(l-lnA)-2 (2.44) 
47r 

of the functions W 2 , W 2 {X), and U 2 corresponding to the 
substitutions 


A(A) 

B{X) 


i(A) = 0, 
i?(A) = A-2, 


(2.45) 


in Eqs. (|2.39)-(2.41). Then the free energy func¬ 


tion (|2.31|) can be decomposed as 

/^~)(t,V) = /A,(t,V)+tri?(H) (2.46) 

into the contribution 

/Ar(t,V) = ^{tr[H(l-lnH)]+ttrV} (2.47) 

Stt 


associated with U 2 and the remainder 
i?(A) = —A(A)A^(1 - InA) -I-^ -I 


VF3(0) 


(2-A). 
(2.48) 

In the derivation of /Ar(t, V) we used Eq. (2.22) and an¬ 
ticipated that IT 3 (A), the analog of the function W 3 (A), 
has the property 


1^3(0) = 1. 


(2.49) 


The latter follows from the explicit expression given be¬ 
low in Eq. (2.52). 

The power series expansions of A and B can be deter¬ 
mined in a straightforward fashion. Explicit results to 
low powers of A can be found in Eqs. (A14) and (A15) of 
Appendix [A| 

The functions ^(A) and IT 3 (A) associated with C /2 (A) 
and W 2 can be obtained via the analog of the relation 

t/3(A) = y" ^C/2[A + 4sin2(p/2)]. (2.50) 


This gives 


C/3(A) =A + 


v/A(4 + A) 


2-h A 


47r 

In 


47r 2 +X+ y/X{A -t A) 


and 


(2.51) 


1E3(A) = C/'(A) = 1 + —In 


47r 2 -t A -f ^A(4-h A) 


. (2.52) 


Eurthermore, the stationarity condition for V corre¬ 
sponding to the omission of 7?(H) in Eq. (2.46), namely 
dfN{t,Y)/dVj = 0, takes the simple form 


(j|lnH|j)=t 

known from the analysis of model A in 


(2.53) 


D. Green’s function reformulation of the 
self-consistency equation 

Being interested in the universal large-length-scale 
properties of the solutions to the above equations, we 
consider their behavior in an appropriate continuum 
(scaling) limit a —> 0. To this end, their reformulation in 
terms of Green’s functions turns out to be helpful. 

Let us introduce the resolvent 

N 

G(A) = (A1 - H)-i = ^ G,.y(A)|j)(/| (2.54) 

ij'=i 

for A € C \ spec(H), where spec(H) = {ey,V^} denotes 
the spectrum of H and a self-explanatory Dirac notation 
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is used on the right-hand side. Equation (2.32) can be 
rewritten as 

±-W3iO) = -£ ^W^2(A)G,.,(A), (2.55) 

where Ci encircles the spectrum in a counterclockwise 
fashion as depicted in Fig. 


ImA 



• ReA 


FIG. 1. Integration paths Ci and C 2 in the complex A-plane. 


Since the function W 2 {X) is analytic in the complex 
A-plane except for the branch cut [—8, 0], we can deform 
the contour Ci into C 2 (see Fig. and recast Eq. (2.55) 
as 


± - tF3(0) = r dAn;2(A)G,-,(A), (2.56) 

J —8 


where W 2 is the spectral function defined in Eq. (2.40). 


It is convenient to rewrite the term W3(0) on the left- 
hand side also as an integral involving W 2 - Using the 
analog of Eq. (12.501) for W 3 gives [52] 


WM= j ^W2{2-2cosp)= 

W2iX) 


' —TT 

= f dX 


-8 v/A(A-4)' 


dA IU2(A) 

0 tt y^A(4- A) 

(2.57) 


This can now be subtracted from Eq. (2.56) to recast 
Eq. (2.56) as 


t 

Att 


= f dX W 2 (X) 
7-8 


Gjj{X) + 


1 


v/A(A - 4) 


(2.58) 


The second term inside the square brackets of 
Eq. (2.58) corresponds to minus the bulk critical analog 


of the first one. Thus it must be given by the diagonal 
element G^’^(A) of the critical bulk Green’s function 


G‘’’^(A) = G(A)|v=o,^=oo. 


(2.59) 


To verify this, note that G for given V and finite N is the 
solution to (A1 — H) • G(A) = 1 subject to the Dirichlet 
boundary conditions 


Goj'iX) — Gpf+ij'iX) — 0 


(2.60) 


implied by Eq. (2.19). Hence G*’’‘’(A) is the solution to 


the difference equation 

(A-2)G^^,;;=,(A) + G^^;^,^.,(A) 

-yb,( 


+ G°’^,^^,iX)= 6 ,,p, (2.61) 

that vanishes as |7 — / | —>■ 00 , namely 
1 




exp[-g(A)|j-/|] (2.62) 


with 


q{X) = 2 arcsinh sJ—XjA. 


(2.63) 


The resu lt proves the above statement made below 
Eq. ( 2.58| ) about its diagonal element . 

Note that all equations given in this subsection re¬ 
main valid upon replacement of the functions W 2 and 
W 3 by their analogs with a tilde, except that the lower 
integration limit must be changed from —8 to —00 in 
Eqs. (jr^-lj^). 


III. SCALING LIMIT 

A. Scaling limit of the self-consistency equation 

We now turn to the study of the appropriate contin¬ 
uum scaling limit a —> 0 of the n = 00 self-consistency 


problem. Equation (2.58) is a convenient starting point 


since the subtracted bulk term eliminates the UV sin¬ 
gularities. To this end, we introduce the dimensionful 
quantities 


Vi = ]ia 
L = Na, 


2/2 = J2a, 


E = Xa 


-2 


z = ja, 
m = ta~ 


(3.1) 


and 


Va{z) = V,a-\ Qa{z,z'-E) = G,-y(A)a. (3.2) 

Note that the variable m, which reduces for t > 0 to the 
inverse of the bulk correlation length becomes negative 
for f < 0. 


Equation (2.58) changes into 

m = Ja{z), 

with 

nO 

Ja{z)=ATT dEw2{Ea'^) 

j- 8 la^ 


(3.3) 


Qa{z,z-,E) 


+ 


y/E{Ea^-4) 


(3.4) 


when expressed in these variables. We are interested in 
the limit a —> 0 at fixed L and m. Since the widths of the 
nonuniversal boundary regions \z\ < a and |L — z| < a 














































shrink to zero in this limit, we must be prepared to en¬ 
counter singular behavior at the boundary planes and 
UV singular contributions with support on the bound¬ 
aries. This means that Ja{z) should b e re garded as a 
distribution and the a —?> 0 limit of Eq. (3.31 interpreted 
accordingly. 

Before we turn to this issue, let us first discuss the 
continuum limits of Va{z) and the Green’s function Qa- 
Dimensional considerations imply that Va{z) can involve, 
on the one hand, o-independent terms that diverge as 
^ z~'^ and ^ (L — as 0 —>■ O-b and 2 —?> L—, respec¬ 
tively, and boundary singularities of the form a~^6'{z) 
and a~^5'{z — L), on the other hand. Let us disregard 
possible UV boundary singularities of the latter kind oc S' 
for the moment by restricting ourselves to the open in¬ 
terval (0,L) when studying the a —> 0 limit of Va{z) and 
assuming that the limit 


lim Vaiz) = v{z), 

a —>-0 


ze (o,L), 


(3.5) 


exists. That is to say, we consider v{z) as a function 
rather than a distribution. For the Green’s function Qa 
such caution is not necessary. Owing to the Dirichlet 
boundary condition (2.60) and the fact that its engineer¬ 
ing dimension is a, it should have a smooth a —>■ 0 limit 
for z and z' in the entire closed interval [0, L], so that the 
limit 

limga(z,z';E)=g(z,z';E), z, z'e[0,L], (3.6) 

a —>-0 

should exist. 

The differential equation that the limiting function 
g satisfies follows in a straightforward fashion. Noting 
that the matrix Hamiltonian (2.21), scaled as H/a^, ap¬ 
proaches the operator 


as a 


n. = -dl 

0 , we arrive at 


;(z) 


[E + dl - v{z)\g{z, z'; E) = 5{z - z'). 


(3.7) 


(3.8) 


We shall also need the c ontinuum analogs of the or¬ 
thonormality relation (2.20) and the spectral representa¬ 


tion (2.55) of the Green’s function. They read 


and 


(Mfi.') = f dzfl{z)f^>{z) = S^y 
Jo 

fAz) ft(z') 


g{z,z';E) = Y, 




E-E, 


(3.9) 


(3.10) 


where the asterisk on f* may be dropped since we will 
again choose real-valued eigenfunctions. 

If v{z) were not singular at the boundaries, we 
would find that both Eq. (3.8) and the corresponding 
Schrodinger equation 

(3.11) 


on the interval [0, L\ must be subjected to standard 
Dirichlet boundary conditions. The latter boundary 
conditions would then ensure the self-adjointness of the 
Hamiltonian Jiy. However, it can easily be seen that the 
self-consistent potential v{z) = v[z\ L, m) must be singu¬ 
lar at the boundary planes z = 0 and z = L. For d = 3, 
it must vary asymptotically as 


v{z-,L,m) ~ 


4 (L-z) 


-2 


z —y L — . 


(3.12) 


An easy way of seeing this is to recall Bray and Moore’s 
exact solution [Ml EH 

u,.(z;0)= 2<d<4, (3.13) 

for the semi-infinite critical d-dimensional case L = oo 
and m = 0. Here we have introduced the notation 


Uoo(z; m) = v{z; oo, m). 


(3.14) 


The solution (|3.13 ) refers to the so-called ordinary sur¬ 
face transition [42H44j . When 3 < d < 4, a second self- 
consistent solution with the same z-dependence oc z“^ 
but a different amplitude exists [MIEZ] , which pertains to 
the so-called special surface transition. Since we will ex¬ 
clusively be concerned with the (d = 3)-dimensional case, 
we shall not consider the latter transition and restrict 
ourselves to the ordinary one. 

Glearly, for general m ^ 0 and L < oo, the behavior of 
this potential v{z',L,m) at distances from the boundary 
planes z = 0 and z = L much smaller than l/|m.| and L 
must comply with the exact solution (3.13). This dictates 


the singular short-distance behavior (3.12). 

The result can also be understood within the frame¬ 
work of the boundary operator expansion (BOE) 

HHj. Since the application of the BOE to v{z) has been 
discussed in some detail in a recent paper [38], we can 
be brief. Its central idea is that scaling operators 0{y, z) 
with scaling dimensions A[C>] can be expanded for small 
distances z from the boundary plane z = 0 as 


=nE 


^a(=)-A[0 ] CoAmz, z/L) dj{y), (3.15) 


where Oj are surface operators with scaling dimensions 
The potential v{z) corresponds to the energy den¬ 
sity operator whose scaling dimension is {1 — a)/v = 
d — 1/iz, which is 2 for 2 < d < 4 when n = oo. The 


asymptotic behavior (3.12) results from the contribution 


of the unity operator 1. 

A similar reasoning can be used to clarify which bound¬ 
ary conditions must be imposed on the eigenfunctions 
fu{z). The leading contribution to the BOE of the order 
parameter <I>(y,z) ~ ^{y) originates from the boundary 
operator <&. Since the associated scaling dimensions of 
these two operators are (d — 2 -|- ? 7)/2 and (d — 2 -|- ryu )/2, 


'Hy'ipi^z) = Etp{z) 
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respectively, one has $(y,z) ~ for z —>• 0, 

where t] = 0 and r]\\ = 1 when n = oo and 2 < d < 4. 
Consequently the bo undar y conditions subject to which 
the Schrodinger Eq. (3.111 must be solved are 


oc 



for z —)• 0+, 
for z —>■ L — . 


(3.16) 


Again, this conclusion is in complete accord with the 
findings of Bray and Moore [MlEZl for the semi-infinite 
critical case. When L ^ oo, the spectrum spec('H„) be¬ 
comes dense and continuous. At to = 0, it is given by 
the interval [0, oo). The corresponding (improper) eigen¬ 
values Efe and eigenfunctions fkiz) can be labeled by a 
nonnegative variable k, rather than by a discrete index 
ly. They read [36l [37] 

h{z) = Vz jQ{kz), Ek = k'^, (T = oo,to = 0 ). 

(3.17) 

The corresponding Green’s function Goo,c{z,z';E) = 
G{z, z'; E;L = oo, to = 0) for A < 0 may be read off from 
Bray and Moore’s result for the pair correlation function. 
One has 


the integral. The result for the semi-infinite critical case 
just mentioned implies that the second boundary plane 
contributes a term —6{z — L)/2. Power counting rules 
out other contributions localized on the boundary planes. 
Hence, to complete the proof of Eq. ( |3.20 |, it remains to 
show that Ja(z) converges to Jsm{z) for all z G (0, L). 
The leading contribution to the integral on the right- 
hand side of Eq. (3.4) for small a comes from the region of 
small values of Ea"‘. Therefore, we can insert the small-A 
expansion of 


»,(A) = - 


A 5A^ 

~ 8 ^ 


0(A3) 


(3.22) 


of the spectral function which one finds from its explicit 


form given in Eq. (AlO) to conclude that the integral 
converges to Jsm as a —)■ 0 provided z G (0, L). 

The upshot is that the self-consistency equation be¬ 
comes 


m + ^[6{z) + 6{L — z)] 


dE 


G{z, z; E) -I- 


2 ^/^_ 


(3.23) 


Io{V^z<)Ko{V^z>), 

(3.18) 

where z< and z> denote the smaller and larger values of 
z and z', respectively. 

Note that for a given potential v{z) which exhibits the 
surface singular behavior specified in Eq. (3.12), the ma¬ 
trix elements 


ifl'Hvg) = / f*{z)ny g{z)dz (3.19) 

Jo 

of the Hamiltonian are well-defined for all complex¬ 
valued functions / and g belonging to the subspace of 
the Hilbert space T^([0, L]) of squar e inte grable functions 
satisfying the boundary condition (3.16). Furthermore, 


we have (flHvg) = {Hy/lg), so that Hy is a symmetric 
operator on this subspace. 

We are now ready to return to the continuum limit 
of the self-consistency Eq. (2.58), (3.3), (3.4). We claim 
that Ja{z) with a —>■ 0 is a representation of the distri¬ 
bution 


J{z) = Jsm(z:) - ]^[5{z) + 5{z - L)], 


whose “smooth part” Jsm('2^) is given by 

1 


ds' 


aW = / 

J — ( 


dE 


G{z,z-,E) + 


2v^, 


(3.20) 


0 < z < L. 


(3.21) 

As has been shown in [SH], Eq. (3.20) simplifies to J(z) = 
—6{z)/2 in the critical semi-infinite case (see Eq. (3.8) 
of [3H|). That JsmA > 0;L = oo,TO = 0) = 0 i s known 
from [Sni [37]. The consis tency with Eq. ( |3.20 ) can be 
verified by inserting Eq. (3.18) into it and computing 


in the a —>■ 0 continuum limit considered. 

Equations (3.8) and (3.11), in conjunction with the 
boundary conditions (3.16), define a Sturm-Li ouvill e 
problem for potentials with the singular behavior (3.12). 
In order to determine the exact n —>■ oo solution of 
the model in the continuum scaling limit specified at 
the beginning of this subse ction, it must be solved self- 
consistently with Eq. (3.23). Finding closed-form analyt¬ 
ical solutions to these equations for noncritical tempera¬ 
tures TO ^ 0 and finite thicknesses L is a major challenge, 
which may well turn out to be too difficult to master. In 
[331133] this problem was bypassed by resorting to numer¬ 
ical solutions of discretized equations. This has yielded 
precise results for the universal scaling functions 0(a;) 
and d{x). 

Regrettably, not much exact information on the self- 
consistent pote ntial v{z;L,m) is available beyond the 
exact solution (3.13) for the critical semi-infinite case 
L = oo, TO = 0. To our knowledge, it is limited to what 
has been provided by us in a recent paper [33] . Before we 
turn to the issue of how the situation can be improved 
through the use of inverse scattering-theory methods, it 
will be helpful to recall these exactly known properties. 
This is done in the next subsection, where we also express 
quantities such as the excess energy density in terms of 
v{z; oo, to). 


B. The self-consistent potential 

The self-consistent potential v{z) is proportional to the 
energy density (<i>^(y,z)), where $(y, z) means the con¬ 
tinuum analog of the lattice field <I>j|| ^ derived from it via 
the scaling z). Therefore, v{z;L,m) 



































10 


must exhibit the corresponding scaling behavior 
v{z-L,m) = 

= z-^Tf{\m\z,z/L) (3.24) 

n=oo, a—6 

on long length scales [521IM] where ± refers as usual to 
m = t/a ^ 0. Although our ultimate interest is in the 
case d = 3, we give here — and where appropriate below 
— both results for general values of n and d G (2,4), 
and for the case of n = oo and d = 3. The latter case 
is special; it involves degeneracies, which imply logarith¬ 
mic anomalies in the leading thermal singularities of the 
surface free energy fs and related quantities [38|. The 
results for general d we are going to present here serve to 
illustrate the special features of the d = 3 case. 

As is discussed in [31] and elsewhere (see, e.g., [3311331 
1381155ll67| l. useful information about properties of the 
scaling function can be obtained from the BOE. For 
the case of the ordinary transition we are concerned with 
here, the leading contributions to the BOE arise from the 
unity operator 1 and the zz-component T^z of the stress 
tensor on the boundary. One obtains 

Xi")(z, 3 )+z"ri")(z, 3 ) + ...' 

Xf)(z,3) + z3yp)(z,3) + ...] , 

(3.25) 

where z and 3 denote the dimensionless distances 

z = \t\’'z/a = \m\z, 3 = z/L, (3.26) 

d—3 

while A™ is a normalization factor ensuring 

X2^)(0,0) = l. (3.27) 

Since BOE expansion coefficients such as and 
are short-distance properties, they are expected to be 
regular in the temperature variable. This implies 

X^^\z,0) = l±ai(d)zi/" + a2(d)z2/" + ... 

z—>-0 

= 1-1- ai(3) mz + 02 ( 3 ) m^z^ -I- ..., 

d—3, n—oo 

(3.28) 

where the value 

ai(3) = -^ (3.29) 

TT^ 

is known from [35] . The coefficients aj{d), j = 1,2,..., 
are independent of the sign of m because they can be 
expressed in terms of derivatives of the energy density at 
m = 0 and L = 00 . The coefficient 02 ( 3 ) is determined 
exactly in Sec. [VI A| We shall prove there that it has the 
value 


Ti")(z,3) =A°/ 

z—^0 


independent of the sign of m, where ^(a;) is the Riemann 
zeta function. 

Returning to the case of general d G (0,4) and n, con¬ 
sider the analogous expansion of y'^‘*^(z, 0 ), 

yi‘'^(z,0) = 6o,±(d)±6i.±(d)zi/" + ... 

z—>-0 

— ^o,± (3) + bi ± (3) mz + .... (3.31) 

d—3, n—oc 


As indicated, the expansion coefficients here must be ex¬ 
pected to have different values for m ^ 0. In fact, the 
term oc 60 ,± yields the leading singularity ^ |mp““ of 
the surface energy density and the ratio bo^+{d)/bo-{d) 
should take the same universal value as its analog for the 
bulk free energy [331 EHl US]- In resu lts fo r d = 3 and 
n = 00 given in the second line of Eq. (3.31) we have uti¬ 
lized the fact that 5o,± (d) does not have a pole at d = 3 
(see the discussion in [35]1. 

As has been discussed in [3S|, the above BOE expan¬ 
sion can also be applied to the L < 00 case at Tc to gain 
information about the distant-wall correction ^ {z/LY 
to the potential. The associated amplitude is propor¬ 
tional to the Casimir amplitude Ac. The proportional¬ 
ity factor is known from |38| for the {d = 3)-dimensional 
case. This led to the prediction 


u(z;L, 0 ) _ ^ 1024 Ac z^ 

u(z;OO,0)d=3 TT 


o[(z/L) 3], 


(3.32) 


which is in conformity with the numerical solution of the 
self-consistency equation [33] . 

Focusing on the case of d = 3, we can expand the 
scaling function T^^(z, 3 ) about 3 = 0 and match with 
the above equations to conclude that 


v{z]L,m) = — 

z/L^O Z^ 


A±{\m\z) + R±(|m|z) -b ... 


= u(z;oo,m)-b ^ R±(|m|z)-b ..., (3.33) 


where the functions A± (z) and B± (z) must vary for small 
z as 




56C(3) , 

TT* 


o{zY 


and 


B±{z) 


z^O 


256 Ac 

TT 


OY), 


(3.34) 


(3.35) 


respectively. In the large-z limit, A±{z) and B±{z) must 
decay on the scale of 1 /|to|. Depending on whether 
m > 0 or m < 0 , the latter length corresponds to the 
correlation length or Josephson length and the decay of 
these functions is exponential or algebraic, where the al¬ 
gebraic decay is due to the presence of Goldstone modes 
in the ordered bulk phase (m < 0,L = 00 ). In fact, re¬ 
sults derived in Sec. |VIlA| yield the asymptotic behaviors 


02(3) 


-^C(3) 


(3.30) 


A_(z) = ^-^+oil/zY (3.36) 

z —>^00 Zz z^ 
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and 

B_{z) = zM_^+o(1/z2) (3.37) 

z—yoo Z 

with 


/32 = C(3)/2 (3.38) 


and the unknown coefficient a 2 whose value we have not 
determined and shall not need below. 

For use below, let us also briefly mention the behav¬ 
ior of v{z;oo,m) for large z. In this limit, v{z;oo,m) 
approaches the bulk value Vhim), 

v{z;oo,m) = v\j{m) = rn? 9{m), (3.39) 

Z—^OO 


given by the square of the inverse correlation length when 
m > 0, and zero in the ordered bulk phase on the coex¬ 
istence curve. The approach to the limiting value Vh(m) 
is exponential ~ algebraic ~ (|to|z)“^, depend¬ 

ing on whether to > 0 or m < 0. Evidently, Vh(m) is 
the continuum limit of the scaled bulk lattice potential 
a~^Vb(t = ma), where 


Vh{t) = lim V)(A^ = oo,t) = lim —trV(A^,t) (3.40) 

j—>-oo N—>-oo jV 


is the bulk potential given in Eq. (2.37). 


C. Excess surface energy and its leading thermal 
singularity 


In Section fill Al we have seen that the continuum limit 


of the self-consistent equation, Eq. (3.23), is well-defined 


and free of UV singularities. However, quantities such as 
the excess energy density, and bulk and free energies still 
involve UV singularities, which must be properly sub¬ 
tracted to gain the desired information about universal 
quantities such as universal amplitude combinations and 
scaling functions. We first consider this problem for the 
excess surface energy density. To this end, we introduce 
the N X N matrix 


V{N,t) = V{N,t)-Vb{t)l, 
in terms of which the excess energy density 


(3.41) 


^s(i) = 


(3.42) 


of our lattice model can be expressed with the aid of 
Eq. (2.33) as 


SJt) = ^ lim trV(iV,<). (3.43) 

Stt N^oo 

For general dimension d G (2,4) with d ^ 3, the leading 
thermal singularity of the surface free energy /g for t -G 
±0 is of the form 


/s - - A^i\d) 


(3.44) 


where = J2k=o is a regular background 

term. As has been explained in [38], the degeneracy 
of the singular and regular terms oc at d = 3 entail 
by a standard mechanism m possible logarithmic tem¬ 
perature anomalies. The amplitudes ^^^(d) and /^’^^(d) 
generically are expected to have Laurent expansions of 
the forms 

(d) = - ^ + A+ 0(d - 3) (3.45) 

and 

A\d) = ^ + f^^l+0{d-3). (3.46) 

These yield the limit 


lim 
d —^3 - 


A±\d)\t\'^ ^ + f2\d)t'^ 


= t 




(3.47) 

which in turn implies the following small-t expansion for 
the surface free energy fs (t) in the three-dimensional sys¬ 
tem d = 3 near the bulk critical point, 

Mt)l=3 = [/o'^(3) + t + + 0(t^) 

+ Ag®^ -In |t| + 0{t^\n\t\). (3.48) 

Accordingly, one obtains for the singular part of the ex¬ 
cess energy density (3.42) 




(3.49) 


where 

AA(^) = A^^V - a[,^L (3.50) 

is a universal amplitude difference [71] . The residue £i 
is also universal and according to [38] given by 


/2,-l 


1 

47r3 


(3.51) 


Upon studying the continuum scaling limit of £s, we 
will now show that this amplitude difference can be ex¬ 
pressed as 


AA 


(s) _ tt-l- + Oi- 


0 


16tt 


(3.52) 


in terms of the following finite integrals involving the self- 
consistent potential Uoo(z;to) of the continuum theory: 

Uoo(z;±l) - Ub(±l) + A 



T — d(l-z) 

TT^Z 


(3.53) 


These integrals are finite because the two last terms in 
the square brackets subtract the boundary singularities 
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near z = 0, while the subtraction produced by the second 
term along with the restriction of the last term to the 
interval (0,1) implied by the theta function ensures the 
integrability at the upper integration limit z = oo. Note 


that the coefficient of the subtracted 


term has been 


chosen in accordance with the result for ai(3) given in 
Eq. (|^|. 

(s') 

To prove the asserted representation of AAq in terms 
ofvoo, we start from Eq. (3.43), express Vj{oo,t) in terms 
of its continuum analog Voo{z;m) = m'^ Voodmlz-jiLl), 
and use the Euler-McLaurin formula 


a'^f{ja)= j dz/(2:) + “/(a)+ 0(0^) (3.54) 

j=i 

to find 

f°° l~ 

in Ss{t) = a\m\ / z;oo(z;±) dz + - Fi(oo,t) + ... 

J a\m\ ^ 

(3.55) 

The second term is regular in t and does not contribute 
to From the integral we split off the contribution 


dz 


2|m| 


1 4 

4z2 TT^Z ^ ' 


4a|m 


1 4 

± ^ ln(a|TO|). 


that diverges as a —)■ 0, obtaining 


(3.56) 


Stt^s = a|m| / dz 

J a\m\ 
4a|TO' 


{ioo(z|±) + T ^ 6»(1 - z) 
4zZ TT^Z 


ln(a|TO|) + reg, 


(3.57) 


where “reg” represents contributions regular in t. Hence 
the leading thermal singularity of Ss is given by 




am 


StT 


a± =F ^ ln(a|TO|) 

TT^ 


(3.58) 


result for given in Eq. (3.52). 


Comparing with Eqs. (3.47) and (3.49) then yields the 


To determine the universal number AHq^^ in an ana¬ 
lytic manner from this equation, one would have to know 
the self-consistent potential Uoo(z;±l) for all values of 
z G (0, oo). Unfortunately, the latter is known only in 
numerical, but not in closed analytical, form. The in¬ 
verse scattering-theory techniques developed in Sec. |IV| in 
conjunction with the semiclassical expansions described 
in Appendix 0 will enable us to obtain an exact ana¬ 
lytic result for a_ and to express the quantities a+ and 
AAq in terms of a single, numerically computable inte¬ 
gral. The results are given in Eqs. (6.20)-(6.22) below 
[see Sec. VIB . 


IV. INVERSE SCATTERING THEORY 

A. Preliminaries 


The aim of inverse scattering theory ESJ [72] is to 
reconstruct the potential of the stationary Schrodinger 
equation v{r) from scattering data. Usually, either 
Schrodinger equations for three-dimensional systems 
with radially symmetric potentials or one-dimensional 
Schrodinger problems are considered. In the first case, 
one must deal with a radial Schrodinger equation in¬ 
volving an effective potential that differs from v(r) by 
a centrifugal term l{l + l)r“^, where I = 0,1,2 ,... are 
orbital angular momentum quantum numbers. Provided 
the effective potential satisfies certain conditions, such as 
absolute integrability of V{r) over intervals (5, oo) with 
6 > 0 and integrability of r|U(r) | over integrals (0, a) with 
a > 0, various classes of inverse scattering problems are 
known to have a unique solutions so that the potential 
can be determined — at least, in principle — from scat¬ 
tering data as solutions to certain integral equations such 
as the Gel’fand-Levitan or Marchenko integral equations 

ESj. 


As we have seen above, the Schrodinger equation we 
are concerned with here involves a self-consistent poten¬ 
tial v(z;L,m) that becomes singular at the boundary 
planes and has a leading near -boundary singularity of 
the form specified in Eq. (3.12). At L = oo and m = 0, 
it corresponds formally to a radial Schrodinger equation 
with angular momentum quantum number I = —1/2. Al¬ 
though there exists some recent work on inverse scatter¬ 
ing problems involving singular potentials ^ Ac/r^ with 
K > —1/4 [32], the case of k = —1/4, corresponding 
to the marginal value below which the particle is sup¬ 
posed to fall into the center according to Landau and 
Lifshitz m, requires appropriate modifications of the 
theory which will become clear as we describe our proce¬ 
dure. 


Let us begin by noting that we are not faced with a 
usual inverse scattering problem here where scattering 
data are given from which the potential is to be recon¬ 
structed. Rather, we are dealing with a self-consistent 
Schrodinger problem. An obvious way to attack the prob¬ 
lem is to exploit the relation of the potential with scat¬ 
tering data to reformulate the self-consistency equation 
in terms of scattering data and then determine the latter 
from it. Upon expressing quantities of interest through 
the scattering data, one can bypass the determination of 
the self consistent potential. This is the strategy we will 
pursue, focusing on the semi-infinite case L = oo and 
considering both the disordered phase m > 0 as well as 
the ordered one m < 0. 


Scale invariance enables us to scale the temperature 
variable m to ±1. Introducing 


V±{z) = Voo{2;±l) 


(4.1) 
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and noting that 

n±(oo) = (5±i,i, (4.2) 

we see that we must study the Schrodinger problem 

[- 8 ^ + v± (z) - (5±i,i]^/>(z, k) = k) (4.3) 

with 

= E- S±ia (4.4) 

on the half-line [0, oo) for potentials that vary as 

v±{z) + 0{z) (4.5) 

and approach the limiting values <5±ip sufficiently fast as 
z —>■ 00 (exponentially or as a power, depending on the 
sign ±) [7^ . 

Since v(z]L,m) corresponds to the energy density 
whose leading thermal singularity in the near-boundary 
region originates from the contribution ^ m?, the terms 
linear and quadratic in m must both be regular in m. 
Therefore the coefficients and must be inde¬ 

pendent of the sign of m. According to Eqs. (3.29) and 
(3.33|)- (|3.35|), we have indeed 


^^-1 = T7a-i(3) = ±- 


and 


= I'D = 


02(3) _ ^ 
4 


(4.6) 


(4.7) 


The spectrum of the associated Sturm-Liouville oper¬ 
ator on [0, oo) should be continuous and equal to 
spec('H„^) = [(5 ±i7,oo). Following a s tand ard approach, 
we introduce two dost solutions of Eq. (4.3) satisfying the 
boundary conditions 


/(z,±k) ~ e 


±ikz 


(4.8) 


where E(k) is a complex-valued function, the dost func¬ 
tion. For real k it can be written as 


F(k) 


Qcr(k)-ir;(k) 


k > 0, 


e‘"(-k)+i’7(-k) k<0. 


(4.12) 


For later use let us also note that F(k) can written as a 
Wronskian of two functions /i and / 2 , 


W[fi{z)j 2 {z)] = /i(z) 


8 / 2 ( 2 ) 8 /i(z) 


8 z 


8 z 


One has 


F(k) = W[f(z,k),^iz,k)]. 


/ 2 (z). (4.13) 


(4.14) 


This known result can be verified in a straightforward 
fashion by substituting (4.11) for Lp(z, k) in the Wronksian 
and evaluating the latter for z —k 00 , making use of its 
independence of z. 

For the sake of notational conciseness, we have re¬ 
frained here from adding subscripts ± to the above quan¬ 
tities ip, F, ..., rj. But it should be remembered that they 
all differ for the cases ± of the self-consistent potentials 
v±. 

In the ordered phase m < 0, L = 00 , the divergence of 
the susceptibility on the coexistence curve implies that 
the Schrodinger Eq. (4.3) must have an F = 0 state (poi 2 ) 
that approaches a nonzero value at z = 00 . It satisfies 


‘Poi^) = v-(z) (po(z). 


(4.15) 


This is precisely the equation for the order-parameter 
profile (4>(y, z)}L=ao,m=-i in the presence of an infinites¬ 
imal magnetic field oriented along a fixed direction. The 
square of the spontaneous bulk magnetization is known 
to be given by |m|/47r [TU]. Hence, if the eigenfunction 
(^o(z) is normalized such that 


lim tpo(z) = 1, 


(4.16) 


then the spontaneous order-parameter profile at L = 00 
is given by 


and normalize the so-called “regular solution” of Eq. (4.3) 
such that 


(p(z,k) = yz[l-kO(z)]. 
The latter behaves asymptotically as 


ip(z,k) ~ 


A(k) 


sin[kz -I- 77(k)], A(k) = e' 




(4.9) 


(4.10) 


which defines the scattering amplitude A(k) and phase 
shift ri(k). The regular solution can be expressed in terms 
of the Jost solutions as 

P(2^, k) = ^ [F{-k) /(z, k) - F(k) /(z, -k)], k > 0, 




\m\ 


^<o = 


(4.17) 


The small-z behavior ipo(z —k 0) ~ const y/z of po is sim¬ 
ilar to Eq. (4.9). At this stage we do not yet know that 


the proportionality constant “const” is exactly 1 and that 
(po agrees with the k —k 0 limit of po( 2 , k), i.e.. 


lirn(p(z, k) = poiz), 

k-fO 


(4.18) 


so that Po indeed satisfies the boundary condition (4.9). 
However, our inverse scattering analysis in Sec. IV C| 


will confirm the validity of both statements, i.e., that 
Eqs. ( 4.16 ) and (4.18) hold and po fulfills the boundary 
condition (4.9). _ 


Owing to the limiting behavior (4.16), po is not square- 


(4.11) integrable. Such a zero-momentum state, which is finite 
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dX z = oo but does not decay sufficiently fast to be square- 
integrable is called “half-bound state” m- It occurs in 
the semi-infinite case we consider here because the lowest 
eigenvalue Ei{L,—l) of the strip tends exponentially to 
zero as its thickness L ^ oo: one has [S21 ESI 


lim L ^ InEi{L,m) = m. 

L—^oo 


(4.19) 


The half-bound state arises for m < 0 from the eigen¬ 
function m) in the limit L ^ oo. Since the latter 

is orthonormalized and fiVL is dimensionless, the nor¬ 
malization (4.16) of tpo implies that 


lim fi{z; L,m) Vl = :poi\m\z). 

L—¥oo 


(4.20) 


The presence of an order-parameter profile for m < 0 
and L = oo entails an additional contribution to the L = 
oo analog of the self-consistency Eq. (3.23). The latter 
reads 


AE 


' —OO 
1 


0±(z,z; E) + 


1 


2^/^_ 


with 




G±{ 2 ,z.;E) = g{z,z-,E-oo,±l). 


(4.21a) 


(4.21b) 


Equation (4.21a) is a well-suited starting point for ap¬ 
plying inverse scattering theory. To this end we pro¬ 
ceed as follows. Let v±{z) be the above specified self- 
consistent potential for m = ±1 and L = oo, g± the as¬ 
sociated Green’s fun ction (4.21bI, and = V’o(z;—1) 
the solution to Eq. (4.15) for m = —1, normalized ac¬ 
cording to Eq. ( |4.16 1 . Hence Eqs. ( |4.21| ) hold for these 
quantities g± and :po{z). We now consider variations 


(z) —>■ v± (z) -I- 5v[z). 


(4.22) 


It will be sufficient and convenient to restrict the varia¬ 
tions 5v to functions with the properties 


(i) 

Sv{z) = regular for z > 0, 

(ii) 

lim 6 v{z) = 0, 

z—>-oo 

(iii) 

lim 5v{z) = Svq, 

z—>-0 

(iv) 

f dzSv{z) = 0. 

Jo 


(4.23) 


Equation (4.21a) simply corresponds to the continuum 


limits of the self-consistency conditions (2.32) and (2.53). 


Since the latter two are nothing else but stationarity con¬ 
ditions such as dfiq{t,'V)/dVj = 0, Eq. (4.21a) may be 
read as the stationarity condition 5f[v\/dv{z) = 0 for a 
free-energy functional f[v], which is the continuum ana¬ 
log of fN=ao{t,'V). Hence the linear variation 


6 f[v,Sv] = 


A A f 


(4.24) 


of this functional must vanish when evaluated at the self- 
consistent potential u*(z), where here and below an as¬ 
terisk is used to mark potentials u* and Green’s functions 
t/* t hat so lve the self-consistency equations. To express 
Eq. (4.24) in terms of the Green’s function, we must sim¬ 
ply multiply Eq. (4.21a I by 5v{z) and integrate z from 0 


to 00 . We then proceed by expressing the result in terms 
of scattering data. 

To this end it is convenient to introduce the energy- 
dependent scattering phase by 

VE = E g i6±i,i,oo), (4.25) 

for the cases m = ±1. This is because both the Green’s 
function and 6 v{z) can be eliminated in favor of ije and its 
variation Sr]E induced by Sv(z) by exploiting the relations 


dr]E 

dE 


dz 


Im0(z,z;E -I- iO) 


1 


2Ve_ 


(4.26) 


and 


POO 

SrjE = / dzdw(z)Im^±,,(z,z;£’-|-iO). (4.27) 

Jo 


The first relation, Eq. (4.26), may be viewed as a spe¬ 


cial case of the Lifshitz-Krein trace formula [ZilZT] 

Tr [h(H) - = -/ dEr]E h'{E) (4.28) 

for the trace of the difference due to a change of Hamil- 
tonians 77*-°^ —>■ H. In the original Lifshitz-Krein for¬ 


mula (4.28) h{E) stands for an arbitrary function h{E). 
Remembering that —Im tr[E-|-iO —77]~^/7r gives the den¬ 


sity of states of 77 one sees that Eq. (4.26) follows formally 
from Eq. ( |4.28| ) if we choose h{E) = s(e — E'). In Ap- 
pendix[B]we give a direct proof of Eq. (4.261. The second 
relation, Eq. (4.27), is a direct consequence of Eq. (4.26) 


and also derived in this appendix. 

Use of the above relations will enable us to derive from 


Eq. (4.24) an equation involving rjE, SrjE, and Svq. The 


first two quantities can be expressed in a straightforward 
manner in terms of the phase shift r](k) and its variation 
6 r]{k). Relating Svq to scattering data is more subtle but 
easily achieved with the aid of the following corollary of 
a trace formula proved in H [78] . 

Corollary IV.1 Letv{z) andv(z) be two continuous po¬ 
tentials on the half-line (0, oo) with the following proper¬ 
ties: 

(i) They vanish sufficiently fast for z ^ oo so that the 
integral 

POO 

/ dz|u(z)|<oo (4.29) 

J Zo 

and its analog for v(z) exist for all Zq > 0. 
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Cl 


1 ^ Rci 


FIG. 2. Integration paths in Eq. (4.321. 


(ii) They have the same singular behavior atz — 0 spec¬ 
ified in Eq. (4.5 ), with identical coefficients V-i and 


V-i though possibly different limiting values vq and 
Vo of their regular parts. 

(Hi) The Schrodinger equation TLytfifz) = etpfz) subject 
to the boundary condition ^/>(z —)■ 0) = 0 (y/z) and 
its analog with v ^ v have no bound-state solutions. 

Then the following relation holds between the difference 
of the latter coefficients and cr(k), the logarithm of the 
scattering amplitude introduced in Eq. (4.10); 


We set w = t; + and a = a -\- 5a va Eq. ( |4.30[ ) and 
linearize in 5a. Taking into account that a{k) is an even 
function of k, we thus arrive at the relation 

4 r°° 

5vo = — / dk (5CT(k). (4-31) 


The implementation of the above strategy requires sep¬ 
arate considerations in the cases to = ±1. We begin with 
the TO = - 1-1 case. 


B. Scattering data for the half-space problem with 

m = 4-1 


We then multiply by 5v(z), integrate o ver z, use the prop¬ 
erties (iii) and (iv) specified in Eq. (4.23), change the 
order of integrations in the double integral over z and 
E, and make use of Eq. (4.27) to perform the integra¬ 


tion over z. Making the change of integration variable 
E = k^ -I-1, we then rewrite the resulting integral dk^ 
as an integral dk, express tje and 5riE in terms of the 
phase shifts 77 (k) and 5r]{k), and take into account that 
the latter two quantities are odd in k. One thus arrives 
at 

/ OO II -I 

— kln(k^-I-1) ( 57 y(k)-I-- (5?;o = 0. (4.34) 

27r 4 


For 5vo we can substitute Eq. lol] ). However, the 
variations 5r]{k) and ( 5g(k) are not independent. To see 
this, recall from Eq. (4.12) that cr(k) and —??(k) are the 


real and imaginary parts, respectively, of the logarithm 
of the dost function, InE(k). Since the Hamiltonian TLv.^. 
has no discrete spectrum, E(k) has no zeros in the upper 
half-plane Imk > 0. Therefore, InE'(k) and (51nE(k) 
are analytic there. Moreover, (51nF(k) can be shown to 
vanish faster than |k|“^ as k —>■ oo in the upper half-plane: 
we show in Appendix that [7^ 


4 

Vq — Vo = — dk\c 

g-2o-(k) _ g-2o'(k) 

(4.30) 


TT Jo 



see Eq. (C19 


5\nF{k) = 0(|k|-3); 

k—voo ^ ' 

lmk>0 


(4.35) 


via the Kramers-Kronig relations 


5r]{k) = HT[(5cr](k) =-V [ dk' (4.36a) 

^ J — OO ^ ^ 


and 


5a{k)=RT-^[5r]]{k) = --r [ dk', (4.36b) 

^ J — OO ^ ^ 

where V denotes the Cauchy principal value and HT[/] 
is the Hilbert transform of / (s ee, e. g., §4.2 of [5D]b 
Another consequence of Eq. (4.35) is the constraint 


dkk5'i]{k) = 0. 


(4.37) 


We choose the plus sign in Eq. (4.21a) and rewrite the 


integral as a contour integral along the path C 2 shown in 
Fig.m obtaining 


d(z) 

2 


dE 


dE 

27ri 


^+,*(z,z, A) -f 
In A ^+_*(z,z, A)-I- 


1 

1 


. (4.32) 


2 v^r^. 

The contour C 2 can be deformed into Ci, which gives 
/■°° dE 


27r 


InE 


Imf/-l-,*(z,z, E-l-iO) -I- 


1 




4 

(4.33) 


This follows from the fact that the integral on the 
left-hand side can be written as —Im dkki51nE(k), 

which vanishes because the integration contour along the 
real axis can be completed to a closed loop in the upper 
complex half-plane by a half-circle of infinite radius. 

We now insert Eq. (4.36a) into Eq. ( |4.31 ). Upon ex¬ 
ploiting the antisymmetry of 5r](k) together with the con¬ 
straint (4.37), one sees that (5?;o can be written as 




dkk^5f]{k.)V / dk' 


g- 2 (T(k') 

k'-k ' 


(4.38) 


The result can be inserted into Eq. (4.34), but we must 


take into account the constraint (4.37), which we do by 
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means of a Lagrange multiplier K/2'k. Equating the coef¬ 
ficient of Sr]{k) to zero then gives us an integral equation 
for cr+(k), namely 


ln(k2 -h 1) 

k 


+ -r 

TT 



g-2cr+(k') 

k' - k 


^=0, (4.39) 


where we have reintroduced the subscript -I- to remind 
us that fT+(k) refers to the m = - 1-1 case. 

We claim that the solution to this equation is 

cr+ (k) = — ^ In (k~^ arctan k), (4.40) 

which implies the remarkably simple result 

A+(k) = y/k/ arctan k (4-41) 

for the scattering amplitude. 

To show this, let us dehne the functions 


The integral V ■ ■ ■ dk' equals the real part of 
• ■ -dk'. We now transform to the variable B = 

J—oo+iO ^ 

arctan kb Then the latter integral transforms into an in¬ 
tegral in the complex /3-plane along a path inhnitesimally 
above the real axis from — 7 r/ 2 -|-i 0 to — 7 r/ 2 -|-i 0 . This path 
can be deformed such that it becomes the sum of two inte¬ 
grals parallel to the imaginary axis, one from (— 7 r/ 2 -|-i 0 ) 
to (— 7 r /2 -I- ioo) and another one from (— 7 r /2 -|- iO) to 
(— 7 r/ 2 -|-ioo), and a vanishing contribution at infinity. It 
follows that 


, 1 1 '°° cothu du 

2 Jq k 2 -I- coth^ u + 7 r 2/4 ’ 


(4.48) 


where u = Im/3. 

Below we shall need the limiting behavior of ri'{k) for 
k —>■ oo. In order to determine this, we split off the 
analytically computable term 


1 e-2CT+(k') 

5^ ^ j_^ dk' ^^ 73 ^ for Im k ^ 0. (4.42) 

The functions < 7 > (k) and (k) are analytical for Im k > 0 
and Im k < 0, respectively, and vanish as k —>■ oo in the 
respective upper or lower complex half-plane. For real k 
one has 

5>(k)-5<(k)=e-2-+W ( 443 ) 


I r 


1 /u 


dw 


2 ln(k 7 r/ 2 ) 


k? + l/vB V? + TT^/d k^TT^ — 4 

_ (4.49) 

from the integral on the right-hand side of Eq. (4.48|). In 


the remaining integral, one can expand the integrand in 
powers of k~^. One thus obtains 


77 +(k^ 00 ) = 


In 


kTT 


J+ 


+ ^+0(k-"lnk) (4.50) 


and 


2 r°° e-2<T+(k') 

2%>(k)+ff<(k)] = -7^y ^dk'^^ 4 -^. (4. 


44) 


Using the latter equa tion to eliminate the principal value 
integral in Eq. (4.39), we are led to 


In(I-ik) A ln(l- 1 - 

-4-! + 2.s>(k) + ^ = —4- 


ik) 


- 2 ig<(k). 

(4.45) 

The left-hand and right-hand sides of this equation are 
analytical in the upper and lower complex half-plane, re¬ 
spectively, and vanish there as k — 00 . Hence both sides 
dehne a bounded entire function which must be identical 
to zero, just as A = 0. Hence 


5§(k) = - 


ln(I =F ik) 
2 ik ’ 


(4.46) 


from which the claimed solution ( 4.40[ ) follows at once 
via Eq. (O^ . As a consistency ch eck, one can insert 
Eq. (4.43) for (T+(k), set A = 0 in Eq. (4.39), and compute 


the principal value integral to conhrm that this equation 
is fulhlled. 


Since the Kramers-Kronig relations (4.36) carry over 


to ? 7 (|_(k) and a'j^{k), we can comput e the derivative of 
the phase shift r]+{k) from the result (4.40). One hnds 


y(k) = 


dk' 


1 

k'^ 


(k'- 


I 

I) arctan k' 

(4.47) 


with 


J+ 


f°° coth(u) — 1/u 

Jo ^ u'^+TT^/d 


(4.51) 


The integral can be shown to have the series repre¬ 
sentation m 


j+ = 


8 


00 


E 


ln( 2 k) 
4k2 - 1 ■ 


(4.52) 


Its value 


J+ = 0.474572605201526198964... (4.53) 


can be determined in a straightforward fashion by nu¬ 
merical computation of either the series or the integral 
using Mathematica [82] . 

In order to determine ry_|_(k), we integrate its derivative 
given in Eq. (4.48), using r]+{0) = 0. Upon changing the 
order of integrations and performing the integration over 
k, we are led to 


7?+(k) = / du 
Jo 


2 arctan(k tanh u) 
4u^ -I- TT^ 


(4.54) 


The remaining integral can be evaluated numerically. 
The result is plotted below in Fig. together with the 
corresponding phase shift 77 -(k) for m = — 1 , whose cal¬ 
culation we present in the next subsection. 
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From the results for cr+(k) and ry+(k), one can also 
derive a helpful representation of the logarithm of the 
Jost function for Im k > 0, namely 


lnF+(k) = r 

Jo 


dk' 


1 


arctanh—. 

W 


' (k'2 + 1) arctan k' 

(4.55) 

We have benefited from it in our analysis of the asymp¬ 
totic X —>■ -boo behavior of the free-energy scaling 
function Y(x), which is introduced and discussed in 
Sec. lyTml 

To obtain this representation, we computed cr(|_ (k) from 
Eq. ( |4.40 ), used the result for r]'^{k) given in Eq. ( |4.47 ), 
inserted these expressions into dlnF_|.(k)/dk = — 

if?+(k), and combined both contributions into an integral 
of the form [...]/(k' — k — iO). If one then chooses 

instead of the wavevector K = k-|-i0 the purely imaginary 
one K = ig, (lnF+(A:))'|k=i 5 becomes real-valued an d can 
be integrated to obtain lnE+(K = ig). Equation (4.55) 
finally follows by analytic continuation of the result to 
complex values of K. 


C. Scattering data for the half-space problem at 

m = —1 


In the case of m = —1, we start again from Eq. (4.21a|. 


Now the spectrum of the Hamiltonian becomes gapless, 
spec('H„_) = [0, oo), and the contribution oc ^Pq{z) must 
be taken into account. The presence of the half-bound 
state ifo{z) implies that the Jost function F{k) now has 
a zero at the origin, F{0) = 0. To proceed, we rewrite 
the integral in Eq. ( 4.21a[ ) as lime_>.o+ dz, transform 
it into 


-hm r^lnF 

e^0+ TT 


Img_(z, z;E + iO) -b 


and take the limit e —>■ 0-b. We then multiply by Sv(z) 
and integrate over z. We thus arrive at 


) ^oo 1 IT' 

dzSv(z) / -—InE 

If) ZTT 


Im t7_(z, z;E + iO) -b 


1 


2 Ve 


1 1 
= 2 Jo (4.56) 

The function Sv{z) is required to have the properties 
(i)-(iv) known from the m = -bl case and specified in 
Eq. (4.23). However, we also require that Sv does not 


shift the eigenenergy Ei = 0 associated with the half¬ 
bound state <b’o('2) to linear order in 6 v. This means that 
the expectation value ((po|<^^^|<Po) must vanish, i.e., the 
fifth property Sv{z) must satisfy is 


(v) 


pOC 

/ dz5v{z)Lp1{z) = Q. (4-57) 

Jo 


Thus the integral on the right-hand side and the contri¬ 
bution from 1/{2'/E) on the left-hand side of Eq. (4.56) 


both vanish. Setting E = k"^ and using Eqs. (4.27) and 
(4.31) along with the Kramers-Kronig relation (4.36a), 


we obtain 


f —OO 
1 


^ Sf]-{k) k In k^ = ^ (5uo 
ZTT 4 


dk' k'2 e-2<T-(k')p 


(4.58) 


Just as in the m = -bl case, we could replace the fac¬ 
tor k'^ by a factor k^ in the principal value integral. 
How ever, the subsequent reasoning that either side of 
Eq. ( |4.45[ ) equals an entire function is not applicable here 
because the integrand involves the function In k^ rather 
than ln(k^ -b 1), which is not analytic at k = 0. 


We assert that Eq. (4.58) has solutions of the simple 
form 


,-2rT_(k) _ 


A_ 


TT 


(4.59) 


where A_ G K is a parameter that needs to be deter¬ 
mined. In order to show this, we substitute this equation 
into Eq. (4.58|) and replace dk' by lim;f_>oo dk'. 


Since (5?7_(k) = 0(k“^) as k —)• oo by Eq. ( |4.35[ ), its 
Hilbert transform should exist and behave such that the 
remaining k' integral is convergent. Our rewriting of the 
k'-integral as a principal value was necessary to enable 
us to interchange the order of the two integrations. This 
leads us to 


-5vo 


/ °° dk 

Sr]-{k) V / dk 
-oo 7^ J-K 


, A_ -b 7r|k'|/2 


k- k' 


The second integral gives 

, A_ -b7r|k'|/2 


(4.60) 


fK 

V / dk' 
J-K 


k- k' 


o A 1 ^ k^ 

= 2A_ artanh-— -b — In - 


K 2 K'^-k^ 

= —Trkln AT-bIn k^-b 0(Ar“^). (4.61) 


Owing to the constraint (4.37), the kin AT term does not 
contribute to dk(5?7_(k).... Hence the limit AT —k oo 
of Eq. ( 4.60| ) indeed yields the left-hand side of Eq. (4.581. 
It remains to determine the constant A_. To do this. 


we compare the solution (4.59) with the asymptotic form 


e—«^=^--^u_-,TO(l/k^). (4.62) 

for large k derived in Appendix Insertion of the value 
of vZi given in Eq. (|4.6[) then yields 


A_ = I. 


(4.63) 
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Hence our results for cr_(k) become 


,-2a_(k) _ j_ 


k 2 2 |k| 


‘'-W = ^" l + t|k|/2 

1 7 rsgn(k) 


= lnH_(k), 


(4.64) 


To compute r]_{k) we first determine r]'_{k) using the 
Kramers-Kronig relation (4.36a). The Hilbert transform 
of the contribution 1 /k to cr!_(k), considered as a distri¬ 
bution, gives HT[l/k](k) = 7 r 5 (k). For k 7 ^ 0, it has no 
support and can be dropped. The Hilbert transform of 
the remaining term, cr'_(k) — 1 /k, can be computed in a 
straightforward fashion. One gets 


J7-(k) = 


ln(7r^k2/4) 
4 — TT^k^ 


(4.65) 


This can be integrated in a straightforward fashion. 
To fix the integration constant, we use the fact that 
a variant of Levinson’s theorem exists for the Sturm- 
Liouville problem we are concerned with which predicts 
that 77-(0±) = ± 7 r /2 in the presence of the half-bound 
state (see Eq. (4.43) of [75]). It follows that for 

real k 


I TT 1 

T]_{k) = sgn(k)<^ 2^^ 
-In 


Li2(-7r|k|/2)-Li2(7r|k|/2) 
2-7r|l 


In- 


2-k7r| 


(4.66) 


where Li 2 (k) is the dilogarithm (polylogarithm Lis(k) for 
s = 2 ) [SMS]. 



FIG. 3. The phase shifts r]-{k) for m = — 1 (red) and r]+{k) 
for m = 1 (blue). The dashed horizontal line shows the 
momentum-independent phase shift rio{k) — 7r/4 for m — 0. 


E. The half-bound state 


Using the above results, we can now also verify that the 
half-bound state (po agrees with the k —>■ 0 - 1 - limit of the 
regular solutio n (p(z , k), as stated in Eq. (4.18). Using the 
analog of Eq. (4.67) for the regular solution with m < 0 


along with our results for the scattering data given in 
Eq. (4.64), we arrive at 


e^r-lk) 

7 ’(z, k) ~ —— sin[kz -h 77 . (k)] 

z—>-oo K 

= (1 -I- k7r/2)“^/^ sin[kz -I- 77_(k)]. (4.69) 

Since 77 _ (k) = 7 r /2 -|- 0(k In k), the sine function becomes 
cos[kz-f O(klnk)] and approaches 1 in the limit k —>■ 0-I-. 
It follows that limk-ko <p(z, k) takes the value 1 at z = 00 , 
satisfies the normalization condition (4.16) of (po, and 
hence must be identical to it. 


D. The phase shifts 77±(k) for m = ±1 and the 
scattering data for 7n = 0 


Fi gure[^ shows our results for r/± (k) given in Eqs. (4.54) 
and ( 4.66| ). To plot 77 +(k), we evaluated the integral on 
the right-hand side of Eq. (4.54) by numerical means us¬ 
ing Mathematica [Hlj. We included the /c-independent 
result rjo(k) = 7rj4 for the critical case to = 0. This fol¬ 
lows from the fact that the exactly known eigenfunctions 
(p{z,k]L = 00 , TO = 0) of the critical potential (3.13) for 
d = 3 |37| behave as 


(p{z,k\oo,Q) = \/zJo(kz) ~ - 1 /— sin(fcz-I-7r/4). 

z-koo V 7rfc 

(4.67) 

Comparing with Eq. (4.10), we see that the associated 
ttering an 

Ho(fc) = e 


scattering amplitude and phase shift are given by 
_ p<To(fc) _ , f^\ 


77 o(k) =-sgn(k). (4.68) 
TT 4 


F. Distinct behaviors of the Jost functions Ut(k) in 
the complex plane 

When we determined above the scattering data for the 
cases TO > 0 and to < 0, we saw that the qualitatively 
distinct spectra of the respective Hamiltonians en¬ 
tail qualitatively different behaviors of the Jost functions 
F±(k). While F'_|.(k = k/\m\) has no zeros in the upper 
half-plane Imfc > 0, F-{k/\m\) has a zero at the origin, 
reflecting the presence of the half-bound state discussed 
in the previous subsection. The difference of the spectra 
for TO > 0 and to < 0 manifests itself also in the dis¬ 
tinct analytical properties of the functions F±{k/\m\). 
The gap of s^ec{%y^) implies that F^(k/m) is analyti¬ 
cal in the complex A:-plane except for a branch cut from 
—Ito to —ioo along the negative Im/c-axis (see Fig. |^). 
The distance to of the end point of this branch cut from 
the origin corresponds to the inverse of the bulk corre¬ 
lation length 1 /to. When to < 0, the branch cut runs 
from fc = 0 to —ioo along the negative Im k axis. From 
Eqs. (4.64) and (4.66) one can see that the Jost function 
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F-{k/\'m\) is analytical on the first (physical) Riemann 
sheet but has successively poles and zeros at ±2|m|/7r 
on higher Riemann sheets, depending on the number 
of 27r turns (see Fig &)• The distance of the pole at 


mulant by 

G^'^\y,z,z';L,m) = lim {4>{y,z) ■ (f){0,z'))/n 

n—>-oo 

= ($(y,z)$(0,z')), (5.1) 


Im k 


(a) 

F+{k/m) 


. —im 



FIG. 4. Branch cuts of the lost functions F±{k/\m\) for m > 
0 (a) and m < 0 (b). The red crosses in (b) indicate the 
poles and zeros that occur successively on the higher Riemann 
sheets. 


and 

= G^‘^Hy,z,z';L,m) - ($(y,z))($(0,/)), (5.2) 

where cf){y,z) and $(y, z) are the continuum analogs of 
(pj and <i>j, respectively, and we have taken into ac¬ 
count the translation invariance parallel to the bound¬ 
ary planes. Since the 0{n) symmetry cannot be sponta¬ 
neously broken at temperatures T > 0 when L < oo, the 
two-point function approaches zero as the separation of 
the two points tends to infinity; i.e., 

lim z, z';L < oo,m) = 0. (5.3) 

|3/|->.oo 

On the other hand, if we first take the limit L —>■ oo 
and subsequently the limit \y\ —>■ oo, then the correla¬ 
tion function factorizes into a product of nonvanishing 
expectation values if m < 0, 

lim G(^)(y,z,z';oo,TO) = ($(0,z))($(0,z')), (5.4) 

|y|->.oo 

so that Gcum(2/, z, z'; oo, m) vanishes in this limit. 

Let us indicate Fourier transforms with respect to y 
by a tilde, writing 

/(y) = /^/»e‘^'^. (5.5) 

Dimensional analysis implies the scaling behavior 


Refc > 0 from the origin sets the scale for the algebraic 
decay of the bulk correlation function and corresponds to 
the Josephson coherence length oc l/|m|, as our results 
for the two-point functions to be presented in the next 
section explicitly show. 


V. RESULTS FOR TWO-POINT CORRELATION 
FUNCTIONS FROM SCATTERING DATA 

A. Relation of two-point functions to scattering 
data 

We will now use the scattering data determined in the 
previous section to calculate the scaling function of the 
L = oo two-point function for m ^ 0. To this end we 
must express the scaling function in terms of scattering 
data. It will be helpful to briefly return to the case of 
finite thickness L. 

We define the two-point correlation function and cu- 


&^um(.PiZ,z';oo,m) 


= |^G(2)^(p/|m|, |m|z, |m|z';oo,±l) (5.6) 


for m ^ 0. 


The analog of this equation for 
and m < 0, given in Eq. (4.17), asserts 


that the half-bound state (po describes the dependence of 
this quantity on |to|z. 

Before we address this issue, let us first consider the 
L ^ oo limit in the simpler m > 0 case. To this end we 
start from the analog of the spectral decomposition (3.10) 
for G^^^ (p; z, z'; L, m), expressing it in terms of eigenfunc¬ 
tions that are normalized such that 


lim z (pi/(z; L, m) = 1, 

z—^ 0 + 

lim (L-z)“^/^VJi.(z;L,m) = (-1)'' 

z^L — 


(5.7) 

(5.8) 


and hence related to the orthonormalized ones, f^, via 

(pv{z\L,m) 


f^iz;L,m) = 


yJxuiL.m) 


(5.9) 
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with 


K^{L,m) = = / [(pj,{z;L,m)]‘^ dz. (5.10) 

Jo 

This gives 


G^^\p,z,z';L,m) = ^ 


ip^{z]L,m) (pu{z'\L,m) 


^ >c^{L,m) + E„{L,m)]' 

(5-11) 

Using Eq. (5.101, we can now determine the limit 
L —>■ oo of x^/L. We can replace (p,^{z; L,m) by its 

Jo dz... for 


asymptotic form (4.10) in the integral L 


all layers z that are sufficiently far from the boundary 
planes. Near the boundary planes, this approximation 
is not justified. However, the integral of the difference 
between the asymptotic form and ‘py is restricted to two 
boundary layers of thickness ~ l/|w|. Therefore, the er¬ 
ror this approximation produces for XyjL is of order 1 /L 
and we have 

lim L-^Xy{L, m) = (5.12) 

L —^oo 2il\, 

Since the spectrum is bounded away from zero when 
m > 0, i.e., Ey{L,m) > m for all v, we can use Eq. (5.12) 


along with L ^ 

considerations to conclude that the L = oo limit of 
Eq. ( |5.11[ ) is given by 

z'] oo, TO > 0) 

1 dk 2 k^e“^‘^+(''^ ip{mz, k; oo, 1) (p{mz', k; oo, 1) 
m Jo TT p^/to^ -|- 1 + k^ ’ 

(5.13) 


L—¥00 


> TT ^ Jo dk and dimensional 


where k = k/m. 

When TO < 0, we must be more careful. Since Ei 
0 as L 


00 according to Eq. (4.19), the contribution 


from the lowest-energy mode i/ = 1 requires separate 
considerations, even though the mode sum 
remaining modes can be treated as before. 

Using Eq. (4.20), we see that the contribution of the 
1 / = 1 mode to {y, z, z'; L, to) for large L can be writ¬ 
ten as 


1 - i /" dP'y (po{\m\z) (fo{\m\z') 

L-s-oo L J (27r)2 + Ei{L,m) 


^ Po(|to|z) po(|to|z') 

L — yoo QiTT L 

= ^</ 5 o(|w|z) <po{\m\z'). 


l\VEi{L,m) 


(5.14) 


Here Eq. ( |4.19 1 for Ei and the asymptotic behavior 
Ko{y —?► 0) ~ — 7 e — ln(y/2) of the Bessel function Kq 
were used, where is the Euler-Macheroni constant. 
The derivation nicely illustrates the emergence of the 
half-bound state ipo as L ^ oo and of the square of the 
spontaneous magnetization in the (noncommuting) limits 


lim|j^l_,,oo limi_>oo of The result confirms Eq. (4.17) 
for the spontaneous order-parameter profile. 

Combining it with the contribution from the mode sum 
yields 


&^Hp,z,z'-,oo,m) - &J 2 ^{p,z,z';oo,m) 
= ‘Po{\m\z)ifo{\m\z') 6 {p) 


(5.15) 


with 


Gi^v!miP,Z,z'-,00,771 < 0) 


1 

ItoI 


dk 




It 


k 2 


X (p(|to|z, k; oo, —1) ip{\m\z', k; oo, —1) 


B. Critical two-point cumulant 


(5.16) 


As a useful first application of Eq. (5.13), let us verify 
that Bray and Moore’s exact result |36ll37| for the critical 
two-point cumulant follows from it. Upon insertion of 
the scattering amplitude e‘^° given in Eq. (4.68) together 


with the critical eigenfunctions (4.67), we arrive at 

k Jo{kz) Jo{kz') 


{p, z, z'; oo, 0) = '/ziz' f 

Jo 


dk- 


Q p'^ + kP- 

= v^/o(pz<) A:o(pz>), (5.17) 


in precise agreement with [i 
tion space, the result 

G(^)(y,z,z';oo,0) 


lEZl- Transformed to posi- 


' zz' 


27r 1^2/2 + (z _ z'Y + ( 2 ; + z'Y 




,gi2) 


zz' 


y^ + {z- z'f 
4zz' 


(5.18a) 


takes the form dictated by conformal invariance under 
conformal transformations that leave the boundary plane 
z = 0 invariant [laiHZ], where the scaling function 
is given by 


G^^Hp) = 


1 


dTTy/p{l + p) 


(5.18b) 


These results mean that G^^^ complies with the famil¬ 
iar asymptotic large-distance behaviors of the form 


G(^)(2 ,,z,z0 ^ + 

r..2 I _ ^/^2l 




d-2 + rju 


d-2 + pj 


(5.19) 


as |£c| —^ 00 along a direction parallel to the surface plane 
or any nonparallel direction, respectively, where d = 3, 
p± = 1/2, and ryu = 1 in our case |55] . 
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The above findings suggest the introduction of the sur¬ 
face operator 


$s(2/) = lim z z) 

z^ 0 + 

and the associated correlation function 

G(l.l)(y _ = ($(y,z) $s(y'))- 


Equation (5.18) then yields 

G(i’i)(y,z) = 




27r(j/2 z^) 


(5.20) 


(5.21) 


(5.22) 


To define a finite two-point surface correlation function, 
we must first subtract the logarithmic singularity ~ In z 
before taking the limit z —?■ 0. We choose this subtraction 
as 


S(fiz) = -jE - ln(/rz/ 2 ) 


and define 


G(°’2)(p) = lim z-i/^G(1’1)(p,z) - S(fxz) 
2—>-0 L 

= -ln(p/^). 


(5.29) 


(5.30) 


for this function and the related surface two-point func¬ 
tion 


G(°’^)(y)^limz-i/2G(i’i)(y,z) = ^. 

y^o ^ 


(5.23) 


Here, the condition y ^ 0 in the definition (5.23) is 


necessary. To see this, note that the Fourier transform 
of Eq. ( 5.21[ ) exists and agrees with the z' —>■ 0 limit of 
(z')-i/^G^(p, z, z') computed from Eq. (5.18), namely 


From Eqs. (5.26) and ( |5.30| ) we see that the Fourier p- 
transform of this function is 

+ (7E - ln2)5(y). (5.31) 


C. Two-point cumulant in the disordered phase 

We next consider the case of m > 0. To determine the 
scaling function on the right-hand side of Eq. (|5.6[ ), we 
substitute the result (4.59) for in Eq. (15.13^ and 


make a change of variable k — >■ k = km. This gives 

2 k arctan(fc/TO) 


~ 1 c\k 

&'^''^'>{p,z) = ^/zKo{pz). (5.24) G(^)(p, z,z';oo,m) = — / — 

m Jo tt 


k^ 


However, G^°’^^(y), when considered as an ordinary 
function rather than a generalized one, does not have a 
Fourier p-transform. This is because in the derivation 
of Eq. (5.23), we made use of the requirement y ^ 0 


X ip{mz,k/m]oo,l)ip{mz',k/m]oo,l) 


(5.32) 


and therefore ignored potential contributions oc S{y). To 
resolve the problem, we should regard the above Fourier 
p-transforms as those of generalized functions in y-space. 

Generalized functions of the form = {yf -|- y|)^/^ 
with negative integer exponents A are discussed in detail 
in inni ■ Let us denote the generalized function associated 
with as 'P 2 y' where the subscript 2 reminds 

us that y € and introduce an arbitrary momentum 
scale p to make the integration variables dimensionless. 
Then V 2 {yy)~^ is defined via its action on test functions 
v{py): 


Note, first, that our results for the critical case 
are easily recovered from the right-hand side by per¬ 
forming the limits limm_j.o arctan(fc/TO) = 7r/2 and 
limm-s .0 (^(toz, fc/m; oo, 1) = y/~zJo{kz). Second, 

we can divide by and take the limit z' —)■ 0 using the 
normalization (4.9) of the regular solution to obtain 


G^^’^^(p, z;oo,to) 

1 dk 2 karctaii(k/m) (p{mz,k/m; 00 , 1 ) 

fm Jn TT 


+ k"^ + 


(5.33) 


V 2 


,V = 


/|y|<i y 


2 „z ^’iy) 


/ ^ y -2 

l\y\>i y 


We next compute the m > 0 analog of defined 

as in Eq. (5.30). To perform the required z —>■ 0 limit. 


(5.25) 

where y = py. Its Fourier p-transform is well defined 
and given by 


—FT 

27r 


r 2 

Lp 2 y 2 ’^ 


1 p 

= -yE-ln-. 


To proceed, note that the limit z —>■ 0 of 
z-^/^G<^P^\p,z) = z-^G<^^\p,z,z) = Koipz) 
does not exist because 

Ko{pz) = - 7 e - ln(pz/2) -|- 0{p^z^). 


(5.26) 


(5.27) 


(5.28) 


we represent the subtraction term as 

S{y-z) = j -f 0[(pz)2ln(pz)], (5.34) 

substitute arctan(A;/TO) by 7 r /2 — arctan(m/fc), use 
lim^_>o( 2 m)“^/^ip(mz, fc/m; oo, 1 ) = 1 , and exchange the 
limit z —>■ 0 with the integration. This gives 


G*'°’^^(p; 00 , m > 0 ) = 


poo 

/ dk 

k 

k 

p"^ + k'^ + m? 

yj + k‘^ _ 

/o 


2 

TT 


,, , arctan(77i/fc) 

/ dkk — -—-- 

'q -I- -I- 


= In 


. -I- a/-I- p^ _ 


(5.35) 
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To determine the Fourier back transform, it is conve¬ 
nient to consider the difference 


For the corresponding Fourier back transforms one 
finds 


G^°'^\y;oo,m > 0 ) — G^°’^\y-oo,0) 


d^p 

( 2 ^ 

dp 

27r 




m 


2 m? + p 2 


pJoipy) In 


- -I- a/ m? -b p^ 


1 


2 TTy‘ 


(^Q-my _ 


(5.36) 


where the last line follows from the integral liste d as 6 .775 
in [^. Upon combining the result with Eq. (5.31), we 
obtain 


G(°’2)(y;oo,m > 0) = e -b (te - In 2) (5(y). 

(5.37) 

If one assumes y ^ 0, then (y; oo, to > 0) can be 

considered as an ordinary function and both the symbol 
V 2 as well as the contributions cx 5{y) be dropped. 


D. Two-point cumulant in the ordered phase 


(y; 00 , TO < 0) = (y; 00 , to < 0) -b 


V 2 


27ry2 2'Ky 47r 


m\ m , , , 

+ u-^ (7 e - ln2)(5(y). 


47r 

(5.42) 


In Table H] we compare our exact n = 00 results for the 
surface correlation function G^°’^^(y; 00 , to) with those 
known for G^‘^\y,z,z' ;oo,to) and the bulk correlation 
function G^\x-,m) = {^{x) $(0))b. 


m > 0 


G^'^\y,z,z') 


G<°'^)(y) G^^\x) 

1 Q~'^y 1 0 “^^ 

27ry^ 47TX 


m = 0 





1 


4Trx 


m < 0 


1 _|_ |m| I 

27ry^ 2'Ky 4 k 


_j, _I IH 

4kx 4k 


TABLE I. Comparison of exact two-point correlation func¬ 
tions G^^\y,z,z'-oo,m), G^^'^\y, oo,m), and bulk correla¬ 
tion function G^\x-, m), where z± = z ± z' and x — (y, z). 


U pon inserting the scattering amplitude from 
Eq. ( |4.64 1 into Eq. (5.16), we arrive at 

G^^umiP^z,z']O0,m < 0 ) 


TT 


2 k /|™| TT 
+ k'^ \ fc 2 


X (p{mz,k/m; 00 , —1) (p{mz\k/m]oo, —1) 

(5.38) 

Just as in the to > 0 case, we can multiply by {z')~^G 
and perform the limit z' —>■ 0 to obtain 


GcumHP:^;00,TO < 0) 


1 


dk 


2 k ip{mz,k/m; 00 ,- 1 ) f\m\ tt 


2 , 
(5.39) 


to| Jo ^ I 

The calculation of the two-point surface cumulant, de¬ 
fined by analogy with Eq. (5.30), is similarly straightfor¬ 
ward. It yields 


Gcum^P; oo,TO<0) = — - In- (5.40) 

p p 

and hence 


VI. RESULTS FOR SOME 
POTENTIAL-RELATED QUANTITIES 

A. The potential coefficient vg 

In order to compute the potential coefficient we 
use the following variant of a trace formula derived in 
the accompanying paper II; 

Theorem VI.1 (Trace formula) Let v(z) and v{z) be 
two potentials on the half-line ( 0 ,oo) that vanish faster 
than z~^ as z —>■ oo and behave as 

v{z -)■ 0 ) = ?;^®®^(z) -b ro -b o(l), 

{)(z-)► 0 ) = ?;^®®^(z)-b fio + o(l)) ( 6 . 1 ) 

with the same singular part 

= + ( 6 . 2 ) 

Let Lpu{z), V = 1,... ,nb, be the regular solutions to the 
Schrodinger equation T-LvPvii^) = SvPvY subject to the 
boundary conditions 

= Vz[l + 0(z)] (6.3) 

z—)-0 

that correspond to bound states (where may be zero). 
Denote by 


oo, TO < 0) = — In — -b 6 {p) —. (5.41) 

p y 47r 





dz 


(6.4) 
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the squares of the L^([0, oo)) norms of these (real-valued) 
functions and by A(k) = the scattering amplitude 

for k > 0. Furthermore, let (fp, ht, kp, and A(k) = 
be the analogous quantities pertaining to the potential 
v(z). Then the following relation holds 


Vo-Vo = - 
TT 


dz k^ 
2 


bound ^ 
states u 


g-2CT(k) _ g-2CT(k) 


bound 
states y 


(6.5) 


We shall use this relation for the (m = ±1,L = oo) 
choices of self-consistent potentials u(z) = u^(z) — (5±ip 


with the singular parts 






„sg,± 


- 4z2 ± 


^ ( 6 - 6 ) 


and vfz) = u®®’^(z). The potentials u®®’^(z) weakly 


violate the required z —> oo behavior v(z) = 
o(l/z), where the “little-o symbol” / = o{g) means that 
lima._>_oo \ f{x)/g{x)\ = 0. If fulfilled, the condition 
would ensure (together with the presumed boundedness 
and smoothness of u for z > 0) that J^q |u(z)|dz < oo. 
The consequence of the slow 0(l/z) decay is known 
from the Coulomb potential: the phase shift rj(k) is not 
well-defined and gets replaced by z-dependent quantity 
ij(k,z) containing a contribution ~ In(zk), as we shall 
explicitly confirm below. Since the scattering theory for 
Coulomb potentials is mathematically well controlled, we 
may trust that Eq. (6.5) carries over to the border-line 
case when v(z) = 0(l/z), except for the mentioned re¬ 
placement of the phase shift by a logarithmically diver¬ 
gent 77(k,z). 


t. The case v = and v = 

For the choices v(z) = u+(z) — 1 and u®8’+(z) the 
Schrodinger equations have no bound states. The spectra 
are continuous with e = k^ € (0, oo). The corresponding 
regular solutions (p(z, k) are given by 

ei7r/4 2 

(p(z, k;K) = —^M_i„,o(-2ikz), k = —, (6.7) 

V2k T^^k 


where M_iK,o(z) is a Whittaker-M function 
asymptotic behaviors at small and large z are 


. Its 


(fi{z,k;±K) = Jz 

z-J-O 


1 ± -1! P-+ 0(.» 


TT^ 


and 


with 


A(k) 

ip{z,k;±K) ~ —^sin[kz-b? 7 (k,z)] 

z—>-oo K 


^(k) = - 

TT 


4 

l + exp(±-j 


( 6 . 8 ) 

(6.9) 

( 6 . 10 ) 


and 


z) = ^ + argT Q ± T ln(2kz). (6.11) 

Here the upper signs appl y to the case of the regular 
solutions specified in Eq. (6.7). The lower signs give 


the analogous results which hold when the sign of k is 
changed. They will be needed in our analysis of the case 
V = vZ, and v = presented below. Note that the 


results given in Eq. (6.11) exhibit the above-mentioned 


logarithmic z-dependence of 7)(k,z). 

We insert these results along with the scattering am- 
plitude_A+(k) given in Eq. dHT] ) into the trace for¬ 
mula (6.5). This gives 


- 1 = -b- 

TT 


dkk' 


7r/k 


arctan k 


1-I-exp[4/(7rk)] k 

( 6 . 12 ) 

Upon making a change of variable k —>■ g = 1/k, the inte¬ 
gral on the right-hand side can be decomposed into a sum 
of two finite integrals Ji and J 2 . The first can be rewrit¬ 
ten as an integral of a series, which can be integrated 
termwise and subsequently summed to obtain 




0 T 

OO OO 


=iy -*’1: 


1 -b 2 

(2j-l)-^ 


56C(3) 


-b 7r'^(2j — 1)2/16 tt^ 


(6.13) 


= -l. (6.14) 


Its contribution to UqMs compensated by the —1 on the 
left-hand side of Eq.(6.12). Hence we arrive at the value 
of given in Eq. (4.7). 


The second integral is 


J 2 = - / ^ 


r dg 

'TT 

/ i 

/o y 

— — q — arccot q 


2. The case v = v and v = u®®’ 

For the potential v{z) = vf^(z), no bound states exists. 
The spectrum is continuous with e = k^ > 0. However, 
in the case of the potential v(z) = vZ(^) infinitely many 
bound states exist. To see this, note that the r egular so¬ 
lutions satisfying the boundary condition (6.3) are again 


given by Eq. (6.7) except that k must be replaced by 
—K = —27r“^k~^. For k/i = k" > 0, the requirement that 
tf{z,\k"-,—K) decays exponentially as z —>■ 00 yields the 
discrete values and associated eigenenergies e,, = —(k/)^, 
namely 


kC = 


-16 


7r2(2zz- 1)’ 


= 


TT^ (2iz — 1)2 ’ 


V = 1,2 ,... , 00 . 

(6.15) 
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The associated eigenstates simplify to 
= (^(z,ik";-K) 


= \/ z exp 


4z 


7r2(2i/ — 1) 


Lv-i 


8z 


■k'^{2v — 1 ) 


(6.16) 


where the L^-i{z) are Laguerre polynomials. 

Using these results, one can compute the required 
norm parameters One obtains 


= - 1 )^, z/ = 1 , 2 ,..., oo. 

64 


(6.17) 


The regular solutions associated with the continuous 
part {e = k ^|0 < k < oo} of the spectrum are give n b y 
(^(z, k; —k) in terms of the functions defined in Eq. 
Insertion of the above results and Eq. (4. 


exp(—2(5-(k)) into the trace formula (6.51 yields 

Trk ^ Trk 


given uy 
Eq. fert . 
4.64^ for 


dk 


_1 + exp[—4/(7rk)] 

1 


- 1 - 


V 

TT^ {2v — 1)3 

56C(3) 112 C(3) 

A ' A 


56C(3) 


(6.18) 


Upon inserting the known values (4.6) of we can de¬ 
termine a± by matching the derivatives 7 ? l|-(k) o f these 
results with the large-k forms given in Eq. ( 4.50| ) or im¬ 
plied by Eq. (4.65). One thus obtains 


and 


a+ = J+ — a_ 


a_ = — [ 7 e - 1 -hln(16/7r)], 


( 6 . 20 ) 


( 6 . 21 ) 


where J+ was defined in Eq. (4.51). 


The universal amplitude difference Azlg follows upon 
insertion of these expresssions for a± into Eq. (3.52) and 
use of the numerical value given in Eq. (4.53). One 
finds 


(s) _ a+ + a- 


= 


167r 


IOtt 

= 0.00944132199. 


( 6 . 22 ) 

(6.23) 


(s') 

This exact value of AAg ^ derived here has recently been 
confirmed quite accurately in |33j by numerical solutions 
of the self-consistency equations for two microscopically 
different 0 {n) (jA models representing the same univer¬ 
sality class (called models A and B there). 


To compute the integral in the first line, we again trans¬ 
formed to the variable q = 1 /k, represented the resulting 
integrand as the series 

OO 

^ (-87r-V")[9" + (M^V4)"]”\ 

and integrated termwise to obtain the first term in the 
last line of Eq. (6.18). The series in the second line can be 


summed in a straightforward fashion. Hence we obtain 
the value of given in Eq. (4.7). 


B. The integrals a± and the universal amplitude 


difference AA, 


(b) 

0 


We next turn to t he calculation of the integrals a± 
defined in Eq. (3.53) and the related universal ampli¬ 
tude difference AAg . To this end we use the following 
strategy. We determine the asymptotic large-k behavior 
of the dost function F(k) by means of the semiclassical 
expansion. Details of this analysis can be found in Ap¬ 
pendix The asymptotic form of this function for the 
half-line case of a potential with a singular part z;®® of 
the form specified in Eq. (6.2) involves u_i-dependent 


analogs of the integrals a±. The corresponding results 
for the phase shifts read 

?l±(k) = [u^i( 7 e -l-ln 8 k) -ka±]-kCl(k-3). (6.19) 


C. A surface excess quantity associated with the 
squared order parameter 

For local quantities of semi-infinite systems whose z- 
dependent thermal average approaches the bulk value 
Qh = Q(oo) sufficiently fast, one conventionally defines 
the associated surface excess as 1121031 


Qs= dz[Q{z) - Q{oo)]. 


(6.24) 


We have already seen in Sec. IIIC that this definition 
entails UV singularities in the continuum limit a —>■ 0 
in the case of the excess energy density because of the 
behavior ^ z~^ of Vao in the regime a z \m\. These 
singularities required the subtractions made to define the 
integrals a± in Eq. (3.53). 

In the case of the squared order-parameter density 
$( 2 )(z) = ($(y,z)) 2 , no problems arise from the behav¬ 
ior in the near-boundary regime z < \m\ when a = 0 . 
As we know from Eq. ( |4.I7[ ), the scaling function asso¬ 
ciated with the density 47 r$(^)(z) is given by the square 
of the semi-bound state (/?o(|m|z). This vanishes linearly 
in z as z —)■ 0 and hence is integrable at the lower limit. 
However, problems arise at the upper integration limit 
L = oo when m < 0 because the existence of Gold- 
stone modes in the ordered phase implies the asymptotic 
behavior ^ 0 -('=*- 2 ) ^t distances much larger than the 
Josephson coherence length |m|“^ [43]. At d = 3, this 
decay ~ z“^ is too slow to ensure convergence at the up¬ 
per integration limit cx). In order to define a finite surface 
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( 2 ') J- 1 

excess $s , the definition (6.24) must be appropriately 
modified. 

Let us begin by explicitly demonstrating the claimed 
slow decay, showing that the regular solution ^pq{z) func¬ 
tion behaves asymptotically as 


~ 1- —-fo(l/z). 

Z—>-oo 4Z 


(6.25) 


To determine this limiting behavior one can solve 
the m = — 1 variant of the Schrodinger Eq. (4.3) for 
z ^ 1. This requires knowledge of the asymptotic 
form of the potential v^(z). In the next section, 
we determine the behavior of the self-consistent po¬ 
tential v{z', L^m < 0) = L~^ v{z/L;l,x) i n the regime 
1 ^ z\m\ |a;|. The result, given in Eq. (7.17), means 
that 


^ ^ z—^oo 2z^ 


(6.26) 


Upon substituting this asymptotic form for the poten¬ 
tial, one can determine the two linearly independent so¬ 
lutions 


(/3o(z —>■ oo) ~ (2z)^/^ Jif •\/2/z] 
ipoi^ —>■ oo) ~ —7r(z/2)^/^ Yi ( \/2/z ] 


of Eq. (4.3), where Jy('u) and Y^{u) are Bessel functions 
of the first and second kind, respectively. We have nor¬ 
malized these solutions such that </ 5 o(oo) = 1 and the 
Wronskian satisfies IU[(^o(z), V'o)^)] = 1. The first so¬ 
lution has the asymptotic behavior ( |6.25 ); the second 
behaves as 


ipo{z -)>oo) = z-f^-kC'o-b o(l/z) 


(6.28) 


with Co = (1 — 27 e + ln 2 )/ 2 . 

The result (6.25) tells us that a finite universal sur¬ 
face excess quantity associated with 47 r $^^)(z)/|to| can 
be introduced as the integral 


J ^ — 


dz 


>pI{z) - 1 -f 


2 z 


(6.29) 


To compute it, we proceed as follows. We expand 
the regu lar solution i^{z,k) of the m = —1 Schrodinger 
Eq. (4.3) to O(k^), writing 


(^(z, k) = (^o(z) + k^ 1^2(2) + 0 {k'^). ( 6 . 30 ) 

The function (^ 2 ( 2 ) satisfies the differential equation 

'Hyip2iz) = [- 8 ^ + v{z)]ip2iz) = ipo{z). ( 6 . 31 ) 
A special solution is 

</52(z) = <^o(z) / dC(/5o(C)V'o(C)-'0o(z) / dC(/3o(C)- 
Jo Jo 

( 6 . 32 ) 


The general solution is the sum of this solution and the 
general solution c,^(/?o(z) + c. 0 ^/>o(z) of the homogeneous 
equation. However, since </J 2 (z) = ^dl(p{z,k — 0), it 
must vanish faster than (p{z, k) for z —>■ 0, namely be 
0(z^/^). Both solutions (po and V'o of the homogeneous 
solution violate this condition. Hence =_c^= 0 and 


the correct solution must be given by Eq. (6.32). 


Integr ating the a symp totic form of (po{z)ipQ{z) implied 
by Eqs. (6.25) and (6.28), one finds that the first integral 


on the right-hand side of Eq. (6.32) behaves as 


[ dC‘^o(C)V'o(C) = 

Jo ^ 


Inz 3 


0 {ln^ 


(6.33) 

To determine the asymptotic behavior of the second in¬ 
tegral, we split off J^p and then compute the integrals for 
large z using again the limiting form of (/jq- This gives 




-r 

dC 

<Po(C) - 1 

[ dz 

1 - 

^(C-I)l 

Jo 


2C J 


g(C-i) 

2C 


= z - - Inz 4--I-0(l/z). 

>-oo z 


(6.34) 


(6.27) Insertion of these results into Eq. (16.32^ then leads us to 


+ |lnz - ^)z-t-0(ln2z). (6.35) 

Since the integral appears in the asymptotic large- 
z form of the regular solution’s k^-term, it should also 
appear in the corresponding k^-term of the dost func¬ 
tion F{k). Hence we should be able to determine by 
computing this function and matching it with our exact 
scattering-data results of Sec. IV C given in Eqs. (4.64) 


and (4.65). To obtain F(k), we first calculate the dost so- 
lution /(z, k) for large z and then use it together with the 
above results fo r the regular solution to compute F{k) as 
the Wronskian (4.14) of these functions. 


The dost solution /(z, k) can be determined at large 
z by means of perturbation theory. We use the familiar 
equivalence of the differential Eq. (4.3) with the integral 
equation [55] 

/(z, k) = + dCu-(C)-— ~ /(C, k), (6.36) 

which we iterate to first order in replacing /((^, k) and 
in the integral on the right-hand side b y e*''^ and the 
asymptotic form of given in Eq. (6.26), respectively. 
We thus arrive at the expansion 

/(z,k) = (^1- ^e-'‘'^[7r-l-iEi(2ikz)] -I-..., 

(6.37) 

where Ei( 2 :) denotes a standard exponential-integral 
function. For the contribution of jth order in the 
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potential, we have the bound (cf. p. 9 of 


|/(^■)(z,A:)| <conste-^“(^"^ 


dC 


Ckoo(C)l 


j\ LJz ^ 1 + I^IC J 

(6.38) 

where “const” is an appropriate (positive) numerical con¬ 
stant. Using once more the asymptotic form of we 
can determine the large-z behavior of the integral. Since 



Ckoo(C)l ^ J_ 

1 -|- |fc|^ z— >-00 2z 


|k| arccoth(l -f 2 |k|z) 


1 

4|k|z2 


+ 0(z-3), 


(6.39) 


the contributions of second order in are smaller by a 
factor of z“^ than the O(v^) term. Note that the latter 
varies as 


approach in the limit x = mL —> —oo those of a one¬ 
dimensional free massless bosonic field theory on the strip 
(0, L) subject to Neumann boundary conditions. This 
was achieved by generalizing previous arguments for the 
semi-infinite case [94] to map the low-temperature be¬ 
havior of the model to a nonlinear sigma model (see Ap¬ 
pendix B of |33j ) and confirming the results via numerical 
solutions of the self-consistency equations. The aim of 
this section is to compute the leading corrections to the 
limiting Neumann eigenvalues and eigenvectors for large 
negative values of x. 

To this end it will be convenient to scale L to unity 
and use the scaling properties of and implied by 
dimensional arguments to express them as 

E^{L,m) = E^{l,mL) = E^{x), 

= ^uiz/E, l,mL) = f^{i,x), (7.1) 


/W(z,k>0) 


—e 


ikz 


8kz2 


+ o(z-^) 


(6.40) 


in terms of the dimensionless variables x and 3 = z/L. 
Their limiting large-cc behavior stated above is 


for large z, in accordance with the bound ( |6.38 ). 

We now wish to compute the Wronskian 
W[f{z,k),(p{z,k)], expand it to O(k^), and take 
the limit z —>■ 00 . To this end it is sufficient to retain 
in /(z, k) only the contributions to 0{v ^) explicitly 
given on the right-hand side of Eq. (6.37) because the 
foregoing considerations imply that the terms with 
j > 2 do not contribute to this limit. Proceeding in this 
manner, we find 


lim E^{x) = E^{-oo) = E^ , 

X—>• —00 

lim f,z( 3 ,a;) = f,,( 3 ,-oo) = f^( 3 ), (7.2) 

x—y — oo 


where 


E^ = kl = [K{iz-l)f, 1 /= 1,2, ...,oo, 
fi(3) = l, f,z>2 =’/2cos(fc,,3), (7.3) 


lim W[/(z, k),(^(z,fc)] 

z—>-oo 


= -ik-h k^ 


iTT ^ 7 e ln(2k) 

4 2 2 


o(k2). (6.41) 


On the other hand , t he ex act Jo st fu nction F{k > 0) 
given by Eqs. (4.12), (4.64), and (4.66) has the expan¬ 
sion 


E(k) = -ik-h k^ 


1 . TT ln(7rk/2) 

2 ^4 2 


o{k?). (6.42) 


Matching these results finally yields the analytically ex¬ 
act and numerical values 


Jv = -^[1 + 7E + ln(4/7r)] 

« -0.9093900700860.... 


(6.43) 


VII. ASYMPTOTIC BEHAVIOR IN THE 
LOW-TEMPERATURE SCALING LIMIT x ^ -00 


denote the familiar Neumann eigenvalues and eigenfunc¬ 
tions. 

Perturbation theory about the x = —00 limit is not 
standard because it involves several unusual, though in¬ 
teresting and challenging, features. Eirst, since Ei = 0 
and El vanishes ^ as x —)■ — 00 , we must deal 

with an asymptotic zero-energy mode. Second, pertur¬ 
bative corrections to the potential must be determined 
self-consistently along with those of the eigenvalues and 
eigenvectors. Third, the perturbation expansions of 3 - 
dependent quantities differ, depending on whether 3 be¬ 
longs to the boundary regions of thickness ^ near 

the two surface planes at 3 = 0 and 3 = 1 or else to the re¬ 
maining inner region. Thus the problem involves aspects 
of singular perturbation theory and bears similarities to 
boundary layer theory in fluid dynamics 

To set up pe rturba tion theory, we start from the scaled 
analog of Eq. (3.23) corresponding to the replacements 
z —)■ 3 = z/L, TO —)■ X = tL, and L —)■ 1, respectively, 
and choose 3 G ( 0 , 1 ) so that the ^-distributions do not 
contribute. This gives 


A. Asymptotic behavior of eigenvalues and 
eigenfunctions 


X = 



^(3)3;^) + 


1 

27 ^]’ 


3G(0,1). (7.4) 


In [32] and [33] we showed that the eigenfunctions 
f J, and eigenvalues E^ of the se lf-con si stent Schrbdinger 
problem defined by Eqs. (3.7), (3.11), (3.12), and (3.23) 


Since the left-hand side diverges ~ x as x —>■ — 00 , so 
must the right-hand side. The latter divergence arises 
due to the behavior InEi ^ — |x|. We therefore subtract 































27 


from both sides of the equation the term 

f?( 3 ,a:) lim \ [ dE -^^ 

%^0+U-oo [E-Ei{x) E\ 

= \nEi{x), 

obtaining 


In p 


X- flii,x) lnEi{x) = P(3,x) 


(7.5) 

(7.6) 


with 
P{hx) 
= lim 

p-i-0+ 

where 


dE 




2^/^ 


- fiii.x) Inp 


rf, j fi(3)fi(30 , fi'(3) f!^(3) 

G[l,l \E) = - ^ -+ 2^ 


E ^ E - El, 

i '—2 


(7.7) 


(7.8) 


As X —>■ —oo, the function fi( 3 ,x) —>■ 1 while Q ap¬ 
proaches the Green’s function ^nn for Neumann bound¬ 
ary conditions at 3 = 0 and 3 = 1. The latter function is 
given by 


1 / nn ( 3 , 3 '; E = -p^) = - 


cosh(p 3 <) cosh[p(l - 3 >)] 
p sinh p 


(7.9) 

It can be determined by means of standard image meth¬ 
ods, which yield it as a sum of contributions from in¬ 
finitely many image charges. Using this representation, 
one finds that the x —>■ —00 limit of the integrand on the 
right-hand side of Eq. (7.71 becomes 


0nn(3!3; e) 


1 




3-2p3 


- CXJ 

+ —J 2 [^ e-2p^' -F -F e-2pO+3) 

j=i ^ 


P=\/-E 

(7.10) 


with 


«'( 3 ,x) = 1 -h 


ip(^) + V’(l - 3 ) + 27e -I- 4 In 2 
4|x| 


-I- o(|x 


-1\ 


(7.13) 

The result tells us that for arbitrary 3 S (0,1), the 
eigenfunction fi becomes a constant in the limit x —>■ 
— 00 , so that fi( 3 , — 00 ) = 1 indeed, as stated in Eq. ( [73t . 

In order for this perturbation expansion to be valid, 
the 0 (I/|x|) must be small compared to the leading 0 ( 1 ) 
contribution. Expanding in 3 gives 




- -41n2 
3 


0 ( 3 ") 


+o(l/|x|). (7.14) 


Owing to the symmetry 3 —> 1 — 3 , a corresponding 
expansion holds about 3 = 1. Hence, for the perturbation 
expansion (7.13) to hold, 3 must not belong to boundary 
layers of thickness ~ l/|x| near either one of the two 
surface planes at 3 = 0 and 3 = 1. In fact, very close 


to the surfaces we know from Eq. (5.7) that fi( 3 ,x) oc 
(^i( 3 ,x) vary ^ and ~ respectively. 

The crossover between the behavior in the boundary 
layers and the remaining “inner region” is expected to 
occur at a distance &/|x| away from the surface planes, 
where b should be of order unity. In terms of the un¬ 
sealed variable 0 this means that for given m and L 
with 6 /|m| <?; T/ 2 , the crossover from the near-boundary 
behavior to the behavior in the inner region occurs at 
0 « 6 /|m| and I — z « 5/|m|. The length 5/|m| agrees 
up to a factor of order 1 with the Josephson coherence 
length [M] . 


The perturbative result for fi given by Eqs. (7.13) and 


(7.14) can be used to obtain the expansion of the po¬ 
tential in the inner region. To this end, we solve the 
differential equation for fi, obtaining 


x( 3 ,x) = 


^3fi(3,a: 

fi(3,a:) 


+ Ei(x), 


(7.15) 


and then substitute the above result for fi. This gives 

^ ^ V’"( 3 ) + V’"(l- 3 ) , n /I l^ 

x( 3 ,x) = - — --l-o(l/|x|) (7.16) 


Upon substitution of this expression into Eq. (7.7), the for 6 /|x| < 3 < 1 — b/\x\ (inner region) and 


integral and limit p —)■ 0-1- that P[i,x) involves can be 
computed in a straightforward fashion using Mathemat- 
ICA [52]. One obtains 

P( 3 ,-oo) = 7 e -f 21n2 -H ^'0(3) + (7.11) 

where ' 0 ( 3 ) = dlnr( 3 )/d 3 is the digamma function. 
Hence we have shown that InUi = —|x|[l -I- o(l)]. 


-’( 3 , 2 ;) = 


1 


+ C(3) + 0(3") 


+ o(l/|x|) (7.17) 


for 6 /|x| ^ 3 1/2. By contrast, in the near-boundary 

regime z = z|to| 1 , the asymptotic behavior of the 
potential w( 3 , x) = x(l — 3 , x) must be in conf ormi t y wit h 
that of the semi-infinite case, i.e., with Eqs. (4.5)-(4.7). 


ai—>• —00 


Se cond , we can solve Eq. (7.6) for fi( 3 ) and then use 
Eq. (7.11) to conclude that 


We next turn to the determination of the asymptotic 
behavior of the lowest eigenvalue Ei{x) for x —>■ — 00 . To 
this end we write the orthonormality relation for fi(z, x) 
as 


fi(3,a;) = 


lx- P( 3 ,x) 
In Ei{x) 


In El (x) 


4 /( 3 , x) (7.12) 


. 1/2 


1 = 2 


= I<{b,x) + I>(b,x) (7.18) 
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with 


and 


pb/lxl 

I<{b,x) = 2 d 3 f?( 3 ,a:) 

^0 


/*V2 

I>{b,x)=2 difli^,x), 
Jb/\x\ 


(7.19) 


(7.20) 


where b^ 1 but b/\x\ < 1/2. To comp ute J> we ca n sub - 
stitute our perturbation results, Eqs. (7.13) and (7.14), 
for fi. This yields 

. 1/2 


/>(6,x) = 


2x 




1 


ln£;i(ai) 

, - 1 /( 3 ) + 1/(1 - 3 ) + 27 e+ 41n2 


2 a; 


a(l/|ai|) 


-1 

In i?i(a;) 


46 


|x| + 7 e - 26 + In — + o(|a;|°) 
\x\ 


(7.21) 


In order to evaluate I< ( 6 , x) for large negative x, 
note that according to Eq. (4.20) ^i{z/L\mL) = 
^J\m\L^i{\m\z-,\m\L,—l) approaches the half-bound 
state ipo{\m\z) in the limit L —)• 00 . Hence fi( 3 ,x) must 
behave as 


fi( 3 ,x) = (/io(z)+o(|x|°) 

x—^ — oo 

and I^{b,x) becomes 

= dz(/jg(z)-ho(l/|x|). 


(7.22) 


(7.23) 


The integral on the right-hand side has been considered 
in Eq. (6.33). Its behavior for large 6 is given by the last 


line of this equation with the replacement z —>■ 6 : 


dzipl{z) = J^ + b- - \nb + o(l/ 6 ). 


(7.24) 


Upon inserting the above results for Iy{b,x) and 


I<^{b,x) into Eq. (7.18), one sees that the 0(l/|x|) terms 
that diverge linearly or logarithmically as 6 —>■ 00 cancel. 
Hence we can take the limit 6 —> 00 m- Matching the 
terms of order |x|. In |x| and |x|° then yields 

InEi(x) = —|x|-|-ln |x|— 7 e— 2 J(.—ln4-|-o(l). (7.25) 

x—^ — oo 


With the aid of Eq. (6.43) our result for Ei{x) thus be¬ 
comes 

lni?i(x) = —|x|-I-In |x|-I-In —-I-o(l/|x|°), 

ai—>■ —00 TT 

Ei{x) = ® |x|e-l"l+°(i/|x|“)^ (7 26) 

x—>• —00 TT 


In order to determine the leading shifts AE^, of the 
eigenval ues E^^i{x) from their x = —00 values E^, given 
in Eq. (17.31), it is natural to use Rayleigh-Schrodinger 


pert urbati on theory in v{i,x). Since we know from 
Eq. (7.26) that Ei{x) remains exponentially small at 
large negative x, the shift 


AE„ = -2 


/. 1/2 

/ d3u(3,x)[f,(3)]' (7.27) 

Jo 


must vanish when v = 1. Upon subtracting 0 = 2Ai?i 
from AEt,yi and remembering that fi = 1 , we obtain 

r o 9'! 

AE^yi{x) = -2j d3u(3,x) |2 - [f^(3)] | 


. 1/2 


= -4 


'b/\x\ 


d 3 u( 3 ,x) sin^[ 7 r(j 9 - 1 ) 3 ] +0{l/\x\) 


= 2[7r(j9- 1)]' 


In I 


■o(i/N), 


(7.28) 


where we used the fact that the contribution from the 
omitted integral d^ is 0(l/|x |). T o compute the 
remaining integral, we inserted Eq. (7.17) for 11 ( 3 , x) and 
kept only the terms |x|“^ln|x|. The result given in 
the last line means that the eigenvalues U,y>i(x) behave 
in the limit x —> — 00 as 1981 


E^>l{x) = 71 ^( 1 /- 1 )" 


1-f 


2 In I 


0(l/|x|). (7.29) 


The result has a simple interpretation: the eigenvalues 
Ev>i agree to the given order with those of the free {v = 
0) system with Neumann boundary conditions up to a 
reduction 

L[l-h2|xrMn|x| -f 0(l/|x|)]”^^^ 

= L [1-Ixp^nlxl+0(l/|x|)] (7.30) 

of the thickness of the slab. 


B. 


Free-energy functional for the slab in the 
continuum limit 


In Sec. III! Al we considered the continuum limit of the 
self-consistency equation for the potential in detail. In 
order to investigate the scaling function 0 (x) of the resid¬ 
ual free energy, it is useful to integrate the corresponding 
self-consistency Eq. (3.23) to obtain a free-energy func¬ 


tional whose extremum yields the residual free energy. 

Since the bulk free energy and the surface free en¬ 
ergy involve UV singularities, appropriate subtractions 
must be made to eliminate these singularities. In the 
translation-invariant case of pbc, this can be done by 
subtracting the Taylor expansion in the mass parameter 
m? to the lowest order required to yield UV finite expres¬ 
sions |61j . In our case of fbc this procedure is impractical 
because the self-consistent potential has a nontrivial de¬ 
pendence on the temperature variable m. This entails a 
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nontrivial m-dependence of the eigenenergies E^,, which 
is not known in analytic closed form. We therefore use 
an alternative, more feasible approach. 

Starting from the self-consistency equation for the con¬ 
tinuum theory, Eq. (3.231, we will construct a free-energy 


functional whose value at the maximizing potential V:^{z) 
yields the free energy per area of the slab up to reg¬ 
ular background terms. To achieve this, we multiply 
Eq. (3.23) by a “test function” Sv{z) with support on 
[0,L], integrate z over this interval, and interchange the 
order of integrations. We thus arrive at 


^ Sv{0) + Sv{L) 
mL Sv -f —^ 


dw 


LSv 

2y/uj 


■ <5x(w) 


where 


(7.31) 


(5x(a;) = f dzSv(z)Q(z,z;—uj), (7.32) 

Jo 


while Sv denotes the average 


Sv = f ^ Sv{z). 


10 


L 


(7.33) 


At this stage, it is sufficient to presume continuity of 
Sv{z) on the closed interval [0, L]. Since this implies 
boundedness and integrability, the average (7.33) exists. 


Let us show that Sxiu) is indeed the variation of a func¬ 
tion x(a;) caused by a modification Sv{z) of the potential. 
Note, first, that the Green’s function Q{z,z';E) can be 
written in terms of the left and right regular solutions 
ipi^r of the Schrodinger equation 

- dlifiiAz] E) + [viz) - E]ipi^riz, E)=0 (7.34) 

on [0, L] which are fixed by the boundary conditions 

¥>i{z;E) = v^[l-hO(z)], 

2—>-0+ 


ipiiz;E) = VL- z[l + OiL- z)]. (7.35) 

z^L — 


The Wronskian of these functions. 


AiE) = W[Mz;E),^iiz;E)], 


(7.36) 


is called characteristic function of the Sturm-Liouville 


problem defined by Eqs. (7.34) and (7.35). It is an entire 


function of the complex variable E with simple zeros at 
the eigenvalues of this problem. Further, for 

any complex E not in the spectrum {E^} one has 


g(z,z';E) = - 


(piiz<;E) ifriz>-,E) 

ME) 


(7.37) 


This shou ld be obvious becaus e the right-hand side of 
Eq. (7.37) is a solution to Eq. ( |7.34 ) for z ^ z' and the 
jump condition dzGj z, z'; = 1 as well as the 

boundary conditions (|7.35|) are satisfied. 


The above-mentioned properties of A(£') imply that 

ln[const A(E)] = tvlniE — 'Ey) =trln^“^(E). (7.38) 

Noting that 5Q~^ = —SH = —Sv, we see that the varia¬ 
tion of the left-hand side yields the negative of the expres¬ 
sion on the right-hand side of Eq. (7.32). Equation (7.32) 


can therefore we interpreted as the variation of the func¬ 
tion 


xiM = -ln[7rA(-a;)], 


(7.39) 


where the prefactor tt was chosen to simplify subsequent 
relations. 

We wish to consider potentials u(z) that are even with 
respect to reflections at the midplane z = L/2 and have 
a Laurent expansion about z = 0 with a fixed pole part. 
Accordingly, we require that 

viz) = viL — z) (7.40) 

and that u(z) can be expanded about z = 0 as 
1 4?7T- 

viz) =^ + uo +UiZ +..., (7.41) 

4Z^ TT^Z 

where the principal part of the Laurent series (pole 
terms) is fixed while the remaining nonsingular expan¬ 
sion terms including uq may vary. Note that we use the 
notation Uj for the coefficients of the varying Taylor se¬ 
ries part of the potentials considered to avoid confusion 
with coefficients such as uq of the self-consistent poten¬ 
tials for m = ±1. Variations Sviz) therefore satisfy the 
relations 


SviL) = SviO) = Suq 


(7.42) 


and do not involve pole terms. Furthermore, the function 
A (A) simplifies to 


AiE)=dzM{z-,E) 


z=L/2’ 


(7.43) 


In order to integrate xiM with respect to w t o de- 
termine a free-energy functional that yields Eq. (7.31) 


as stationarity condition, we must know the behavior of 
x(w) at large and small w. The leading asymptotic be¬ 
havior of xiM for large a; is a bulk property and easily 
determined. Setting m = 0 and subtracting the bulk 
term at w = 0 to eliminate UV singularities, we obtain 
x(w —>■ oo) ~ — (L/tt) dfcln[(A:^-I-w)/A:^] = 

Since J dzz;(z) has momentum dimension 1, we can con¬ 
clude that the next-to-leading terms are of order 
mod Inw. Furthermore, x(^) can be expressed in terms 
of its derivative as 


nOO 

xiuj) = ~ E\/uj — / dw' 

J ijj 


X'M) + 


2VM_ ■ 


(7.44) 


In Appendix|^we extend the above asymptotic expan¬ 
sion of xiM by computing the two next-to-leading terms 
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by means of the semiclassical expansion, obtaining 


X(a;) = -Ly/uj- 


2y/uj 


4-777 

L7i{[v]-,L, m) H-^ ln(a;L^) 




(7.45) 


where 7^([z;]; L, to) denotes the functional 
7^([^;];L,TO) = ^ + (xe + ln4) + w” 


(7.46) 


while 


_ 2 

u"® = — / dz 

-b ./n 


v{z) 


1 4to 


(7.47) 


is the average of the nonsingular part of the potential. 

In the limit a; —>■ 0, the function x(w) approaches a con¬ 
stant x(0)! whose bulk contribution —Lm follows from a 
calculation analogous to the one performed above. The 
value x*(0) the extremal point Vi^^z] L,m) that maxi¬ 
mizes the functional can be computed in a straight¬ 
forward manner. To this end, we integrate the self- 


consistency equation (3.23) over the interval (0,L), 
terchange the integrations over z and w, and use 


x'(a;) = [ dzg{z,z;-uj) 
Jo 


(7.48) 


to conclude that 


toL -I- 1 = 


dw 


L 

2y/u] 




= -x40). (7.49) 


Here the contribution from the upper integration limit 
vanishes by Eq. (7.44), and the asterisk indicates evalu¬ 


ation at the maximizing potential v^{z; L,m). 

We are now ready to introduce a free-energy functional 
whose variation subject to the constraints specified in 
Eqs. (7.40)-(7.42) is in conformity with the stationarity 
condition (3.23) and gives back the implied Eq. (7.31). To 
integrate we subtract from it a contribu tion w ith 

the same large-w behavior as specified in Eq. (7.45), so 
that the integral of the difference dw converges at the 
upper integration limit. To avoid convergence problems 
at the lower integration limit, we make the replacement 
w —-b w in the last one of the terms explicitly dis¬ 


played in Eq. (7.45), where /i > 0 is an arbitrary mo¬ 


mentum scale, and introduce the free-energy functional 


Tf,{[v]-L,m) = L,m) + L,m) (7.50a) 

with 


= + (7.50b) 

OTT 


and 




X([u],a;) -b Ly/uj 
2m\n(ujL^) L 


ir'^y/uj 

2rn? 

7r2(/r2 _|_ 


2y/uJ 


n{[v\] 1 , 771 ) 


(7.50c) 


where we have explicitly indicated the functional depen¬ 
dence of X on the potential v. 

Since 5TZ(\v\,L,'m) = Sv, the stationarity condition 
(5J^^ ([u], L , m)\v=y^ = 0 is satisfied as a consequence of 
Eq. (7.31), where it should be remembered that only the 
nonsingu ar part of v{z) (including the parameter uq) is 
to be varied whil e the singular one is fixed at vl^{z), as 
indicated in Eq. (7.41). 

The logarithmic derivative of with respect to fi at 
fixed V, L, and to is given by 


2 

TTl^ 

= (7.51) 

Under a change ^ —?> it transforms as 

2 

Tfl^ 

= filL, 

(7.52) 


where the last line follows from dimensional consider¬ 
ations. We denote the value of this functional at the 
maximizing potential u* as 


/^(L,to) = max T^{[v]-,L,m) = ([?;*]; L, to). 

yns—v — v^^ 

(7.53) 

Choosing i = (fiL) ^ in Eq. (7.52) then gives us the 
relation 


TO 


f^{L,m) = L Y{mL) - — ln(/rL) (7.54) 


with 


Y{x) = 


(7.55) 


To relate the scaling function Q{x) of the residual free 
energy to f^{L,m), we must determine the bulk and 
surface free energy contributions Lf\^(rn) and 2/s (to;/!) 
contained in this free energy. The former is readily de¬ 
termined using 7?.([u*]; L, to) = rn?9(m) -b o(l/L) and 
x(a;) = —Ly/w -b rri^9{m) -b o{L). One obtains 

fh{m)= lim L~^f{L,m) 

L—¥00 



(7.56) 
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in accordance with Eq. (2.38). 


To determine the surface contribution to ([u*]; L, m) 
we must compute the limit 

9 

TTi 

2fs{m;fj,)= lim [j'f^{[v^];L,m)-Lf^,{m) + —lii{fiL)]. 

L->-oo TT'^ 

(7.57) 

In Appendix [E] we show that 


=9{m) 


L 


n 2 r 




3 7C(3)1 

m? 


4 TT^ 

2tC’ 


2L\m\ 


m 

2tt^ 


ln(/rL) + 0(1/L). 


(7.58) 


This yields 


/^Ka^)=£3 


3 7C(3) 1^^2|m| 

4 TT^ 2 TT/r 


+ 9{m). 


(7.59) 


The above results enable us to identify those parts of 
the scaling function Y (x) that correspond to bulk and 
surface contributions to the free energy per area. It fol¬ 
lows that the scaling function 0(a;) of the residual free 
energy is related to Y (x) via 


Y{x) = 


(s) 


2 TT 


2x^9{x) H-T 


3 , 7C(3) 

4 TT^ 


-I- Q{x). 


(7.60) 


We are now ready to turn to a discussion of exact 
properties of the scaling functions ©(x) and 'd{x) of 
the residual free energy and the Casimir force. In the 
next subsection, we recapitulate the qualitatively dis¬ 
tinct asymptotic behaviors of these functions in the high- 
temperature, low-temperature and critical-temperature 
limits x —>■ - 1 - 00 , x —>■ —oo, and x —>■ 0±, recall their 
known exact properties, and then focus on the deriva¬ 
tion of their limiting forms for x —>■ —oo. For the latter 
purpose we shall need the following property of the free 
energy introduced in Eq. ( 7.53|), 


i9/^(L,m) _ T77.([u*];L,m)-I-2uom 

dm ~ 87r 7r3 '' 

(7.61) 

This relation may be viewed as the continuum analog of 
Eq. (2.33). Since the term L7^([u»]; L, m) agrees with 
Lw"®, the “trac e” of u"®, up to an L-dependent constant, 
its form (7.61) is plausible. However, its derivation is 


not entirely trivial because of the singular parts of the 
full self-consistent potential u*. From Eq. (7.54) we see 


that what needs to be proven to establish Eq. (7.61) is 
the implied relation for the scaling function Y{x), namely 

Y'{x) = ^TZ{[v^\,l,x) + Cix (7.62) 

87r 


with 


^^0 




I4C(3) 


= 0.054992529830367.... (7.63) 


In Appendix [ f| we give a direct derivation of Eq. (7.62) 
from the free-energy functional (7.50). For the constant 
Cl we obtain the result 


Cl = 


7e 


In 8 




I 


dM^^ii{2z) 

dn 


K—0 


+ 


Ko{z) (7e +ln4) 


^ dW^,fi{2z) 

2z 


dtt 


K—0 


In z, 


(7.64) 


where dtiMi^^o{z) and 9 kITk,o(. 2) are derivatives of Whit¬ 
taker M and W functions. 

We have not been able to evaluate the integral on the 
right-hand side analytically. However, upon computing 
the corresponding expression for Ci by numerical inte¬ 
gration via Mathematica m we verified that the an¬ 
alytical value given in Eq. (7.63) is reproduced to 13 dig¬ 


its. Moreover, we have shown by independent arguments 
that Cl has the exact analytical value given in Eq. (7.63): 


We have determined the asymptotic x —oo be havior 
of Y{x) in two different ways, one that uses Eq. (7.62) 
with an unspecified constant Ci , and a second one based 
on inverse scattering results, which involves the poten¬ 
tial parameter vq (see Appendix [G|). Matching the corre¬ 
sponding results yields the analytical value Ci = vofAir. 

Let us also note the following useful relation between 
the scaling function Y{x) and the value of the functional 
l,x) at the self-consistent potential v^iz] l^x): 


^d^^ - r(a;) = l,x) + (7.65) 

ax TT'^ 


which follows immediately from Eqs. (|7.50|), (|7.55|), 
(lA^, and (lA^. 


C. 


Asymptotic behaviors of the scaling functions 

0(x) and i}(x) 


The scaling functions 0(x) and 'd(x) are known in nu¬ 
merical form from the results of [52] and [55] . They ex¬ 
hibit qualitatively distinct asymptotic behaviors in the 
high-temperature, low-temperature, and critical limits 
X —)■ -l-oo, X —)■ —00 and x —)■ 0±, respectively. Although 
some discussion of these issue can be found in [52|, [55] , 
[38] . and elsewhere (e.g., in [16][ST][24] for general values 
of d), it will be helpful to recall the necessary background 
and some known results. 

For a: > 0, we are dealing with a massive theory con¬ 
fined to a film with free surfaces. Therefore, both scaling 
functions must vanish exponentially; one has 


©(cc —>• oo) 
^{x —>■ oo) 


3—2a;+0(ln x) 


(7.66) 
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where the factor 2 in the exponential is specific to free 
boundary conditions [99]. 

In the critical limit a: —>■ 0, 0(x) approaches the 
Casimir amplitude Ac [cf. Eq. (2.13)], whose exact ana¬ 
lytical value is not known but for which the numerically 
precise result 


Ac = -0.01077340685024782(1) 


(7.67) 


was obtained in [33] . 

Since the free energy must be regular in m when 
L < oo, the singular bulk and surface free energy contri¬ 
butions must be canceled by corresponding contributions 
~ x^, ^ x^, and ^ x"^ In |a;| of 0(a:). As discussed in [38] . 
this implies the limiting behavior 


- 

x—>-0± 

1 


AA, 


(s) 


4871- 


2a;^ 9{x) 


+ X In |a:| -I- Qix + Q 2 X + Q 3 X 


2 Tr^ 

+ o(x^), 


(7.68) 


where AAq is the universal amplitude difference (|6.22 ). 


Here ^jX^ are regular contributions and the o{x^) 

terms involve regular contributions of orders x'^ and 
higher. Setting d = 3 in Eq. (2.12) to determine the 
consequences for ^(a;) yields 


^(x) = 2Ac + Qix- ^ 

s—>-0± Z'K'^ 

- 0337^ -I- o(x^), 


247r 


9(x) 


(7.69) 


from which one recovers the previously mentioned result 


d"(0) = -I/tt^. 


(7.70) 


In the low-temperature limit x —> — 00 , 0(a;) 

approaches the exactly known Casimir amplitude 
—C(3)/(167r) of a free field theory subject to Neumann 
boundary conditions [16] . as follows from the mapping 
to the nonlinear sigma model 


The results given in Eqs. (7.29) and (7.26) can be com¬ 


bined in a straightforward fashion with the known result 
—C(3)/(167rL^) for the residual free energy of a massless 
free field theory subject to Neumann boundary condi¬ 
tions to gain information about the asymptotic behavior 
of the scaling function 0 ( 37 ) for x —)■ — 00 . The reduction 


of L specified in Eq. (7.30) gives us the leading correction 


37 ]“^ In |37| to the limiting value 0(— 00 ). Upon includ¬ 
ing also the next-to-leading term ~ di/|x|, we arrive at 
the limiting form 


0(^) = ~ 
x^-oo iOTT 


^ + 21 n \x\ . ,, 

1 +-^^+o(l/|a7|) 


(7.71) 

This limiting behavior has previously been stated in 
[32] and [33]. The exact value of 2 of the coefficient in 
front of the 1 37 ]“^ In 1 37 ] — here obtained as a consequence 


of Eqs. ( |7.29 | and (7.30) — was derived there by match¬ 
ing the low-temperature behavior of our model to that 
of a nonlinear sigma model. Unfortunately, the latter 
analysis left the exact value of the universal coefficient 
di undetermined. 

We now turn to the direct derivation of the asymp¬ 
totic form (7.71) within the framework of the n —> 00 


analysis of our 0{n) model. For this purpose we must 
extend our results for the free energy and its scaling func¬ 
tion Y{x) given in Eqs. (7.58) and ( |7.60 ) by working out 
their asymptotic contributions in the limit x —>■ — 00 up 
to 0(l/|x|). Let us start by determining those of these 
contributions that originate from the free-energy term 
l,x). We need the asymptotic x —)■ — 00 be¬ 
havior of 77.*(l,x) to o( l/|x|). This can most easily be 
comp uted via Eq. (7.62). Upon inserting the x —)■ —00 
form (7.711 of 0(x) into Eq. ( |7.60 ), we can compute Y'{x ) 
and solve for 77.*(l,x). Taking into account Eq. ( fL^ , 
we find 


7^*(l,x)=-^[l-ln(2|x|/7r)] 

TT^ 


+ C(3) 


1 — di/2 — In jxj 


+ o(l/|xn, (7.72) 


i.e., no terms of order l/|x| appear in the asymptotic 
X —>■ — 00 form of this quantity. 

The calculation of jF^\[v^]-,l,x) to o(l/|x|) is more 
involved and relegated to Appendix [G] It c onfirms the 
asymptotic form of 0(x) given in Eq. (7.68) and yields 
for the coefficient c?i the result 


7 e 


In 


- 1-2 


C'(3) 
C(3) 

= 0.96720564466060109.... 


di = 2 


(7.73) 

This value is consistent with the result di = 1.0(1) 
obtained in [33] by numerical solutions of the self- 
consistency equations. 

The X —>■ — 00 behavior of the Casimir force scaling 
function implied by Eq. (2.12) is, of course, the same as 
given in [33] . namely 


d{x) = ~ 

x^-00 » 7 r 


1 -b 


3di/2-l-t31n| 


■o(l/|x|) , 
(7.74) 


except that we now have the exact analytical value (7.73) 
for di. 


VIII. SUMMARY AND CONCLUSIONS 

In this paper we explored the potential of inverse scat¬ 
tering theory to gain exact information about the n = 00 
solution of the 0 {n) cp'^ on a three-dimensional film of size 
00 ^ X L bounded by a pair of free surfaces. A wealth of ex¬ 
act results could be determined. Since their derivation re¬ 
quired a combination of various tools and lengthy calcu¬ 
lations (most of which are presented in the appendixes), 
it will be helpful to briefly summarize them. 
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The essence of our strategy was to eliminate the poten¬ 
tial v(z) of the self-consistent Schrodinger problem that 
the exact n = oo solution involves in favor of scattering 
data. This enabled us to reformulate the self-consistency 
equation in terms of scattering data. 

For the semi-infinite case L = oo we succeeded in the 
exact determination of the scattering data pertaining to 
the self-consistent potential v(z', L = oo, m) for all values 
of the temperature variable m at, above and below criti¬ 
cality ( m = 0 ). The co rrespondi ng re s ults ar e pre sented 
in Secs. IVB and IVC see Eqs. (4.68), (4.41) and ( |4.54 |, 
and (4.64) and (4.66) for the cases m = 0, m > 0, and 


directly. Our results for the latter are given in Eqs. (7.71) 
and (|7.74|). They agree with the asymptotic forms ob- 


m < 0, respectively. The obtained phase shifts r]o{k) and 
r]±{k) are depicted in Eig.|^ 

The knowledge of these scattering data enabled us to 
obtain exact results for the two-point boundary correla¬ 
tion function for all m ^ 0. The results are described in 
Secs. IV BTIV DI and summarized in Table HI 

We then exploited these scattering data in conjunc¬ 
tion with information about the potential obtained via 
boundary-operator expansions, on the one hand, and var¬ 
ious semiclassical expansions and a trace formula, on 
the other hand, to derive a variety of exact results for 
the case of a film of finite thickness L. In Sec. EH we 
computed several quantites related to the self-consistent 
potential, among them the potential coefficient vq = 
m~^ lim 2 _,.o oo, to), the universal amplitude differ¬ 
ence and an excess quantity associated with the 

squared order-parameter density. The results are given 
in Eqs. (4.7), (6.22), and (6.43), respectively. 


A major issue to which we then turned was the asymp¬ 
totic behavior of the spectrum and the scaling func¬ 
tions of the residual free energy and the Casimir force 
in the low-temperature scaling limit x —>■ — oo. Our 
findings for the asymptotic behavior of the spectrum 
in this limit are presented in Sec. |VII A| Our result 
for the exponential vanishing of the lowest eigenvalue 

The remaining energy lev- 
imiting values associated with 


El is given in Eq. (7.26). 
els Ej^>i approach their 


a free massless theory subject to Neumann boundary 
conditions algebraically modulo logarithmic corrections; 
see Eq. (7.29) and m- We also obtained the asymp¬ 
totic behavior of the self-consistent potential in the inner 
region of the film; our results for the scaled potential 
v{zfL] l,a; < 0) are presented in Eqs. (7.16) and (7.17). 


To derive the asymptotic behavior of the scal¬ 
ing functions 0(a;) and we extended our anal¬ 

ysis by constructing a free-energy functional (cf. 
Sec. VIIBI whose stationarity condition yields the self- 


consistent Schrodinger equation involving the potential 
(see Sec. VIIB|. This enabled us to obtain an integral 
representation for the scaling function Y{x) of the free- 
energy density in terms of the Green’s function of the 
continuum theory and the self-consistent potential. 

Using semiclassical expansions in conjunction with the 
L = oo scattering data and our other results, we then 
determined the asymptotic x —> —oo behavior of this 
function Y{x), from which those of 0(a:) and -d^x) follow 


tained in |32| and |33) and corroborated there by numer¬ 
ical results. 

The analytic results of the present work for the asymp¬ 
totic X —> —00 behavior of the scaling functions go in 
several ways beyond those of |M] and [331: In the latter 


two papers, both the asymptotic forms (7.71) and (7.74) 


and the exact values of the amplitudes of the logarithmic 
anomalies ^ |x|“^ln|a:| of 0(a:) and 'd{x) were obtained 
via the mapping to a nonlinear sigma model. Here, we de¬ 
rived these results directly from our 0 {n —>■ oo) model 
in the scaling limit. In addition, we calculated the exact 
analytical value of the coefficient di which governs the 
amplitudes of the x —)■ —oo contrib ution s oc of 0(a;) 


and ^{x). This value, given in Eq. (7.73), is in conformity 


with the numerical estimate of [33 . 

Note that in our analysis in this paper we started from 
a lattice model to avoid UV singularities, but then fo¬ 
cused quickly on its continuum scaling limit (specihed 
at the beginning of Sec. Ill A). The rationale of this ap¬ 


proach is easily understood: Of primary interest are the 
universal features encoded in the scaling functions Q{x) 
and {^{x). Focusing on the scaling limit is a convenient 
way to eliminate the less important microscopic details 
associated with the behavior on microscopic scales. If 
one determines the self-consistent potential and the free 
energy of the n —)■ oo solution of a given 0 {n) lattice 
model representing the universality class of our model by 
numerical means, as was done in |32] and |33j for two 
families of such models, care must be exercised when ex¬ 
tracting the universal scaling functions 0(a:) and ^{x) 
from the numerical data. The latter include the full in¬ 
formation about the nonuniversal behavior of the model 
on microscopic scales. The determination of the scaling 
functions requires the elimination of such nonuniversal 
small-scale features. To this end, the appropriate scaling 
limit TO —)■ 0 and L —> oo at fixed x = uiL must be inves¬ 
tigated, where corrections to scaling should be properly 
taken into account. As x takes on large negative values, 
it becomes increasingly more difficult to determine reli¬ 
able results from the numerical data. Nevertheless, suffi¬ 
ciently large values of —x could be reached in [33] to see 
the asymptotic x —> —oo behavior. Moreover, the numer¬ 
ical results of [33] for the scaling functions 0(a;) and 'd(a;) 
are consistent with all exactly known results, including 
the value of the coefficient di obtained here. The preci¬ 
sion with which the latter has been numerically checked 
is admittedly still modest. However, other exactly known 
properties such as ^?"(0), the universal amplitude differ¬ 
ence or the amplitude of the |a:|“^ln|a:| contri¬ 

bution to the asymptotic x ^ oo behavior of 0(a;) have 
been confirmed with much greater numerical accuracy. 

In a recent paper [lOOj . Dantchev et al (DBR) stud¬ 
ied the exact solution of a mean spherical lattice model 
on a film of size oo^ x L bounded by two free surfaces 
with separate constraints imposed on the average of the 
squares on the spins in each layer 2 . This model is known 
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to describe the n —>■ oo limit of a fixed-length spin 0{n) 
spin model [lOlj . DBR investigated, on the one hand, its 
exact solution by numerical means and derived, on the 
other hand, a low-temperature expansion for finite L. A 
discussion of how their work relates to the one of Diehl 
et al (DGHHRS), presented first in the earlier letter (32] 
and subsequently in more details in [33] , has already been 
given in m- Hence we restrict ourselves here to a few 
clarifying remarks. 


DBR overlooked the fact that 0{n) (j)'^ models with 
quadratic and quartic interaction constants r and g turn 
into fixed spin-length models in the limit g —> oo at fixed 
T/g. Therefore they did not realize that DGHHRS’s work 
includes the exact solution of their mean spherical model 
as special case and that the numerical data DGHHRS 
presented for the g —>■ oo limit of DGHHRS’s model 
B directly apply to this mean spherical model. 

DBR’s numerical data appear to be in conformity with 
those of DGHHRS. However, differences occur in the 
analysis of the numerical data. In the near-critical regime 
(small to) DBR recovered the scaling behavior reported 
by DGHHRS. On the other hand, they had problems to 
observe scaling for larger values of —x = —mL. As is 
discussed in |102j . this may be attributed to their inclu¬ 
sion of data well outside the scaling regime in attempted 
scaling plots and an insufficient study of the scaling limit. 
In their low-temperature analysis, DBR focused on the 
behavior for low temperatures at fixed L rather than on 
a proper determination of the scaling functions for large 
values of —x. They found that their low-temperature 
data at fixed L are well described by the analytical low- 
temperature expansion they determined. 

Note that the asymptotic x —>■ —oo behavior of the 
scaling functions 0(a;) must of course comply with the 
low-temperature expansion. Inspection of this expansion 
by DBR shows indeed that the limiting value 0(—oo) 
and the coefficient of the |a;|“^ln|a;| contribution [cf. 
Eq. (7.71) and the related one for •d{x), Eq. (7.74)] may 
be recovered from it. On the other hand, the coefficient 
di of the correction term oc cannot be derived from 
it. Unfortunately, the information that can be extracted 
from the low-temperature expansion is rather limited. 
The reasons are twofold. First of all, the expansion is 
not in inverse powers of the scaling variable x (modulo 
logarithms). Second, nonuniversal microscopic contribu¬ 
tions are still contained in this expansion. 

As can be seen from the lengthy and, in part, tedious 
calculations we had to perform, the precise determina¬ 
tion of the asymptotic x —>■ —oo forms of the scaling 
functions to the given order in l/|a:| is a quite demand¬ 
ing task. By taking the continuum limit, we could elim¬ 
inate the nonuniversal behavior at microscopic distances 
from the boundaries. For local densities such as the en¬ 
ergy and magnetization densities distinct nontrivial be¬ 
haviors in the near-boundary regimes of distances L — z 
and z < l/|w| <?; L from the surfaces and the inner 
regime 1 /|to| <C z, A — z < L/2 occur, both of which gen¬ 
erally contribute to quantities involving integrals over the 


film such as the free energy or the Gasimir force. 

Although we have not been able to determine the de¬ 
sired scaling functions 0(x) and d{x) exactly in closed 
analytic form, the wealth of exact analytic results we have 
been able to obtain by means of inverse scattering theory 
indicates that the strategy we have pursuit — eliminating 
the potential from the self-consistency equation in favor 
of scattering data — has considerable potential for suc¬ 
cessful applications to analogous self-consistency prob¬ 
lems. Interesting examples of such problems where simi¬ 
lar strategies might prove fruitful are studies of stochas¬ 
tic dynamic models and quantum quenches in the large-n 
limit. 
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Appendix A: Properties of the functions W 2 and U 2 


In this appendix we give the explicit analytic expres¬ 
sions for the functions IU 2 (A) and 172(A) in terms of spe¬ 
cial functions and briefly discuss some of their properties. 
The function H 2 (A) can be expressed exactly as 


IU2(A) = 


K 


(Al) 


7 r(A-|-4) \^A-|-4^ 

in terms of the complete elliptic integral of the first kind 

da; 


K{\) = r 

^0 


0 a/( 1 - a;2)(l - A^ a;2) 
1 1 


F 


2 ’ 2’^'’^ 


(A2) 


where pFq denotes the generalized hypergeometric func¬ 
tion. 


and 


]) 


Its antiderivative [/ 2 (A) defined in Eq. ( 2.25 ) can be 
determined in closed form by integrating Eq. (A2). This 
gives (cf. Eq. (48) of 


C/2 (A) = 


-2 


(A + 4)2 
-I- ln(A -I- 4). 






(A3) 


From the definition of IU 2 (A) in Eq. (2.17) one can read 


off that it can be analytically continued into the complex 
A-plane except for the branch cut [—8,0]. To determine 
the form of W 2 (A) for small A, we use the expansion 

r2a 


= Q) 

j=0 '' L 


[ln{l/k)+dij)] (A4) 
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of K(k) near the singularity at fc = 1 (see Eq. 19.12.1 in 
where 


(A5) 


k = yi - fc2 

is the complementary modulus, while 

{a)j = a{a — l)(a — 2) • • • (a — j + 1) (A6) 

denotes the Pochhammer symbol and d{j) can be ex¬ 
pressed in terms of the digamma function ijj{z) as 




(A7) 


This expansion converg es fo r 0 < |fc| < 1 [MUSS]. Upon 
substituting it into Eq. (Al I and expanding in A, we con¬ 


clude that W 2 (A) can be written as 

fU 2 (A) = —W 2 {X) In A -f i? 2 (A), 


(AS) 


where W 2 (A) and i? 2 (A) are both analytic for sufficiently 
small |A|. The former function, which for A G (—8,0) 
measures the discontinuity across the branch cut, is the 
spectral function introduced in Eq. (2.40). It can be com¬ 


puted for real values of A in a straightforward fashion ei¬ 
ther from Eq. (Al) or the integral representation for W 2 
implied by Eq. (2.17), namely 


-I-4sin^ (pi/2)-I-4sin^ (P 2 / 2 )]. (A9) 

One finds that 


W2{X e M) = ( ^/— 


27r2 


(V-A(A + 8)/16) e{8 + X)0{-X), 

(AlO) 


where K{k) is the elliptic function (A2). 


Integration of Eq. (A8) implies that U 2 {X) can be writ¬ 
ten as 

C/2(A) = MiA^[l + AA(A)] + B(A), (All) 

where A (A) and B{X) are analytic for sufficiently small 
|A| and related to the functions W 2 {X) and R 2 {X) via 


and 


u;2(A) = —[1 + 2AA(A) + A2A'(A)] (A12) 

47r 


i?2(A) = A[yi(A) + AA'(A)] + B'(A). (A13) 

47r 


Their power series expansions 

,,,, 1 5A 7A2 

^ ' 16 768 


169A3 


and 

B{X) = — 

TT 


8192 1310720 


5Aln2 A2(3-101n2) 


47r 


1287r 


O(A^) (A14) 


0(A3) (A15) 


can be determined in a straightforward fashion. Here 

(A16) 


G = V ^ = 0.9159655942 ... 


is Catalan’s constant. The value AG/tt of B{0) = 1 / 2 ( 0 ) 
follows directly from Eq. (A3) by evaluating the right- 


hand side of this equation at A = 0. 


Appendix B: Derivation of equations for the 
scattering phase ri{E) 


We start by giving a direct proof of Eq. (4.26) for the 
derivative drjE/AE of the scattering phase. 

Let Q±{z,z'; E) = {z\{E — 'Hy^) ~^\z') be the Green’s 


function introduced in Eq. ( |4.21b ) for the cases m = 
±1. It is defined for all complex values of E that do not 
belong to the spectrum spec('H^j_) = [(5±ip,oo) and has 
a branch cut along the positive real axis starting at (5±p. 
Eor energies E -|- iO = 4i,i + -|- iO on the upper rim of 

the branch cut, it can be represented as 


g±iz,z';E + i0) = - 


y’(u k)/(z',k) 
F(k) 


for z' > z, (Bl) 


where (p{z, k) and /(z, k) denote respectively the regular 
and dost solutions of the Schrodinger Eq. (4.3), while 
f (k) is the dost funct ion in troduced in Eqs. (4.11) and 
(4.12). Note that Eq. (4.14) for E{k) has been used. 

To relate the numerator of the expression for the diag¬ 
onal element of the Green’s function given on the right- 
hand side of Eq. (IBII) to a Wronskian, we consider these 


regular and dost solutions for two different positive val¬ 
ues k and k'. We multiply the differential Eq. (4.3) for 


ip{z,V!) by /(z, k), that for /(z, k) by (p(z, k'), and then 
subtract the results from each other, obtaining 

(k'2 - k2)^(z, kO/(z, k) = d,W\^{z, k'), /(z, k)]. (B2) 


This can be integrated to yield 


dzip{z, k')/(z, k) = 


W[(^(z,k'),/(z,k)]| 
k'2 - k2 


(B3) 


where z^ > 0 is a finite upper integration limit whose 
limit z„ —>■ (X) we intend to take. 

We first evaluate the Wronskian at the lower limit 
0-I-. In addition to the regular solution (/?(z, k), the 
Schrodinger Eq. ( |4.3 ) has a linearly independent second 
solution, which we denote as p(z, k). Unlike the small- 
z behavior (4.9) of (/?(z, k), the latter function varies as 
£i(z, k) ~ z^T^Tnz for z —>■ 0. More precisely, insertion of 
the ansatz 


g{z,k) = (p(z, k) lnz-t-p(z, k) 


(B4) 


into the differential Eq. (4.3) shows that 

1 " 


[-8^ + v±- k^]p(z, k) = 


z 


V5(z,k) = 0(z 

(B5) 
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which in turn implies that 


p(z, k) = (^(z, k) Inz + 0(z^/^). (B6) 

z-fO 

We can now expand the Jost solution as 

/(z, k) = A ¥>( 2 , k) + Bk p(z, k) (B7) 

and compute the limit 

lim W[¥>(z, k'), /(z, k)] = lim W[(p{z, V!),q{z, k)] 

z—>-0+ z—>-0+ 


= Bk lim W[^/z,^/zhlz] = = —F{k), 

z—>-0+ 


(B8) 


where the last equality follows from the fact that the 
result is independent of V! and Eq. (4.14). 

To compute the Wronskian for large z^, we use 
Eq. (4.11) for (fi{z,k'), substitute the Jost solut ions 
/(z, ±k') and /(z, k) by their asymptotic forms (4.8), in¬ 
sert the resulting value of the Wronskian at z„ together 
with Eq. (B8) into the right-hand side of Eq. (B3), and 
then take the limit k' k. We thus arrive at 


dz (z, z; E-I-iO) = 


2kF'{k) + F{-k) e 


1 

4k2 
— 1 — 2ikz,, 


2 ikzi, 


F(k) 


(B9) 


The (oscillating) term cx vanishes as z„ —)■ oo 

owing to the positive imaginary part of k implied by the 
iO of the energy. Proceeding to the limit z„ —^ oo yields: 


dz 


^±(z,z;£’-ki0) -k ^ 


1 d 

■inF(k)' 

“ 2kdk 

ki/2 


(BIO) 

For real positive k, the term oc (InE)'(k) becomes 
\u'{k) — i77'(k)]/(2k). Taking the imaginary part of the 
above relation finally gives 


Im 


dz 


g±{z,z]E + \Q) F 


rm 

2k ’ 


(Bll) 


which is Eq. (4.26). 

We next turn to the proof of Eq. (4.27). Let us con¬ 
sider cl^ges of the potential v±{z) —>• v±{z) + 5v{z) in 
Eq. (4.26), where the var iation 5v(z) has the properties 
(i)-(iv) specified in Eq. (4.23). For the implied linear 
variation Stje of the scattering phase we obtain from 
Eq. 


dSrjE 

dE 


pOO 

/ dzlm5^(z, z; + iO) 

Jo 

{E + iO- Hy^y^Sv {E + iO- 


= —Imtr 


= dE'"*' 


Sv {E + io-nyy 


-1 


(B12) 


This can be in tegra ted with respect to E from 0 to E 
to obtain Eq. (4.27). The integration constant 6 r]\E=o 
is zero because condition (iv) of Eq. (4.23) ensures that 


the values r7_E|£;=o = 0 and ?7 _e|£;=o-i- = 7’'/2 for the cases 
m = 1 and m = — 1 are not changed by dv(z). 

Note that the above analysis can be extended to 
all complex k in the upper half-plane. In particular, 
Eq. (BIOI holds for all complex k with 0 < argk < tt. 


Appendix C: Semiclassical expansions 


1. Half-line case 


Here we use semiclassical expansions to gain informa¬ 
tion about the asymptotic k —>■ oo behaviors of the Jost 
solution, the regular solution, and the Jost function of the 
Schrodinger problem on the half-line specified in Theo¬ 
rem |VLl] 

Following a standard procedure (described on pages 
35ff of |103] 1. we make the ansatze 

/(z, k) = 

ifiz, k) = yl(k) + cc (Cl) 

for the Jost and regular solutions, where “cc” means 
complex conjugate. They lead to a Riccati equation for 
i9z3^(z, k), namely, 

2 kaz 3 ^(z, k) + [d,y(z, ky-id^y{z, k) + v{z) = o. (C 2 ) 

The function 3^(z, k) has for k —>■ ck) the asymptotic ex¬ 
pansion 


3^(2,k) = '^yyz)k F 


(C3) 


i=i 


The coefficients 3^j(z) follow from the equations that re¬ 
sult upon iteration of Eq. (C2). For and 3^2 one 


finds the differential equations y'i{z) = —v{z)l2 and 
3 ^ 2 ( 2 ) = i3^"(z)/2, respectively. Solving these with the 
boundary conditions 3 ^ 1 ( 00 ) = 3 ^ 2 ( 00 ) = 0 gives 


3 ^ 1 ( 2 ) = ^^ dsz;(s), 
iv{z) 


3 ^ 2 ( 2 ) = - 


(C4) 

(C5) 


The large-k expansion (C3) holds in the regime z ^ 


1/k, but cannot be used in the narrow boundary layer 
0 < z < 1/k. For the latter, we will give an alternative 
large-k expansion below. Before we turn to it, let us first 
investigate the behavior of 3^i(z) for z <C 1. 

Owing to the singular near-boundary behavior of the 
potential [see Eqs. (6.1) and (6.2)], divergences linear in 
1/z and logarithmic divergences appear in the asymptotic 
form of 3^1 (z) for z —>■ 0. To uncover these divergences, 
we add and subtract appropriate contributions to 3^i (z), 
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rewriting it as 


where W(<;, k) behaves as 




1 


V-1 


Inz 


= +0(z^) (C6) 


with 


a = dz 
Jo 




(C7) 


We thus arrive at the following two-parameter asymp¬ 
totic expansion 


ln/(z,k) = ikz- — 


^ -I- v-i Inz - a -I- uoz -t- O(z^) 


1 r 1 

4k2 422 




- Uo -k 0(z) 


0(k-3), 

(C8) 


which holds for k —oo with 1/k <C z 1. 

In order to develop an alternative large-k expansion 
that holds in the near-surface regime, we express the reg¬ 
ular solution in terms of the variable c = zk, defining 


k) = yk(^(^/k, k), 


(C9) 


and expand z;(<j/k) in inverse powers of k. We thus arrive 
at the differential equation 




‘^(k, k) = 0, 

(CIO) 


k) = - — - In - iM. In . + 2 


m T^k 
' 4 2 




+i 

+ + 0 {<^ 


(C15) 


for 2rkS^ 3> 1. The result can be expanded in inverse 
powers of k. This gives 




+o(v“) + ^y 


17E 

2 

1 


iln(80 1_ JL+o(,-3) 

2 4 ^ 32^2 + ) 


ivo , Wo , 

-T^ + T + ^ 

1 /I 


+ jE (2.i,- |1) 


16c 
0(c-^ 


+ 0(k-"). 


(C16) 


Returning to the variable z = ck and the function (/j(z, k), 
we arrive at the result 


ip{z,k) = 


1 


V^k/2 

with 

W(z,k) = lT(zk, k) 

iTT 1 f —i 

= ikz-h - - 

4 k \ 8z 




(C17) 


W-l 


17E 

2 


TT 

4 


is to be solved subject to the boundary condition i ln(8zk) 

izwo' 


r-1 

+ 2 

2 J 

1 ^ 4k2 1 

4z2 z 


(p(c ^ 0, k) = v^[l -kO(c)]. 


By analogy with Eq. ( |6.7| ), the solution is again given 
by a Whittaker M-function; one has 

(^(c,k) = ?^Pi^M_i^„o(-2irkc)[l+0(k-3)], (Cll) 
V2wk 

with 


■ Wo 


-kO(k-^lnk). 


(CIS) 


In the regime 1/k <C z <C 1 (i.e., 1 <C C <C k) both this 
expansion as well as the one implied by Eqs. (C8) and 


(Cl) for the regular solution hold. Hence we can match 


them using Eq. (4.11). Upon substituting our result (C8) 


and 


a/ 1 - Wo/k2 

(C12) 

W-l 

' “ 2krk' 

(C13) 


Using the asymptotic expansion of the Whittaker M- 
function given in Eq. (13.19.2) of [51] yields 


for the Jost solution into this equation a nd eq uating it 
to Eq. (C17) with lU(z, k) given by Eq. (C18), one can 
determine the functions F{k), a{k), and the phase shift 
r]{k) in a straightforward fashion. One obtains 


^(c,k) = ^e’^C,k)+cc, 

V ZTT 


(C14) ^(k) = iln(2k/7r) + ^ + ^^+0(k-3), (C20) 
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VW = i + ln(8k)] + 0(k-3), (C21) 


and 


Proceeding as in the previous subsection, one finds that 
the result for the regular solution in the region 1/k ^ z < 
L\m\/2 can be written as 


In the m = — 1 case, where cr(k) = CT_(k) and V-i = vZi, 
the last equation becomes Eq. (4.62). 


2. The finite interval case 

An analogous calculation can be made for the singular 
Sturm-Liouville problem 

-a^(^(z, K) + v{z)ip{z, K) = kl ip{z, k^), 1 / € N, 

on the finite interval 0 < z < |m|L with an even potential 
v{z) = v{L\m\ — z) that has a Laurent expansion of the 
form 

u(z) = -^ + —+U 0 + ..., 

near the boundary plane z = 0. The usual boundary 
conditions for the eigenfunctions (/^(z, k^) are implied: 

ip(z,k„)= (C23) 

/Vz[l + 0(z)] forz^O, 

I [1 + 0(L|m| — z)] for z —>• L|m|. 


While the result (C5) for 3^2 (z) remains valid except that 
z now is restricted to [0 ,L|to|], the analog of Eqs. (C4) 


and (C6) becomes 

f*L|m|/2 


1 

yiY)= 2 j dzz;(z) 


Llml— 1 


4 " ' 4L|to| 8z 

+ 0{zY, 


1 v-i 2z vqz 

2 L|m| 2 


(C24) 


where 


^L|m|/2 


L|m| 


dz 


K^) 


1 

4z2 


f-1 

z 


Instead of Eq.(C8), we now have 


In /(z, k) = ikz + ^ 


1 

4z 


2L[ 

I 

4k2 


+ v-i In , + uqz + O(z^) 


L\m\ 


1 

4z2 


V-1 


- Uo + 0(z) 


+ 0(k-3) 
(C25) 


^) = \l cos[Q(z, k)] (C26) 


with 


k) - ^ + ^ + ^ + 0(k 3), (C27) 

Q{z, k) = kz — Y — TP ( -pi—r + L|m| u" 


4 4k 
~ + ln(4kL|m|)] 

1 /•^l”i|/2 

+ ^J dCu(C) + 0(k-3lnk). (C28) 

The discrete values of k^ corresponding to the eigenen- 
ergies = k^ follow from the Neumann and Dirich- 
let boundary conditions that hold at the midplane z = 
L|m|/2 for odd and even values of v = 1 , 2 ,..., oo, re¬ 
spectively. Noting that u'(L|m|/2) = 0, one concludes 
that the k^ for all G N are given by the solutions to the 
equation 


Q{L\m\/2, kY = 7r(i/ - l)/2. 
Solving this in an iterative manner yields 
1 r 1 _ 2u_i 


(C29) 


k;, = kl, H-^ 


2k^ {{L\m\Y L\m\ 

-I- 7e + In4] I -I- 0(k“^lnki,), 


[ln(ki,L|m|) 

(C30) 


where kj, = Tr{v — l/2)/(L|m|). 


Appendix D: Asymptotic behavior of xiY 

In order to derive the asymptotic large-w behavior of 
x(a;) asserted in Eq. (7.45), we use again the semiclassi- 
cal expansion, proceeding along lines analogous to those 
followed in Appendix 0 Since we are now dealing with 
the case of negative energy E, we write q = Y—E and 
make the ansatz 




(Dl) 


for the left regular solution, analogous to Eq. (Cl). The 


implied Riccati equation for dz2(z, q) becomes 

-2qdzZ{z,q)-[dzZ{z,q)]‘^-d^^Z{z,q)+v{z) = 0. (D2) 

We expand Z{z,q) in inverse powers of q and separate 
the contributions from even and odd powers, writing 


for the logarithm of the Jost solution. 


Ziz,q) = Zo{z,q) + ZYz,q) 


(D3) 
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with 


subject to the boundary condition 


2o{z,q) ='^b2j+i{z)q 
i=o 

OO 

2 e{z,q) ='^b2j{z)q~'^^. 
i=i 

For the two lowest-order coefficients one finds 

rL/2 


1 f ' 

bi{z) = --J Az'v{z') 


and 


b 2 {z) = -v(z)l4:. 

The function Zg{z,q) must satisfy the equation 
2d:,Ze{z,q)[q + d:,Zo{z,q)] = -dlZo{z,q). 


(D4) 


(D5) 


(D 6 ) 


(D7) 


This can be integrated in a straightforward fashion to 
obtain 

Ze{z,q) = + d^Zo{z,q)/q], (D 8 ) 

where the integration constant has been chosen such that 
the consistency of the result with the expansions (D4) 
and Eqs. (D5) and (D 61 is ensured. 


Inserting Eq. (D 8 ) into Eq. (D3) and the resulting ex¬ 
pression in turn into the ansatz then gives 

^liz, q) = A, q^/^[q + q)]-^'^ (D9) 

To determine the consequences for the functions 


A{E = —q^) and we use Eq. (7.43), taking into ac¬ 
count that Zi{L/2,q) and q)\z^L /2 vanish since 

Zi{z,q) is an odd function with respect to reflections 
about the midpoint i/2, i.e., Zi{z,q) = —Zi{L — z,q). 
We thus arrive at 


A{E) = 2qAle‘^^\E=-g2, 
Xi.q^) = -qL - \TL[ 2 TTqAl]. 


(DIO) 

(Dll) 


Hence the problem is reduced to the calculation of the 
large-g asymptotics of the amplitude Aq, which is fixed 
by the boundary condition (7.35) at z = 0 for the left 
regular solution ipi. 

As before, the above semiclassical expansion cannot 
be used in the near-boundary region with 0 < z < 1 /q 
because of the singularity of the potential at z = 0. By 
analogy with the analysis in Appendix the remedy is 
to go over to the rescaled variable q = zq, to introduce 


‘fl{^,q) = y^^liz,q), 


(D12) 


expand v{(;/q) in inverse powers of q, and solve the re¬ 
sulting differential equation 


„2 1 , 4to 

— d — —^ -|- 1 H- q — 


vq m 

„2 


+ 0 {q-^) 


‘fi{^,q) = 0 

(D13) 


^ii<i,q) =V‘^[i + 0(>r)]. 

^->■0 


(D14) 


The solution is 


<fi{<;,q )=2 i/ 2 M_ 2 ^^_ 2 ^-i_o( 2 cK 9 / 9 )[l-kO (9 ^)] 

(D15) 

with 


= ^/q^ + Vom'^. 


(D16) 


Using the asymptotic expansion (Eq. 13.9.1 of [M] l 




e-^x-^ ^ ( 1/2 - k)j{1/2 -H k)j 


x-s-oo r(l/2 — k ) 


E' 

3=0 


j\x^ 


(D17) 


of the Whittaker function, where {a)j is the Pochham- 
mer symbol introduced in Eq. (A 6 ), and expanding sub¬ 


sequently K = — 2 rmT ^/Kq in powers of 1 /g, one is led 
to the following two-parameter expansion for large and 

g: 


Hipi{q, q)]=-^ ln(27r) + ^ + ^ + ^ + 0{q 


-3\ 


-k 


2 m 


ln(2<^) -I- 7 e -I- ln4 - ^ -I- 0(c 


1 

H—^ 


T 

+ 0 {q-^ 


+ 0 ( 1/0 

4 TT^ 


(D18) 


In the crossover regime l/g <C z <C 1 both Eq. ( |D18 ) 
and the original semiclassical expansion are valid. In this 
regime, the coefficient 6 i(z) is given by 

1 1 O 1/*-^/^ 




Vq m?Z 


rL /2 


dz' 




I 


8 z 4L 


+ ^ln(2z/L) + 0(z2). 


(DI9) 


Upon substituting this small-z expansion along with that 
of 62 (-z) into Eq. dm^ , we arrive at the expansion 

\ii[ipi{q/q,q)] = In \q^^^ifi{q,q) 


— In(Aq) -|- 


1 


1 


1 


1 2m Am 
2L 


4m , 2c ,, 

+ ^ In — - / dz'u“(z') 

qL Jo 


+ 0(c-^) 

+ 0 {q-^) 


1 


uom2(2C - 1) 


+ o(k-^) 


(D20) 
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in the crossover regime. This can be matched with the 
two-parameter expansion (D18) to conclude that 


1 


1 


8m . 


\n.Aq=- - ln(27rg) -|- — L, m) H-^ In(gL)] 


m 

Tv'^q^ 


Aq 
0 {q-^) 


(D21) 


where 7?.{[u*]; L, m) is given by Eq. (7.46). Substitution 
of this result into Eq. (Dll I finally yields the asymptotic 
form stated in Eq. (7.45). 


Appendix E: Calculation of /a(m;/i) 

In this appendix we compute the surface free energy 
/s(m;^) from the free-energy functional (7.50). We be¬ 
gin by comp uting t he contribution ([u»J; L, m) de¬ 
fined in Eq. (7.50b). Upon adding and subtracting to 


the nonsingular part u"® of the self-consistent potential 
Vif the bulk potential Ub(m) = m? 6 (m) and inserting 
uot{L, m) = Vo we arrive at 


=—Vbi'm) + -find) 




mL . 


+ -§7 + 0{1/L). 


(El) 


The integral LuJ® — vo can be decomposed as 
— Wb =2|m|Q;± — 2|m 


Asgn{m) 


— 2|m| / dz 

J\m\L/2 

- Ub(sgn(m)) -f ^ 


u(z; 00 , sgn(m)) 


■0(1/L). (E2) 


In the second integral, we can replace w(z; oo, sgn(m)) — 
Ub(sgnm) by the limiting forms ^ e“^^ and —(2z^)“^ for 
m > 0 and m < 0, respectively, to see that it yields 
contributions ~ Q-mLl 2 ^ L~^. Dropping these, we 
arrive at 

L - Vh = 2|m|[a± =F 47r“^ ln(|TO|L/2)] -b 0{l/L), 

(E3) 

which upon addition of the remaining contributions to 
and inser tion of our results for a± given 
in Eqs. ( 6.20[ ), and (6.21) yields 


L,m) = 6 {m)m’^ 


L^+4A4-> 


Vo m 
Sir 


+ 


1 - In 


2L| 


m 


+ 0(1/L). 


(E4) 


( 

Turning to the calculation of Jyi t[u*];L, m), we first 
consider the case m < 0. In this case one has Vh = 0, 
together with the relation 


LTl^.{L,m) = 


8 m 


) +0(1/L), (E5) 


following directly from the above analysis. After sub¬ 
stitution of (E5) into (7.50c) and transformation of the 
integration variable w = m^y, one obtains 

_ 2 poo 

J dyI{y,mL) 


(E6) 


where 


^(j/;a;) = x( 2 /;T = |a;|,m = -I) -b |a:|v^ (E7) 




I -bln 




Tr^{y -b 1) 


To obtain the contribution of x(y; |a;|,—1) = x(a;;L,m) 
to I{y;x), we have set L = oo in the regular solution 
(b’(z|m|,/c/|m|; L, —1) and used its relation 

yp(z,k;oo,-l) = e^v^[l-b 0(l/z)] (E8) 

z»i Z^y 

to the dost function for large z, obtaining 

X{y,\x\,-'i^) = -\x\^ + \n{2y/y/Tr) 

x—^ — oo 

-21nE_(iy^) + o(|x|°). (E9) 


To compute the integral over I{y;x) in ( |E6[ ), one can 
integrate by parts, using the fact that the boundary 
terms vanish since I{y;x) = 0 {y~^^‘^ \ny) for y —>■ oo. 
For the required derivative of the analytically continued 
logarithm of the dost function one finds 


dlnF(iyy) ^ ^ ^ 


dy 


2y 2(4 -b TT^y) 


ln(7r^r//4) tt^ 

Vv ^ 


(EIO) 

Substituting this, one can compute the integral in a 
straightforward manner using Mathematica [S2]. One 
obtains 

/.oo /.oo Q 

J dyliy]x) = -j dyy—I{y;x) 

= 27r-2[l _ 21n(2/7r)], (Ell) 

which yields 

Tj^\[v4-,L,m) = ^ (^-l + 21n^^ +0(1/L). 

(E12) 

The right-hand sides of Eqs. (|E4|and ( |E12[ ) can now 

/e at 

2|m|L 


be added and inserted into Eq. (|7.50a|) to arrive at 

-.2 


r ^ T \ -T- /r 1 r \ 1 ■ 

/^(L,m) = J-^([u*];L,m) = —In- ^ 


^-£^ln(^L) + 0(l/L). (E13) 
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tical with Eq. (7.581 when m < 0. 


As a consequence of Eq. (6.18), the result (E13) is iden- tional l,x), 


At the end of Subsection VIIB we had announced a 


check of the asymptotic behavior of the scaling function 
Y(x) for X —>■ —00 by an independent calculation. The 
foregoing result provides such a check. It enables us to 
prove that the coefficient Ci has the analytic value given 
in Eq. (7.63). In fact, straightforward differentiation of 
Eq. ( |E13p gives 

df^{L,m) 

dm 

m m, 2|m|L mvn mln(uL) , , 

= ^ ^ In ^-^ + 0(1/L) 

TT'^ TT 47r TT'^ 


L'JZ{[v^]\L^m) mvQ mln(//L) 


Stt 


47r 


0(1/L), 


(E14) 


where the second equality follows via Eq. ( |E5[ ). S etting 
/r = 1 in the latter and comparing it with Eq. (7.62) 
yields the analytic value of Ci given in Eq. (7.63). 


To obtain the result (7.58) for m > 0, one must sim¬ 
ply add the bulk contribution oc 6 (m) along with the 
term m?6{m) involving the universal amplitude differ¬ 
ence. Note that this to > 0 result means that the scaling 
function Y (x) for x > 0 can be written as 


+ &{x), 


7C(3) 3-21n(2x/7r) 

TT® 


J+ 

87r 


(E15) 


where 0(x) vanishes ~ e for x —>■ oo [cf. Eq. (7.66)]. 


We checked that the two terms in the first line of 


Eq. (E15l can be derived directly from the free-energy 
functional (7.50b) by showing that Y{x) agrees asymp¬ 


totically for X —^ - 1-00 with these terms to 0(x). In this 
calculation, we used the representation of the logarithm 
lnF_|.(k) of the dost function given in Eq. (4.55). As a 
by-product, we obtained the following alternative repre¬ 
sentation of the integral defined in Eq. (4.51): 


= - 


4 ln(7r/2) 


dt 1 


TT I -I- t 


I 


arctan-\/t 7r_ 


(EI6) 

We also verified by numerical integrations that the inte¬ 
grals on the right-hand side of this equation and in the 
original definition (4.51) yield the same values up to an 
error < 10“^^, even though we did not work out an ana¬ 
lytical proof of their equality. 


Appendix F: Derivation of the relation 17.621 


To prove the relati on (|7.62 ) from which Eq. (7.611 fol¬ 
lows by means of Eq. ( 7.54[ ), we set fj, = L = 1 and to = x, 
and consider the derivative d/dx of the free-energy func- 


T- /r 1 1 ^ xn{[v^];l,x) uo 

= -g;;-+ ^ 

f°° dw r , , ^ 2xlnw 

/ — x4a;;I,x) -b Vca-b 

Jo 87r L 


(FI) 


7^([u*];I,x) 2x^ 

2yAJ 71^(1 -I- w) - 

The latter has an explicit dependence on x and an 
implicit one through the x-dependence of the maximizing 
potential u»(z) = I,x). Hence 

d (5J'i([u*];I,x) au*(z;I,x) 

97~i([u«];I,x) 

dx 

Since [9z;(l®/9x]dx is a variation 5v that does not mod¬ 
ify the singular part of u* and we have shown that 
I, x) = 0 for such variations, the nonsingular 
part u"® of u* does not contribute to the first term on 
the right-hand side. The only contribution to this term 
originates from the x-dependent singular part 4x7r“^2:“^ 
of Moreover, we can decompose this part into 

a sum of a component that is singular at z = 0, 
namely 4x7r“^z“^0(zo ~ z), and the nonsingular one 
4 x7r“^z“^0(z — Zq), where Zq G (0 ,1/2) can be taken 
to be arbitrarily small. Since the latter corresponds to 
the derivative of a nonsingular variation Ju, it also does 
not contribute to the first term on the right-hand side of 
Eq. (F2). This decompostion of the term ex z“^ of u must 


be consistently made also in other u-dependent terms 
such as 77.([u*]; I,x) to determine those of their contribu¬ 
tions to dJ'i([u*]; I, x)/dx that must be kept. In the case 
of 7?.([u*]; I, x), we may replace u* by 4x7r“^z“^0(zo — z). 
In conjunction with the subtracted I/z term in u*®, a con¬ 
tribution 4x7r“^z“^0(z—zo) results whose implied contri¬ 
bution to d77.([u*]; I, x)/dx must be retained. One thus 
obtains 

d 7^* x(7e + ln4) 

= ^ +-^^5- 

dz I 

TT'^ ./^o n->-oo 


X 

As 


where 7?.* = TZ{[v^]; I,x) while Iq denotes the integral 

r2, 


Iq= dw 
Jo 


>0 


Svt(z) TT^Z 


+ 


4x 


, , r^^dz\ 21na; 

7e + ln4 - / — M- 


z y TT^-yo; 


77^(1 -I- cj) 


(F4) 


one 


with X* = x([u*];a;). From Eqs. (7.38) and ( 7.4 
sees that the functional derivative i5x*/(5u*(z) yields the 
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Green’s function Q{z,z-,—u}), i.e., 


We can now compute the contribution of the first term 


6 v^{z) 


= g{z,z;-uj;l,x) =-{z\{ui+ 'Hv,) \z). (F5) 


Since Zq is arbitrarily small, we need the behavior of this 
function asymptotically close to the boundary plane z = 
0. For such values of z, we can safely replace the Green’s 
function ^(z, z; —w; 1, a;) by its semi-infinite (L = oo) 
analog g(z, z; —uj; oo, x) and omit o(z) contributions to 
u*. This gives 

^|z)l,=oo 

=-(p{z,p,K) fj,{z,p,K), (F6) 


where p = y/uj + vqx^, xr = 4x/tt'^, and (p{z,p,>c) and 
p{z,p,k) denote the two linearly independent solutions 
of the differential equation 

- dli^iz) + + 7 +P') V’(^) = 0 (F7) 

on the half-line 0 < z < oo that are fixed by the boundary 
conditions 

(p{z 0+,p, xr) = yi [1 -b 0{z)], 

p,{z 0-|-,p, xr) = —v/z [Inz -b 0(1)], 

/x(z —>■ oo,p, xr) —)■ 0. (F8) 

Explicitly, these two solutions read 

p[z,p, >c) = -^= M_^/(2p),o(2pz), (F9) 

Y 

p.{z,p, k) = W^-x/(2p),o(2p^), (FIO) 

where M^^y{x) and Wk„u{x) denote Whittaker functions 
[55H55] . We expand these solutions at fixed pz to first 
order in x’, 

ip{z,p,K) = [lo{pz) -b x:zhi{pz) -b O(xr^z^)] , 

p{z,p, x) = \fz [Kq{pz) -b x-z h 2 {pz) + O(x-^z^)] , 

(Fll) 


and introduce the two functions 

hi{z) = -(2z)“^/^5„Mk,o(2z) 

and 

7e + In 4 


K = 0’ 


(F12) 


h 2 {z) = -Ko{z)- 


2z 


-(2z)-3/2V^a„W«,o(2z) 


K,—0 


(F13) 


Substitution of the expansions (Fll) into ( |F6[ ) then yields 

~ -zIo{pz)Ko{pz) - ^z^[Ko{pz) hi{pz) 
ov»(z) 2 <i nz 


■ loipz) h2{pz)] 


on the right-hand side of Eq. (F14) to the first term on 
the right-hand side of Eq. (F4) for XI ^ oo and small 
zo > 0: 

8 dz 

- 2 / —z loipz)Ko{pz) 

Jo Jo z: 

= - 7 e - ln(8zo) - i InO] -b 0 {XI~^^'^) 

TT Z 


do; 


dz^ 

2 Inw 

L 

7r2^w_ 


0(0 


-1/2n 


(F15) 


The result shows that this contribution is canceled by the 
sum of terms in the second line of Eq. (F4). It follows 
that 


lim 7 q = —— lim 

>-oo TT"^ Q—>-oo ^ 

rzo 


dw 


1 

1 -b W 


8 f ° 

-b ^ / dz z[Koipz)hiipz) + Ioipz)h 2 ipz)] 
'x Jo 

—dr ( r^oVd ip 

lim \ / du ^[Koiu)hi{u) 

TT^ n->oo ( 7 o TT^ 

+ Ioiu)h2iu)] In (zo's/O/it) -b InO > -b O(zo). 


(F16) 

One of the required integrals can be computed exactly, 
namely 

r°° TT^ 

/ duu[Ko{u)hi{u) + Io{u)h 2 {u)] = . (F17) 

Jo ^ 

One way to obtain this result is to express the integral in 
terms of the hypergeometric function 2 F 3 and the Meijer 
G-function G^’^, thereby showing that it agrees with the 
negative of the integrand of the integral denoted Ho{ 0 ) 
in [35], whose value tt ^/8 was compute d th ere. 

Inserting the above results into Eq. (F3) finally yields 


the desired relation (7.62) with Gi given by Eq. (7.64). 


(F14) 


Appendix G: Asymptotic behavior of jz^\\v\t\l,x) 
for X —00 and computation of d\ 

The results for the free-energy functional J^i([u]*; 1, a;) 
derived in Sec. jVII B| and Appendix [^ yield the contri¬ 
butions to the scaling function Yix) that do not van¬ 
ish as I a; I —)■ 00 . Upon exploiting the relationship 
of J^l^l([u*];l , x) to Y'ix) along with the asymptotic 
X —)■ —00 form ( |7.71 ), we have also been able to obtain 
the asymptotic behavior of J^l^l([z;*]; 1, x) for x —>■ —00 
to order l/|x|, expressed in terms of the still unknown co¬ 
efficient di- What remains to be done is the calculation 
of 1, x) for x —?> —00 to 0(l/|x|). 

The contributions of order |x|“^ln|x| and l/|x| arise 
from two sources. First, the self-consistent potential 
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|x|, — 1 ) differs from the asymptotic one oo, — 1 ) 
by a distant-wall correction. As will be become clear be¬ 
low, the latter affects the dost function in a nontrivial 
fashion. Second, the limiting x —>■ —oo form given in 
Eq. (E9| cannot be used for values y < \x\~^ because 
contributions from the second dost solution ^ e“^^ are 
not negligible. 

In order to separate these two types of corrections, 
we split the integral dw on the right-hand side of 

Eq. (7.50c I as f^doj + dw, where is an arbitrary 
number such that 


1 < x‘^{27T/bf < < 


(Gl) 


( 2^ 

Thus J"! 1, a:) becomes 


+ J<,2(a;;fd) -i- J^{x-fi), 
(G2) 

where the first two integrals are given by 


J<Ax-,n) = - j — [xiuj]l,x) + yA] 


(G3) 


and 


j<AxA) = - 

Jo 


dw 


Stt 
2a;2 
7r^(l -I- w) 


77.*(l,a:) 

2y/iJ 


2x Inw 
2 


TT 


/W 


(G4) 


respectively. The third integral, J>(a;;n), stands for 
— ■ of the sum of the integrands of J<p(a;; Id) 

and J<, 2 ix; ft)- 

The calcu lation of J <^2 to the required order is easy. 
Using Eq. (7.72), one can conclude that it is suffi¬ 
cient to substitute the contributions ~ \x\ for 77.*(l,x) 
since the remaining terms yield contributions of order 
0{x~‘^ In |a:|). The integration can be done in a straight¬ 
forward fashion, giving 


7 7 0\ W1 , ON , |a;|V^, 

J<Ax\ fd) = —^ ln(l -I- Id) H-^ In ■ 

47r'^ TT'^ 

-I- 0{x~^ In |x|). 


2\x\ 


(G5) 


We next turn to the calculation of J<p(a;;fd). To de¬ 
termine its X —>■ —oo asymptotics, we need the regular 
solution l,x) of the differential equation 

i-d'^+ v{i;l,x)+u;)(p{i,iA.x) =0, w = (G 6 ) 

both in the inner region 6 /|a;| < 3 < 1/2 and the bound¬ 
ary region 0 < 3 ^ b/\x\ (cf. Sec. VIIA). 

In th e inn er region, w, is small ~ l/|a;| and according 
to Eq. (7.16), can be written as 




AAA 

4\x\ 


0{\xA) 


(G7) 


where w{i) is defined as 

w{i) = lnr( 3 ) - lnr(l - 3 ) - V'(1/2)(23 - !)• (G 8 ) 


We treat u* as a perturbation and subtract from the in¬ 
tegrand of the resulting integral its Taylor expansion to 
0 ( 3 '^) about 3 ' = 0 to ensure its convergence at the lower 
integration limit 3 ' = 0. We thus obtain 


cosh(3l) 1 


|a;|i/2 

ri 


4|a;|3/2 


-b 




/o 


' 1/2 
sin [(3 - 


di'w"'A)P2{A2,i) 

3')I] cosh(3T) 
i 


-P2{Ah^) 


■ 4Ai cosh( 3 t) 


-b 4i?i sinh( 3 l) I-b o(|a;| 


(G9) 


with 


P 2 {AdA) = i ^sinh( 3 t) - 3 'cosh( 3 l) -b 13 '^ sinh( 3 t), 

(GIO) 

where the coefficients Ai and Bi still need to be deter¬ 
mined. 


Before doing this, we first integrate by parts to rewrite 
Eq. (G9) as 


AhiA,x) = + ^^|u;'(3)cosh(3t) 

— 2w(^) sinh(3l) t 

+ Auj f d 3 ' wA) cosh [(3 - 23 ')?] 

Jo 

+ AAi cosh( 3 l) -b 4i3i sinh( 3 t)| 

+ o{\xAA- (Gil) 


This can be inserted into Eq. (7.43) to compute 
A(— gj; 1, x ). From the latter, x(a;;l,x) follows via 
Eq. ( 7.39| ). We thus arrive at 


xA; l,a;) 

^ 1 ’’■1 , - 2 (\ „ Ai-b i?i coth t 

= -e - In - In (I - e - 2 

2\x\ I ’ l^l 

2 g2 l•^/2 ^ 

--p-r / d 3 'r(;( 3 ')[e“^^ *-b (I - cothI)sinh( 23 'I)] 

\x\ Jo 


o{u 


|- 2 \ 


(G12) 


In order to determine Ai and i?i, we compute the reg¬ 
ular solution in the boundary region and match the result 
with the one for the inner region given in Eq. (Gil) at 
the matching point (mp) 3 nip = h/\x\. For 3 « z^p <C I, 
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Eq. ( |G11[ ) simplies to 

l,a;) 

tV _ J_ 

2 4|x|3 


= 1 + 




+ 


- 7 e - ln( 43 ) 


+ ^[(^1 + In2)fi2 - C(3)] + 0(3^) + 0{x-^). 

(G13) 

In the boundary region, the potential x* cannot be 
treated by perturbation theory. To determine the regu¬ 
lar solution in this region, we use the scaling properties 
r;*(z; |x|, - 1 ) = v*{y, l,x)/x^ and 

(/3(z, k = e/|x|;|x|,-l) = ^/\x\ip{i,l-,l,x), (G14) 


where z = |x| 3 , to draw from Eqs. dTTTl ), dT^ , and 
(G7| the following conclusions: 


1 . the potential z;*(z; \x\, — 1 ) can be expanded about 
X = —oo as 


i;*(z; |x|, - 1 ) = ii*(z; oo, - 1 ) -|- 


U3(z) 


■ 0 ( 1 - 


|-4\ 


where M 3 (z) behaves asymptotically as 

W3(0) = 0, 
z—^0 

U 3 (z) = -C(3) -/ 32 Z~^ -kO(l/z2), 

Z—>-oo 


(G15) 

(G16) 

(G17) 


by consistency with Eqs. (3.36) and (3.37). 


Since this remains to be proven, we will not make 
use of it here, taking P 2 to be an unknown number. 

2 . the function ip{z, k; |x|, — 1 ) is the solution to 

[-5^ -kx,(z;oo,-l)](^(z, k; |x|,-l) 


M3(z) 


+ 0(x-") 


(/ 7 (z, k; |x|,-l) (G18) 


subject to the boundary condition (6.3). 


The contributions to the potential describing the de¬ 
viation from u*(z; 00 , — 1 ), which are given on the right- 
hand side, can be treated perturbatively because k^ = 
0{x~‘^). The term U 3 {z)/\x\^ can be replaced by its value 
at z = 00 , namely —C(3)/|xp. The reason is that the dif¬ 
ference U 3 (z) — ^ 3 ( 00 ) should be expandable in inverse 
powers of \x\ (modulo eventual logarithms) and hence 
vary for z —>■ 00 as z~^ with k > 1. In fact, k is ex¬ 
actly unity since we shall show below that the potential 
coefficient /?2 is given by Eq. (3.38) and hence does not 
vanish. 

Once rt 3 (z) has been replaced by its limiting value, it 
can be absorbed by introducing the shifted wavenumber 
kx through 


k^^k^-C(3)/N' 


(G19) 


To determine the coefficients Ai and Bi, we will match 
the pertur bative result for :p(z, k; |x|,—1) obtained from 
Eq. (G18) at the right edge of the boundary region, 
namely for z=|x|3«53>l, with the result given in 
Eq. (G13). For such values of z the term cx u^^z) is 
smaller than the second one by a factor of It does 
not contribute at the order of the terms involving the 
coefficients Ai and Bi, vanishes in the limit 5 —^ 00 , and 
hence may be ignored here. 


The associated perturbative solution of Eq. (GI 8 I can 
be read off from previous ones derived in Sec. |VI C[ Not¬ 
ing that the quantity k^ in Eq. (6.30) corresp onds here 
to k^ — C(3)|x|“^, we can conclude from Eqs. (6.25) and 


(6.35) that the regular solution at z = |x |3 
be written as 


5^1 can 


(p(z, k; |x|,-l) 




k2_^ 


z / 3 4z 


(G20) 


Matching this with Eq. (G13) at z = |x |3 then shows that 
the coefficients Ai and Bi are given by 

Ai = — In 2, 




■ln(|x|/7r)]Vw. 


(G21) 

We can now insert Eq. ( G12[ ) into Eq. ( |G3[ ). For those 
parts of the integrand of the integral dw that are ex¬ 
ponentially decaying, such as the ones proportional to 
In (1 — 1 — coth-^/w, and we can safely 

replace the upper integration limit by 00 . Using 

poo 

-C(3) = 2/ dw In (1 - e-^'^) 

Jo 

roo 

(1 — coth -x/o;); 


(G22) 


one arrives at 
>^<, 1 ( 3 :; U) = 


- :r+ln 




47r|x 


In 2 — 


2|x 

C(3) 


U3/2 

127r|x| 


1-kln^ 

TT 


87r 


1 1 -Hn(|x|/7r) 

2 |x| 


+ >7<,3(x; U) -|- T<^ 4 (x) -)- 0{x 


(G23) 


with 


1 /-i/z 

J< 3 (x;U) = -—r / dujuj dz w(z) 
47r|x| Jo Jo 


47r|x 

1 

27r|x 


. 1/2 


dzw{z) R{z\ U) 


(G24) 


and 


I ri/2 roo 

{x) = —:—r / dzw(z) / du) uj (l 

^Ax\Jo Jo ^ 

— coth(-\/a;)] sinh{2z^/u}) 

_ di,i 
87r|x| 


J<A 


(G25) 
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Here 

and 


3 _ g-2zni/=(3 g^f|i/2 g^2^ 4 z3h3/2) 


8z4 


(G26) 


o ^1/2 

rfi,i = -J^ dz w{z) [z-4 - (z)/6], (G27) 

where i(;(z) denotes the function 

w{z) = Inr(z) — lnr(l — z) — (2z — l)ip(l/2), (G28) 

which varies as 

w(z) = Wi,s(z) + O(z^) (G29) 

z-fO 

with 


Wa.s(z) = - 7 e - ln(4z) + 4z In 2 - -((3)z^. (G30) 

To further evaluate J<^ 3 (a;;H), we divide up w(z) 
into its a symptotic contribution Was ( 2 ) introduced in 
Eq. (G301 and a remainder, writing 


J<, 3 (a^; n) = n) + 4%(x; H)] (G31) 


( 2 ) 


87r 


with 


4 

=-—r / dz Wa,s(z) E(z; fi) (G32) 
fI Jo 


and 


4 

J^l(x;fi) = --r—r / dzlw(z) - Wa.s(z)]B(z;Q). 
fI Jo 

(G33) 

In the last integral, we can replace R{z] H) by its 
asymptotic large-z form 3 /(8^4) to obtain 


H) = |a;r 1 [di ,2 + , (G34) 


with 


3 

di ,2 = -2 J dzz ‘^[w{z) - Wi,s{z)]. (G35) 

The large-il asymptotics of the remaining integral 
can be determined in a straightforward fash¬ 
ion. One finds 


J^%(x; fi) = 


l + ^\n {2/Vn 


- 2 H In 2 


- ^ -47E + 81n2-bC(3)^7E- y 
+ In + O . (G36) 


The combination of Eqs. (G23)-(G36) then yields 




1 J C(3) [S 2^^7rH4/2 




1 


FI 

7rHi/2 


2 ^""“21x1 


- 3 - 47e + 8 In 2 


In 




9 3 

1 

\x 

|a:| 

„-2 


2|i 


di,i + di 2 


7e 


+ 0{x 4ln|a:|). (G37) 


It remains to calculate the x ^ —00 asymptotics of 


integral J>(a;; H) in Eq. (G2) to the order ~ |a;| Upon 


making a change of variable uj ^ y = x ‘^oj, it becomes 


J>(x;H) = 


J Qjx‘ 
TZ*{1, x) 
‘^\x\y/y 


dy 

Stt 


x{y\ kl.-i) + \x\^ 

2\n{yx^) 2 




\y- 


P-2'1 


(G38) 


where 77.*(l,a;) can be read from Eq. ( |7.72 ). The inte¬ 
grand in the square brackets is the (rescaled) analog of 
the one we encountered in Appendix]^ in our calculation 
of F^\[v^]\L,m) to o(L). There we could express the 
asymptotic from of the function x{yxn^] L, m) for large L 
in terms of the Jost function F{iyJy) for the semi-infinite 
case L = 00 . Here we must take into account the effects 
of the finite film thickness. The differential equation from 
which the reg ular s olution (p(z, k; |a;|, — 1 ) is to be deter¬ 
mined is Eq. (G18). Since we shall need the value of the 


potential coefficient P 2 introduced in Eq. ( G17l, we begin 
by proving our result stated in Eq. (3.38). 


To this end, we shall compute the contributions of 
72.* (1, a;) prop ortion al to x~^ mod In |a;| directly from 
the definition ( 7.46 ) of 7^([u*]; 1, a:) and match them with 
those in Eq. (7.72). These terms originate from the con¬ 


tribution u)!®. Setting L = 1 and m = x, we therefore 
split the corresponding integral into a sum of contribu¬ 
tions from the boundary and inner region. Making the 
change of integration variables 3 —>■ z = \x\i in the near¬ 
boundary integral, we arrive at 


uj® = |a;| / dz 
Jo 
/■ 1/2 

+ L d3 

''b/\x\ 


, , , . 1 4 

4 '*(z;|a;|,-i) + —y-b-y- 


4z2 

, , . 1 4|a:| 

>■(3; 1) x) + -I ^ 
432 7r23 


. (G39) 


In the first integral, we insert the large-a: expan¬ 
sion (G15) of u*(z; |x|,—1). For the potential in the sec¬ 


ond integral we can substitute the corresponding expan¬ 
sion 


u*(3; l,x) = 


G(3) , hii) 


(G40) 
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where 


^ 1 ( 3 ) = 


V’"(3) + V'"(l-3) 


(G41) 


for b/\x\ < 3 < 1/2 according to Eqs 

( 

7.16 

) and 

Q 

00 

Furthermore, consistency with Eqs. ( 

3.36 

), (3.37) 

, and 


(G17) requires that 

^2(3 —>■ 0 ) = —Q !2 3 ^ ~ 102 3 


+ 0(1). (G42) 


The contributions linear in |a;| or proportional to 
|a;|ln|a;| of the first term in Eq. (G39l originate from 
the leading x = —00 term of the integrand. Their b- 
dependent part is canceled by analogous ones produced 
by the antiderivatives of the last two terms of the second 
integral’s integrand at the lower integration limit. All 
contributions linear in |a;| and proportional to |a;|ln|a;| 
must add up to produce t he co ntribution to TZ^,{l,x) 
given in the first line of Eq. ( |7.72 ). Since the latter equa¬ 
tion is already established, there is no need to check this 
contribution again. 

The contributions to 77.*(1, x) of order x~‘^ mod In |a:| 
arise from the sum of the integrals 

*1/2 

Ti = 2\x\-^ / d3ti(3) = {2\x\)-^w"{};)\y^ 

Jb/\x\ '' ' 

^ N , 2&C(3) 

2b 


+ 0(|c 


|-4\ 


T, = 2 |a; 


-2 


T. = 2|cc| 


*1/2 

/ djtzCa), 

lb/\x\ 

/ dzM 3 (z). 


Assuming that 6 > 1 and defining 
^ I'T'i I “2 /32 

) + pi + 7 j 


L2{i) = / dj' 

L 3 (z) = r dz' [U 3 (z') + C(3) + h eiz' - l)/z '], (G44) 
Jo 

we can rewrite T 2 + as 

T2+r3 = ^|L3(^)+A2(V|x|)-/32ln^ 

+ y(8-S)-k(3)} (G45) 

by adding and subtracting appropriate terms to their in¬ 
tegrands. We can now let 6 —>■ 00 and x —>■ — 00 such 
that b/\x\ —>■ 0 , i.e., we set b = const |x|'‘ with 0 < k < 1 
and let x —>■ — 00 . Since the contributions oc of Ti 

and T 2 + T 3 cancel each other, no contribution ^ 
to 77.* results and we obtain 


7^*(l,:r) = -^ 


r , 21x11 

2 

l-ln-^ 

H— 2 

TT 



- 132 ln(|x|/ 2 ) + -02 


7.2(0) + 7.3(00) 
o{l/\x\^). (G46) 


Gomparison of this result with (7.721 yields, besides the 


value of (32 given in Eq. (3.38), the following relation 


C( 3 )^—\^= 7.2(0)++3(00) + ^ 02. (G 47 ) 

Until now we have not taken into account the far¬ 


boundary correction oc M 3 (z) of Eq. (G15) other than 
through its asymptotic value 113 ( 00 ). In our calculation 
of T>(x; n) we must be more care ful an d include also its 
correction oc z“^ specified in Eq. ( |G17 ). 

To this end, we express k in the Schrodinger equa¬ 
tion (G18) again in terms of the shifted wavenumber k^. 


defined in Eq. (G19). As a result, the potential term 


M 3 (z) gets replaced by the shifted one 

Am3(z) = ■u3(z) + C(3) ~ -/32z~\ (G48) 


which vanishes at z = 00 . Let us introduce the notations 


/+(z, k^-,x) = 

F+{K) = FiiK), 


(G49) 


for the Jost solution and dost function on the right-hand 
sides. Exploiting Eq. (4.11 ), we can express the regular 
solution introduced in Eq. (G14) as 


(^(z, k; |x|,-l) 

^ F+ik^)f+iz, K]x) - F+{-k^) f+{z,-k^;x) 

2k„ 


(G 50 ) 


(G43) Here /+(z; |a;|, — 1) is a solution to 


+ v{z; 00 ,- 1 ) + + kl f+{z,k^; x) = 0 

(G51) 

(G52) 


that behaves asymptotically as 

f+{AK;x) ~ exp k.*z- 

z— >-oo 


132 


Inz 


2ka: 

where the Inz in the exponent is due to the slow 
Goulomb-like decay of Au 3 (z) m- 

To determine /+ (z, k; x) in the boundary region for 
1 < C z < b, we replace u*(z; |a:|, — 1 ) and Au 3 (z) in 
Eq. (|G51 ) b y their large-z expansions — (2z)“^ — a 2 Z“^ + 


0(z ’”) and (|G48|, respectively, and use the ansatz 


f+{z,k,,-x) = A(ka;,a:)exp |k 2 ,z 

132 


+ 0(z-^) 
Solving for / 2 ,..., /4 yields 
/2(kx,x) = 


2 k., 


i+E” 

7=2 








4k3|x|3 8k4|a;|' 


(G53) 


(G54) 
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f3{K,x) = 


and 


1 


132 


+ 


PI 


Sk? Sk^lxP 8k^|a;|® 32k®| 


(G55) 


f4{kx,x) = 


8k3 


, r, 1 + Q(l/ka;) ^ 

6k? 24ki|x|3 


1-6 


)• 

(G56) 

We now express the exponent of the exponential of 
/+(z, kx,x) in terms of k, expand in powers of l/|x| and 
1/k, dropping the terms of higher than the order s k“^, 

I a; I “3^ and z“3^ and insert the result into Eq. (G50). The 


contribution from /_|_(z, — k^;,a;) can be omitted because 
it is exponentially small. We thus obtain 


(^(z, k; \x\ 


1 -j ^+(k3;,a:) 

/32 In z 1 


1 


8k z2 

P2 


Oi2 

6k z3 


2k|a;|3 

zC(3) 

2k|a;|3 


8k2z3 4k2|a:|3z 


(G57) 


Having determined the regular solution in the bound¬ 
ary region, we now turn to its calculation in the inner re¬ 
gion. Since we need the asymptotic behavior of x(w; 1, x) 
for large oj = , we can again use the semiclassical ex¬ 
pansion to determine the behavior of l,x) in the 

inner regime. We start from the analogs of Eqs. (D3|- 

dD^, 


r 1 

<f( 3 ,t;l,x) =A(i,x) expj^t 3 -—J dj' 


^ 1 ( 3 ') 


^ 2 ( 3 ') 

1^12 


1 

'^ 1 ( 3 ) ^ 2 ( 3 )' 

412 

[ |x| |x|2j 


0(i 


- 3 \ 


(G58) 


perform the integral involving ti(y) using Eq. (G411, 

^xpres; 

C(3)3 


make a Laurent expansion in 3, and express the integral 
of ^ 2 ( 3 ^ i'l terms of ^ 2 ( 3 )- This gives 


<p(3,t;i,x) _ L 1 
A{i,x) 


1 


8|x|32 
-^ 2 ( 3 ) +/32ln(23) 


1 


2|x|2 

1 

[8|x|33 
02 


m 

4\x\ 


4|x|234 


0(3^ 


2|x| 

02 

6|x|2 

P 2 

4|x|23 
1 ,- 3 ' 


0(3^) 


+ 0(r3)^. (G59) 


Matching this with Eq. (G57l at 3mp = Zmp|x| = h/\x\ 
then yields for the amplitude the result 




21 

4«|x|2 2 


C(3) 


4o2 


4fi2|x 
£ 2 ( 0 ) 

2«|x|2 


3«|x|2 
+ o(x- 2 ) 1 .(G 60 ) 


where l/k^, = k~3 exp[^(3)/(2fi^|x|3) -|- 0(|x|“®)] was 
used. 

The result must be inserted into the analog of 
Eq. (Dll), namely 


x({^; 1, x) = —t — In [27rlx)^]. 


(G61) 


To this end, we expand lnE'+(ka;,x) in its second argu¬ 
ment about In El), (k, 00 ) and subsequently its first argu¬ 
ment to linear order in k^^ — k = 0(|x|“3)j obtaining 


lnE+(ka;,x) = lnE+(k:j,,oo) -I- 


Xsikx) 


+ 0(|a 


|- 6 \ 


= In E+(k, 00 ) — 


C(3) i91nE4.(k, 00 ) 


2k|x| 


elk 




(G62) 


where the term oc X 3 {kx) accounts for the change oc 
|x|“3 induced by the far-boundary correction Am 3 (z) to 
u*(z;00,-1) [cf. Eq. (G15)]. 

We can now substitute these results into J>(x, fl), 
change to the integration variable y = ixijx^^ and use 
the fact that the contributions to the integrand resulting 
from the first term inside the curly brackets in Eq. (G60) 


and the term involving clkE+(k, 00 ) add up to a derivative 
dy of a function that vanishes at y = 00 . It follows that 
J>(x, H) can be decomposed as 


J>(x, D) = J>_i(x, fl/x^, 00 ) -|- J>_ 2 (x, D) -I- 


where J>p(x,yi,y 2 ) denotes the integral 

J>,i{x,yi^y2) 


Q(H/x^) 

|x| 

(G63) 


2 py2 f 

= dy|21nF+(v^,oo)-ln 


2Vy 


’■Vy 


1 -hln 


^Vy 


+ 


7r^(y+ |x|-2) 


(G64) 


while 


J>,2(x,H) = - 


C(3) 

47r|x| 




d!,|-In T 


dy 


C(3) E+(y3/2,oo) 


.1/4 


J/ = n/|£c |2 


(G65) 


and 


47r|x| (2/7r)3/2yi/4 

Q(yo) = ^ / dy[2X3(^)-L3(oo)y"^/2]. (G66) 

'^yo 

The integral J> i(x, D/x^, 00 ) can be decomposed into 
the ^-independent term J>_i(x,0,oo) and the remain¬ 
der — J>.i(x, 0, D/ x^). The value of the former can be 
gleaned from Eq. dElI] ), 


J>,i(x,0,oo) = — 


-1 -H21n 


21x1 


(G67) 
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Since the upper integration limit of the latter integral 
J> 4 (a;, 0, is small, n/a;^ <C 1, we can replace the 

logarithm of the Jost function in its integrand by the 
asymptotic small-y form 


In F+{y/^, oo) =\n^/y- ^ 


1-ln 




0{y\ny), 
(G68) 


which follows by integrating the analog of Eq. (ElO) for 


E+(y^, oo) with respect to y and expanding the result 
for small y. One obtains 


3 


TT^ 

n 

ri 

, 7rHi/2] 



— In —^—— 

87r 

2 

2|a:| . 

H3/2 

1 1 ^ ttI 

jxjTT 

9 12 ^ 2 


101/2 

In 


01/2 


2x 


+ 0(^x ^ln|a;|), (G69) 

and substitution of Eq. ( |G68[ ) into Eq. (G651 yields 


J>,2(a;, O) = 


C(3) TT^O 


167r|x| 4|x 


+ 0(a:-2ln|a:|). (G70) 


Eor large negative x, the last term on the right-hand 
side of (G63) behaves as 


Q(0/a;2) 

kl 




)• 


(G71) 


Hence, to the required order in l/|a;|, all 0-dependent 
terms of J>{x, O) are contained in 0)-|-J> 2 ( 3 ;, O). 

Substitution of our results for this sum along with <5(0) 
and those for J<_i(a;,0) and J<,2(a;,0) into Eq. ( |G2[ ) 
yields 


.FP'(Wa,»)=|p(-l + 21 n?M 


Stt 


+ 2 In I 


+ o(|a:|-i), (G72) 


with 


di — + ^ 1.2 ~ g ~ 47 e -|- 81n2 


71-2 C(3) 


■7E + in—- (G73) 

As one should expect, all 0-dependent terms cancel in 


(G72). 


Insertion of Eq. (G72) into Eq. (7.65) then gives a dit 


ferential equation for Y (x), which can be solved to deter¬ 
mine the asymptotic large-a: form of this function. The 
result is [cf. (7.59)]. 


y(x) = - 




7C(3) 1 2|x| 

TT^ 2 TT 


0(a 


(G74) 


where 0(a;) has the asymptotic series expansion (7.71). 


To establish the result for di given in Eqs. (7.73), we 
first prove that 


g(o) = 0 , 


(G75) 


if /I 2 and 02 take their self-consistent values, given in 
Eq. (3.38) and unknown, respectively. Subsequently, we 


shall complete the derivation of di by computing the sum 
of integrals appearing in -|- di_ 2 . 

The first step of the proof is to show that 


Lsjoo) 

2k 


-^3(k)= r 

Jo 


dz 


0oo(z,z;-k^) + 


1 


Au3(z), 


(G76) 

where ^oo(z,z'; —k2) = —[k2 +is the Green’s 
function for L = 00 and m = — 1. To do this, we 
can proceed along lines similar to those followed in Ap¬ 
pendix Let /±(z, k) be the Jost solutions associated 
with the Schrodinger equation /± = 0 on the half¬ 


line that vary ~ e*^^ as z —> 00 . By analogy with 
Eqs. ( |B7| ) and (B8), one finds that /-(z —)■ 0, k) ~ 
—E+(k, oo)-yzlnz. Using the analog of Eq. (B2), namely 
(k'2 - k2)(^(z, k)/_(z, k') = one 

can derive 


dz 


lim _ i 

k'-s-k F+(k,oo) 2k 


= d \^2 In 


■E+(k,oo)' 

poo 

- - dz 

[ Tk 

Jo [ 


5oo(z,z;-k2) -K 


2k 
(G77) 


in a straightforward manner, where the last result follows 
both from Eq. ( |B10 ) or else from the fact that the limit 
in squar e bra ckets is nothing other than —5oo(z,z; —k2) 
[cf. Eq. @1. We now replace u* by u* -I- Sv and de¬ 
note the implied change of lnF+(k, 00 ) as i51nF3_(k). 
A calculation analogous to the one in Eq. ( |B12 ) yields 
—d \^2 tr(0oo dv) for the variation of the Green’s function 
term on the right-hand side of Eq. (G77). If 6v{z) de¬ 
cayed faster than as z —>■ 00 , we could integrate with 
respect to k2, taking into account that the integration 
constant must be zero, to conclude that 


poo 

SlnF+{k) = - dz5oo(z,z;-k2)(5u(z). 

Jo 


(G78) 


However, we are interested in the variation 6v{z) = 
Au 3 (z)/|xp which has the slower Coulombic de¬ 
cay (G48). We therefore replace U 3 (z) by ur{z) = 
0{R - z)u3(z) in 6v. The associated regular solution 
varies as ifji{z,k) ^ exp[kz — /32(2k|a:p)“^ Inz] and ^ 
e''^ for 1 ^ z < i? and 1 ^ i? < z, respectively. 
Taking into account that (pn and its z-derivative must 
both be continuous at z = i?, one sees that the cor¬ 
responding Jost function must behave as E+. fl(k, x ) = 
F+(k, a;) exp[—/32(2k|a:p)“^ Ini?]. Using Eq. ( G78| ) for 
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(51iiF+,_r, we get 


(51ni^+(k, x) = lim 

R^oo 


i51nF+,_R(k,x) + 


^2 


2k|x 


In R 


= lim < — 


i?—>-oo 


dz 


Goo{z,z;-k^) + — 


3 

Am3(z) 


+ 


2k|a 


J dz[AM 3 (z )+/32 6 '(z- l)/z]|. (G79) 


The limit i? —>■ oo of the second integral is ^ 3 ( 00 ), and 
the left-ha nd si de is A 3 (k)/|xp according to Eq. (G62). 
Thus Eq. (|G76|) is established. 


To proceed, we multiply the self-consistency Eq. (4.21) 


by Ait 3 (z) and integrate it over the half-line. The result 
enables us to express the right-hand side of Eq. (G77) 
in terms of the square of the semi-bound state ipo{z). It 
conjunction with Eq. (G 66 ), it yields 


Q(0) = ^|^-k^ dzAu 3 (z) [(pg(z) - 1]| (G80) 


with 


nOO 

J 3 = / dz Ausiz) [(filiz) - 1], (G81) 

Jo 

where Au 3 ( 0 ) = ^(3) was used. Hence we must show 
that 


J3 = -C(3)/2, 


(G82) 


By adding and subtracting appropriate terms to the 
integrand and choosing a finite value R of the upper in¬ 
tegration limit whose limit i? —>■ 00 must be considered, 
one can rewrite J 3 as 

^3 = C(3) J^-E3(oo) 

Jq ^^'“3W‘^o(z)|-(G83) 

To evaluate the integral on the right-hand side, we 
can exploit the exponential decrease ( |7.26 1 of the lowest 
eigenvalue Ei{x)/L‘^ in the limit x —>■ —00 [cf. Eq. (7.1)]. 
This means that x )]/fi( 3 , x) = 0(|x|e“l^l) as 

X —>■ —00. Erom Eqs. (4.20) and (4.16), we know that 
the associated eigenfunction fi(z/|x|,x) approaches the 
semi-bound state fo{^) when x —>■ — 00 at fixed z so that 
fi( 3 ,— 00) = (/ 5 o(oo) = 1. The leading large-jxj correc¬ 
tions to the latter limiting value result from the poten¬ 


tial terms given on the right-hand side of Eq. (G40) and 


hence are of the same orders |x|“^ and |x|“^. Allowing 
for a change of the normalization of the eigenfunction, 
we therefore make the ansatz 


fi(3,x)= [l-k 6 i/|x|-kO(|x| 2)]fi(3,x) 

= 1 + 91 ( 3 ) ^ 1 / 2 ( 3 ) 


(G84) 


for the modified eigenfunction fi( 3 ,x). Eor convenience, 
we choose bi such that 

51(3) = [V’( 3 ) +^(1-3) +27e]- (G85) 

This choice implies that 

gi(3) cs - 1-^32 + 0 (/), (G86) 

a-fO 43 2 

ensuring the absence of a term in this expansion. 
Substituting of the ansatz (G84) and the large- 


jxj form (G40) of into the differential equation 
= C/dxle”!"^!), we can solve for g" and g'^ 

to obtain 


51 ( 3 ) = G(3), 

52(3) =^1(3)51(3) +^2 (3)- 


(G87) 

(G 88 ) 


We now integrate Eq. ( |G 88 ) over the inner region 
[ 6 /|x|, l/ 2 j and add and subtract terms to the integrand 
to obtain finite integrals over [0, l/2j. This gives 


9iib/\x\) = [1-02 

with 


363 


-L2(0)-di,3 (G89) 


.1/2 


di.3 = 


dz 


9i ( 3 ) 51 ( 3 )-^-^ 


(G90) 


The latter integral can be written in terms of the pa- 
rameters and di .2 introduced in Eqs. (G27) and 


(G35). Upon expressing the integrand of the combina¬ 
tion (ii ^3 -I- (di^i -I- (ii, 2)/4 in terms of the function w{z) 
and its derivatives, one can determine its antiderivative 
and compute the integral. The result yields the relation 


<^1.3 = 


di.i + di,2 


7e + In 4 - Q -I- 7e -I- In 2^ . 

(G91) 


Turning to the rescaled analog of fi, 
fi(z,x) = fi(z/|x|,x). 


(G92) 


we note that the leading large-jxj correction to its 
X —>■ —00 limit is due to the |x|“3 term of x*(z; |x|,— 1 ) 


in Eq. (G16). We therefore write 
fi(z,x) = <fo{z) + 


(G93) 


and substitute this ansatz into the differential equa¬ 
tion for fi. Taking into account the large-z behavior 
x*(z,oo,x) ~ —(2z3)-i — a 2 Z“^_^plied by Eq. (3.36lj 


and that of U 3 given in Eq. (G17l, one concludes that 
6-1 (z) behaves as 

61 (z) ~ z^-f Oiz -I- 0(ln^ z), (G94) 

z—^oo Z 
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where the absence of a contribution oc z In z requires that with 


/?2 takes its self-consistent value (3.38). 


To determine Oi, we match the limiting behavior of 
the derivative dzfi{b,x) for large b that follows from 
Eqs. (G93l, (6.25), and (G94) with the behavior of 
fi( 3 ,x) at the matching point jmp = b/\x\. One finds 
that the a;-independent term cx b~‘^ and b~^ as well as 
the ones cx 6 jeep of both expressions are consistent and 
that the coefficient of the 6 -independent term is 

given by 


oi — -(1 — 802) — ^2(0) — di_ 3 . 


(G95) 


We are now ready to express the integral in Eq. (G83) 
in terms of this coefficient. To this end we equate 
the terms of 'H„.fi(z, x) to zero. This yields 

u^{z) V3o(z) = 6,1 (z). Multiplying the result by 

integrating by parts, one is led to 

[ dzu 3 (z) <pI{z) = [ipo{z) h'^{z) - c/?^(z) /Ii(z)]^ 

■3 


= C(3) 


8 


- R 


-|- ai -|- 0(1/i?), 

(G96) 


where we have again used Eqs. (6.25) and (G94). Note 


that the contribution from the lower integration limit z = 
0 vanishes because <po(z) ~ and hi{z) ~ z®/^ as z —>• 0 , 
respectively. To see the latter, note that U 3 (z —0) = 
0 ( 1 ) and hence does not contribute to the pole part of 
u*. However, inserting the ansatz fi(z,x) = y/z[\ + riz + 
Tiz?' -)-...] into the differential equation for fi, one finds 
that the regular parts of f* do not effect ri, which implies 
the stated small-z behavior of hi. 

Gombining Eqs. (G83), ( G95| ) and (G96l gives 


J 3 = C(3) 


= C(3) 


Jip + 


- ^2(0) - ^3(00) - rfl .3 + 


Jv) ]—fxla2 


— di,3 + X 


1 — 802 
3 

(G97) 


To obtain the second line, we eli minate d L 2 (0) -I- T 3 (0) -I- 
802/3 by exploiting the r elation ( G47[). 

We can now use Eqs. (G80) and (|G 8 l|) to eliminate 


(5(0) on the right-hand side of Eq. (G73) in favor of J 3 . 
This gives us di expressed in terms of di .i, di_ 2 j and J 3 . 
Upon inserting the result into Eq. (G971 and using the 


exact result for given in Eq. (6.43) and Eq. (G91), 


we can solve for J 3 . All contributions involving di i and 
di ^2 cance l out. One obtains the result for J 3 given in 


Eq. (G82), which implies that (5(0) = 0. Inserting this 
result into Eq. (G73), we combine the integrals that the 


coefficients di ^2 and fii ^2 involve into a single integral 
and integrate by parts using lf{z)f"'{z) — f'(z)f''(z)]' = 
/(z)/W( 2 ) - [/"(z)] 2 , with fiz) = Inr(z) - lnr(l - z). 
This gives |105j 


<^1 = + 2[7e + ln( 2 / 7 r)] - 


3C(3) 


4C(3) 


(G98) 


/• 1/2 

J 4 = / dz{[tjj'{z)z)f - s{z)}, (G99) 

Jo 

8C(3) 


where 


■sW = 71 - 


(GlOO) 


The square of the difference of "0' functions in the inte¬ 
grand of J 4 can be decomposed into the square of the sum 
of these functions, which reduces to tt^/ sin'^( 7 r 2 ;), and a 
remainder. As a consequence, we arrive at 


J4 = 


r.1/2 


dz 


Lsin‘'(7rz) z4 


3z2 


- 2 J 5 


where 


with 


- (2 - 1 - TT^) — 2 J 5 , 

(GlOl) 

. 1/2 

:/ dz [■0'(z)V)'(l - z) - Si(z)] 

Jo 

(G102) 

, ^ TT^ 2 , . 

(G103) 



It is straightforward to rewrite the last integral as 

-,r2 

(G104) 


J 5 = Je + y+4C(3)ln2, 


where 


Jq = f dz\il}'{z)'il}'{1 — z) — si{z) — si{\ — z)]. (G105) 
Jo 

The integral Jg, in turn, can be represented as the limit 
e —)■ 0 - 1 - of the contour integral 


Jp, = — lim 




e->o+ 27rie 


(G106) 


with 


’/ 7 (e) = j^dz (^-^) i)'(z)ip'{l-z), (G107) 

where the integration contour C runs around the branch 
cut [0,1] of the integrand in the clock-wise direction. Re¬ 
casting the integral Jj^e) into the sum of residua of the 
integrand at the poles located in the exterior of the inte¬ 
gration path C, one obtains 

27ri 

00 

= E 


n—2 

n — 1 
n 


n—1 


?/’"(n) -f e'(/''(n) ( -- 1- 

n n—1 


- ^"(n) -f eV''(n) ( -- 

' 71 71—1 


(G108) 
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After substitution of this result into (GIO 6 I and proceed- ing to the limit e —t 0-I-, one finds 


Je = 2^'{2) + 2Y, 


— 1 + 2 ln(n — 1) 


n—S 


{n — 1)3 


= --2C(3)-4C'(3). 


(G109) 


In deriving the first line of Eq. (G109I we have 


rearranged the summation on the right-hand side of 
Eq. ( G108| ), and used the first and the second derivatives 
of the well-known relation — ip{^ ~ 1) = ~ ^)- 

Subst itution of this result into Eqs. (G104), ( |G101[ ), and 
(G98) finally leads to the result given in Eq. (7.73E 
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